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THE HYDRAULIC JUMP AT AN ABRUPT DROP! 


Walter L. Moore,* A.M., ASCE and Carl W. Morgan,** A.M., ASCE 
(Proc. Paper 1449) 


SYNOPSIS 


The hydraulic jump may form at various locations relative to a low abrupt 
drop in a rectangular channel. The role of the drop in determining the form 
of the jump and in stabilizing its position is clarified by analysis and experi- 
ment. An example illustrates the application of the results to the analysis of 
a stilling basin. 


INTRODUCTION 


In many types of hydraulic structures a considerable portion of the kinetic 
energy of the flowing water must be dissipated to prevent erosion of the chan- 
nel downstream from the structure. Various methods of energy dissipation 
have been used to achieve tranquil flow conditions as the flow enters the down- 
stream channel. The primary purpose of all such methods is to convert as 
much as possible of the kinetic energy of the flow into turbulent energy and 
ultimately into heat. The required energy dissipation is frequently most ef- 
fectively accomplished by means of a hydraulic jump which may assume 
several different forms depending on the geometry of the channel boundaries. 
In an effort to reduce the length of stilling basins and to insure that the jump 
will not be swept from the stilling basins and to insure that the jump will not 
be swept from the stilling basin under varying flow conditions, an abrupt drop 
or rise in the channel bottom may be introduced. 

For economic reasons, there has been a tendency in practical design to 
make the length of the stilling pool considerably less than the length of the 
jump. This paper presents information concerning the characteristics of the 


Note: Discussion open until May 1, 1958. Paper 1449 is part of the copyrighted 
Journal of the Hydraulics Division of the American Society of Civil Engineers ,Vol. 
83, No. HY 6, December, 1957. 
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jump at an abrupt drop in the channel bottom with particular emphasis on the 
characteristics when the drop is located within the length of the jump. 


Previous Study 


Basic information concerning the characteristics of the hydraulic jump on 
a level surface is published in many books on hydraulics and fluid mechanics. 
Mr. Edward A. Elevatorski(1) has compiled a bibliography of published infor- 
mation concerning the hydraulic jump under various channel and flow condi- 
tions. This publication lists some 504 references arranged by years starting 
in 1819 and extending to 1954. 

Information concerning flow problems closely allied with those discussed 
in this article have been published. Moore(2) made an analytical and experi- 
mental study of a free overfall and determined the energy, velocity, and sur- 
face profiles for various conditions of flow. Forster and Skrinde(3) studied 
both theoretically and experimentally the stabilizing effect of an abrupt rise 
at the end of a jump and developed a composite graphical summary of results 
by which the downstream depth after the rise could be computed. Valuable 
contributions to the flow analysis were also made by the several persons dis- 
cussing this paper. Rouse, Bhoota, and Hsu(4) investigated the stabilizing 
characteristics of an abrupt drop in a channel bottom. They found that two 
different types of flow may occur at a given drop, the transition from one to 
the other being characterized by an undular wave. The type of flow which will 
form is dependent upon whether the downstream depth is above or below that 
which produces the undular standing wave. The conditions which produce the 
undular stage (critical transition zone) cannot be foretold by equations utiliz- 
ing only the momentum and continuity relationships but must be obtained by 
experimental measurements. Rouse, Bhoota, and Hsu presented such data 
which verified their analysis and indicated a systematic transition between 
the two regimes of flow depending on the Froude number. To simplify the 
analysis for the two types of jumps, the pressure on the face of the drop was 
assumed to be hydrostatically distributed and assumed to be dependent on the 
downstream depth if the jump is upstream from the drop and dependent on the 
upstream depth if the jump is downstream from the drop. Experimental data 
was not presented for intermediate positions of the jump which were not sub- 


ject to a simple analysis based on the two assumed hydrostatic pressure dis- 
tributions on the face of the drop. 


Analysis 


The different forms of a hydraulic jump which will occur at an abrupt drop 
in a rectangular channel are shown in the definition sketch, Fig. 1. For con- 
venience the different forms are identified as A, W, and B as indicated on the 
Y 
sketch. The relative height of a stationary jump ¥ will be dependent upon 
1 


AZ V 

the dimensionless parameters — » > and the Froude number 
¥ 

1 1 


entering flow. As the jump changes from one form to the other, the value of 


of the 
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y,= Depth upstream from drop 
Yp= Depth downstream from drop 
4zz Height of drop 

Np= Piezometric head on drop 


FIG. |. FORMS OF HYDRAULIC JUMP AT AN ABRUPT DROP. 
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the hydrostatic head, hp, will change causing the pressure on the face of the 
drop to change from a maximum value of yhp = y Y2 when the jump forms up- 
stream from the drop to a minimum value of yhp< y (4Zo + Y1) when the 
jump forms downstream from the drop. Photographs illustrating various 
forms of the jump are shown as Fig. 2. 

If the pressure distribution on the face of the drop is assumed to be hydro- 
static under a head of hp, if the shear stress along the solid boundaries be- 
tween section 1 and 2 is neglected, and if the momentum coefficients 8; and 
Bo are assumed as unity, the momentum equation for sections (1) and (2) is 


? ? 


YY, AZ, YY," (1) 
~ 
+ AZ. QV, ) 


in which 7 is the specific weight of the fluid and q is the discharge per unit 
width of the channel. If the continuity equation, 


q = V,Y, = 


1%1 (2) 


and Eq. 1 are combined and reorganized, the following general equation relat- 
ing the several variables is obtained. 


Plots of Eq. 3 for different values of the variables are shown as Figs. 3, 
4, and 5. 


It can be shown that for the condition hp = Y2 (Fig. 1A) Eq. 3 reduces to 
2 
2 Y 7 4 


A plot of this equation gives the upper limiting curve (plotted as a dashed 
line) for Figs. 3, 4, and 5. 


AZ 


For the condition = Y, +A Z_, that is — 1+ (Fig. 1B) the 
1 


general equation reduces to 


>) 
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FIG.3 RELATIVE DOWNSTREAM DEPTH AS A FUNCTION OF 
FROUDE NO. AND PIEZOMETRIC HEAD AT THE DROP. 
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FROUDE NO. = 
Vg 


FIG.5 RELATIVE DOWNSTREAM DEPTH AS A FUNCTION OF 
FROUDE NO. AND PIEZOMETRIC HEAD AT THE ODOROP 
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> 2 
(%) 
= -[1l+ 
1%) 


Equation 4 and 5 are equivalent to the equations given by Hsu, et al,(5) for 
the two values of hp which were assumed to determine the characteristics 
analytically. 


It is noted that a reduction of YY causes a shift to the right in the curves of 
1 
Figs. 3, 4 and 5. That is, for a given condition of the entering flow, the jump 
can be held in place with a reduced downstream depth Yg. The converse is, 


of course, equally true. 


Experimental Apparatus 


Experimental observations were made on two drops which were con- 
structed of 3/4 inch waterproof plywood and each in turn mounted in a one- 
foot wide glass walled flume. One drop was 0.400 ft. deep and the other was 
0.200 ft. deep. Each was 3 ft. long with an elliptical section at the upstream 
end and a vertical section at the downstream end. A brass sluice gate adjust- 
able by an attached screw device could be positioned to 0.001 ft. The sluice 
gate had a carefully machined 1/2 inch rounded edge which gave only a slight 
contraction in the jet. 

Water was supplied from a constant head tank through an 8 in. pipe to the 
forebay of the flume. The discharge was regulated with an 8 in. gate valve, 
and the flow rate was measured by a 5 1/4 x 8 in. Venturi-Orifice Meter with 
a 4 3/8 in. sharp edged orifice plate inserted in the throat section. The dif- 
ferential head on the meter could be read to ¢ .001 ft. which made it possible 
to set the rate of flow with an error less than 1%. 

Surface elevations were determined with an electric point gage mounted 
upon an instrument carriage that could slide along rails above the channel. 
The vertical position of the gage was read with a direct vernier reading to 

0.001 ft. The depth of water in the downstream channel was controlled by a 
wicket type gate, which caused a minimum disturbance to the flow pattern up- 
; stream. It was located 17.5 feet downstream from the last section at which 

| measurements were taken. 

Velocities were measured with a standard Prandtl pitot-static tube that 
had a coefficient of 1.00. The pitot tube was fitted to the same instrument 
carriage which held the point gage and could be positioned to 0.001 ft. Velo- 
cities were measured near the bottom of the channel at several positions 
downstream from the drop. Distances were measured on a section of engi- 
neers steel tape with an error of no more than ¢ 0.002 ft. 
| Values of hp were measured by means of piezometers attached to 
piezometric holes drilled in a 1/2 in. x 1 in. brass strip set into the face of 
the drop so as to be flush with the vertical surface. The openings were set at 
0.073 ft., 0.200 ft. and 0.327 ft. from the bottom of the 0.400 ft. drop and at 
0.100 ft. from the bottom of the 0.200 ft. drop. 


| 
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The following flow conditions were chosen for observations. 


Z Fy Type of Jump 
0.400 0.100 4 4,6 B,W,A 
0.400 3 2,4 B,W,A 
0. 200 0.067 3 6,8 B,W,A 
0.200 0.100 2 2,4,6 B,W,A 


Values of the depth Yo, the piezometric head at each piezometer, and 
velocity along the bottom were read for each of the flow conditions listed. 


Results 


It was found that the pressure on the face of the drop was hydrostatically 
distributed under a piezometric head, hp, which varied with the location of 


the pump relative to the drop. The measured values of 


=— are shown on 


Figs. 3, 4, and 5 at the values of = and F, at which they were observed. 
1 
The decimal point indicates the plotted point. These observed values confirm 
the analytical calculations within the limits of accuracy of the observation and 
clearly indicate the variation of hp with the change in the flow conditions. 
The limiting values of + within which each type of jump would form at 
1 
given entering flow conditions are shown in Figs. 6, 7, and 8. The conditions 
under which Jump A will form are shown as a very narrow band, but actually 
Jump A may form at values of . greater than the upper bound being limited 
1 
only by the free distance upstream and the slope and resistance of the chan- 
nels. The undular wave will form at values of the variables which fall within 
the middle band. The height of this wave may exceed the value of Yo by as 
much as 50 per cent and might cause considerable damage unless adequate 
training walls were provided. Jump B will form when values of the variables 
are those within the lower band (dotted). The narrow gap between the inter- 
mediate and lower bands in Figs. 6, 7, and 8 represents a transition zone 
where there is uncertainty as to whether the wave or Jump B will form. For 
= 


values of Z and F, which would lie below the lower band, the jump would be 
1 


swept downstream from the drop. Perhaps the most important characteristic 


we for the entire 
range of entering Froude numbers covered by the experiments. This might 


be anticipated from the previous analysis because the hydrostatic head, hp, 
should vary continuously between its limiting values. Visual observation of 
the gradual change in curvature of the main stream of the flow as it crossed 


presented is that the jumps can form over a wide range of 
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FIG.7 FORMS OF THE HYDRAULIC JUMP AS A FUNCTION 
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the drop also indicated a continuous variation of hp. As indicated by Eq. 3, a 


continuous variation of hh implies a continuous variation of 
D 


between its 
1 
extreme values. 


To demonstrate the effectiveness of the different types of jumps in reduc- 


ing the velocity along the channel bottom, the velocity near the channel bottom 
is plotted against distance from the drop in Figs. 9, 10, and 11. Of particular 


V 


Vo 


interest are the high velocities produced by Jump B which may give ratios 


of 5 or 6 at = of approximately 2. The wave and Jump A give values of 
2 


Vv 


v from 2 to 3 in the upstream direction near the drop. The return flow re- 
2 


sults from the reverse roller produced by both the wave and the Jump A. 
Jump A gives the smallest values of velocities along the bottom in the down- 
stream direction. This condition would be expected since the greater pres- 
sure on the face of the drop causes the entering flow to be directed slightly 


upward and the violent eddy motion to occur in the upper levels of the down- 
stream channel. 


Application of Results 


The investigation has clarified the action of an abrupt drop in stabilizing 
the position of a hydraulic jump. The results obtained are directly applicable 
to the design of a stilling basin incorporating an abrupt drop. 


If the downstream depth, Y9, is required for given entering flow conditions, 
it may be read from Fig. 6, 7, or 8 using the chart with the applicable relative 


AZ 
= 


drop height. For example, given Y, = 2 ft., Froude No. = 4 and 


1 


Then, from Fig. 7, the following information could be obtained: 


Y 
Range of —— 


1 


Type of Jump 


Range of Y> 


Jump B 


6.8 
7.0 - 7.6 
7.6 - 7.8 


10.2 - 13.6 ft. 
14.0 - 15.2 ft. 
25.7 


Wave 


Jump A 


The results may be applied to analyze the performance of a stilling basin 
for a range of discharges and tailwater elevations. For example, assume a 
spillway with a crest 100 ft. long at elevation 460 and a coefficient of 4.0. Let 
the apron be at elevation 400 with a drop of 8 ft. to elevation 392. Fora 
series of values for the head on the spillway, discharge may be determined as 
well as the depth and Froude number at the base of the spillway; i.e., Y, and 
Fy. (For convenience, the effects of resistance and air entrainment were neg- 
lected in the example.) By use of Figs. 6, 7, and 8, and by interpolating for 
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AZ 

intermediate values of = , the range of values of Y»g for each type of jump 
1 

can be determined. These results are shown in Fig. 12 in which the depths 
Y9 have been converted into elevations. Also shown in Fig. 12 are curves of 
the tailwater elevations at sections 2 and 3 for a rectangular downstream 
channel 100 ft. wide with its bottom at elevation 398.5, a slope of 0.0016, and 
a Manning n of 0.03. This is referred to as Channel I. In making the compu- 
tations it was assumed that the depressed apron (elevation 392) was 125 ft. 
long which would allow tranquil flow to develop at section 2. If the depressed 
apron were shorter so that high velocity flow would impinge on the step at the 
end of the apron, dynamic pressures would develop on the face of the step 
adding an additional upstream force component tending to stabilize the jump. 
Since it is not yet possible to evaluate the effect of this force, it was elimi- 
nated from consideration by choosing a long apron. 

From Fig. 12, it is apparent that the type A jump will occur at discharges 
above about 20,000 cfs with the form changing to the wave at lower rates of 
discharge. There is no danger of the jump washing downstream for these 
conditions and from Figs. 9, 10, and 11, it is apparent that the velocity near 
the bottom would be negligible while the type A jump forms; i.e., Q > 20,000 
cfs. At lower discharges when the wave forms the velocity near the bottom 
at the end of the depressed apron would not exceed 3V9 or about 3 x 7 = 21 
ft./sec. 

Figure 13 shows what would happen if the same spillway discharged into a 
channel 100 ft. wide with its bottom at the same elevation as the depressed 
apron (392) and with a slope of 0.0010 and a Manning n of 0.035. This is re- 
ferred to as Channel II. From Fig. 13, it is apparent that type A jump would 
form at the maximum discharge with the form changing successively through 
the wave to type B as the discharge decreases. As the form of the jump 
changes, it would move progressively down over the drop but, even at the 
lowest tailwater, it would remain close to the location of the drop. 

These examples illustrate the type of analysis made possible by the re- 
sults of these experiments. By making it possible to determine the jump 
characteristics for a wide range of Froude number and tailwater conditions, 


the effect of both upstream and downstream geometry on the jump character- 
istics can be determined. 


CONCLUSION 


1. An abrupt drop in the bottom of a rectangular channel is effective in 
stabilizing the hydraulic jump over a broad and continuous range of values of 
the relative downstream depth. 


2. As the relative downstream depth is reduced from its upper limit, the 
jump shifts progressively downstream resulting in a corresponding continuous 


decrease in the pressure on the face of the drop as the jump passes through 
three different forms. 


3. A momentum analysis which includes a term for the pressure force on 


the vertical face of the drop gives results which are consistent with the ex- 
perimental observations. 
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4. The maximum velocity attained near the channel bottom, which may be 
considered as an index of the scouring power of the flow leaving the jump, 
varies over a wide range and depends mainly on the form of the jump. For 
the type B jump (Fig. 1) this velocity may attain a value several times greater 
than the downstream velocity, but for the type A jump, the maximum velocity 


near the bottom exceeds only slightly the average velocity in the downstream 
channel. 
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ABSTRACT 


Transition from laminar to turbulent flow was experimentally observed 
during flow establishment in a smooth straight pipe. Turbulence originated 
as spots. The downstream face of the turbulent spots was transported with 
the centerline velocity of the preceding laminar flow. Turbulent spot forma- 


tion was in qualitative agreement with the Tollmein-Schlichting theory of 
small disturbances. 


INTRODUCTION 


A study of flow establishment in a pipe was performed in the Hydraulics 
Laboratory, School of Civil Engineering, Georgia Institute of Technology. 
The experiments were performed with a smooth, straight, circular pipe which 
was horizontally aligned. The upstream extremity of the pipe was a well- 
rounded inlet which was located in a large reservoir. The downstream ex- 
tremity of the pipe was unobstructed so that a liquid jet was formed in the at- 
mosphere downstream from the pipe. The flow was established by rapid 
(effectively instantaneous) release of a disk valve placed against the square- 
edged downstream end of the pipe. Velocity-time data were obtained from a 
motion picture record of the jet. Pressure-time data were obtained at 
selected stations along the pipe. These data were recorded as an oscillo- 
graph record of the output signals of pressure transducers. 

The experiments were performed with a systematic variation of the inde- 
pendent dimensionless variables in the equations: 
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and 


The letter symbols are defined as follows: 


D - inside diameter of pipe; 

g - acceleration of gravity; 

h - piezometric head in the pipe, referenced to the pipe centerline; 

hy - piezometric head in the reservoir, referenced to the pipe center- 
line; 

L - length of pipe; 

Vv - kinematic viscosity of the fluid; 

t - time since opening of the valve; 

V - mean velocity in the pipe; and, 

Xx - axial distance from the inlet. 


Results obtained from ten experimental runs (Runs 20-29, inclusive) are re- 
ported in this paper. The value of L/D was 380.3 in Runs 20-24 and was 
475.3 in Runs 25-29. The other a controlled variable, gh )D3/Lv2, 
was established at values of 6(106), 12(106), 18(106), 24(106), and 30(106). 

In this paper only the phenomenon of the transition from laminar to turbu- 
lent flow is considered. The specific objectives of this study are (1) to deter- 
mine the position of turbulence formation, (2) to determine the manner of ex- 
pansion of the turbulence, and (3) to determine the requisite flow conditions 
for the formation of the turbulence. 


Experimental Equipment and Technique 


Figure 1 is a schematic drawing of the equipment. The description of the 
equipment is presented from the upstream section to the downstream. 

The reservoir was a sealed cylindrical steel tank with a diameter of 42 
inches. The air pressure in the reservoir above the water was controlled to 
any desired magnitude by means of an air overflow line. The air overflow 
was allowed to bubble into the bottom of a standpipe. Thus by controlling the 
water level in the standpipe, the air pressure in the reservoir could be main- 
tained so as to obtain the desired magnitude of the reservoir piezometric 
head, h,. No water was added to the reservoir during a run, and no run was 
made without allowing the reservoir to be quiescent for at least four hours. 
Since the tank cross section was large and the volume of air above the water 
was large, the piezometric-head change was negligible (A ho/hy < 0.06%) 
during the short time required for a run. With this arrangement, the value of 
ghoD3/L v2 was determined within an estimated accuracy of + 0.4% and the 
magnitude was controlled within + 1.0% of the desired value. 

The pipe leading from the reservoir consisted of coupled sections of extra- 
heavy smooth brass pipe. In order to vary the pipe length easily and 
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systematically, five sections approximately 95 diameters in length were ma- 
chined and fitted to insure close abutment and proper alignment when forced 
into the machined-brass sleeve couplings. Each pipe length was carefully 
aligned along a straight horizontal axis. 

The extruded brass pipe proved to be extremely uniform in diameter. 
Micrometer measurements were taken on orthogonal diameters at each end 
section. The mean diameter of these twenty measurements was 0.544 inch 
with all readings being within + 0.001 inch of this value. The average diame- 
ter of each pipe length was also checked by weighing the water that would fill 
each length. 

The rounded pipe inlet was located in the reservoir. The inlet piece was 
attached to the first pipe length by a machined connection similar to the 
sleeve coupling. The machined brass inlet piece was made with a radius of 
rounding equal to 1.13 pipe diameters. A band of non-uniform sand roughness 
approximately 1/8-inch wide was cemented to the inlet piece slightly upstream 
from the tangent point of the inlet roundness. 

Wall piezometer holes were placed 1/2-inch from the downstream end of 
each pipe length. Four piezometer openings were drilled on orthogonal diame- 
ters with a No. 45 (0.0820 inch) drill. The holes were carefully inspected, 
and any burrs caused by drilling were carefully removed. The four piezome- 
ter openings were connected by a piezometer ring which was machined in the 
pipe sleeve coupling. 

Pressure-time measurements were obtained by means of pressure trans- 
ducers, the signal of which was electrically amplified and recorded. These 
differential-pressure transducers were mounted on a transducer collar in or- 
der that the water passage from the pressure transducer to the piezometer 
ring was large in cross section and short in length. The pressure transducers 
were of the unbounded-resistance-wire type. The transducer impulses were 
amplified and recorded on an oscillograph. The oscillograph contained an in- 
dependent timing mechanism with the result that the oscillograph record con- 
sisted of the timing pips as well as the output record of the two pressure 
transducers. Because of the cost of the transducers and the recording equip- 
ment, only two transducers with a two-channel recorder were utilized. Runs 
with four and five pipe sections (Runs 20-29, inclusive) were duplicated in or- 
der to obtain pressure-time records at all piezometer stations. 

In order to simplify the boundary conditions associated with the flow es- 
tablishment, a downstream valve was required which could be instantaneously 
opened and which in the opening process would not influence the resulting 
fluid motion. The physical requirements of such a valve are (a) in a negligi- 
ble time after valve release the piezometric head, h, at the pipe outlet must 
drop from the value, hj, to zero, and (b) the valve must be moved in the same 
direction and with a faster motion than the fluid particles issuing from the 
pipe outlet. The type of valve selected was a flat disk, which was positioned 
flush against the downstream end of the pipe. The force of a stressed tension 
spring was utilized to remove rapidly the valve disk from the oncoming fluid. 

Velocity determinations were made from individual frames of a 35 mm 
motion-picture camera record of the jet. The exposure time for each frame 
was 0.0015 second. The camera was run at a speed of approximately 36 
frames per second. Rear lighting was used in order to eliminate the shadow 
of the jet. Also, a clock image was recorded in each motion picture frame. 
Time could be determined to the nearest 0.001 second. A horizontal and a 
vertical grid scale were also photographed. These scales were positioned in 
the plane of the jet to avoid parallax. 
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Transformation of Data into a Form Suitable for Analysis 


With this experimental equipment, the time-dependent data were available 
in the form of a 35 mm motion-picture film strip of the jet and the oscillo- 
graph records from the pressure-transducer circuit. The reduction of the 
pressure data into a useable form was straightforward. The reduction of the 
film strip data into a useable form was not straightforward. In order to ex- 
plain the process of data transformation, the data of Run 29 is used through- 
out the following description. In Run 29 all five pipe lengths were utilized, 
that is, the value of L/D was 475. The other independent variable, 
gh oD3/L v2, was established at the value, 30(106). 


Pressure Transducer Data 


In each run the pressure data as recorded by the oscillograph was re- 
plotted in the form of the dimensionless piezometric head ratio, h/ho, as a 
function of dimensionless time, t Vgho/L. Time, t, was measured from the 
instant of opening of the valve. The very rapid change in piezometric head, 
h, from the initial value, ho, at the instant the valve was opened could be 
easily discerned on the oscillograph record. Time was measured from this 
point on the record by the pips made by the independent timing stylus. 
Figure 2 is a graph of h/hg at each of the five measuring stations as a func- 
tion of t V gho/L for Run 29. The position of measurement is given by the 
dimensionless distance from the inlet, x/D. Runs 29 and 29a were duplicate 
runs in order to obtain pressure measurements at four stations with the two 
pressure transducers. The pressure at the outlet, x/D = 475, is zero for t>0. 


Film Strip Data 


Six frames of the film strip of Run 29 are shown in Fig. 3. 

Two assumptions were used to determine fluid velocity from the photo- 
graphically recorded jet. First, air resistance was considered negligible. 
Second, the jet was considered to be a succession of fluid masses composed 
of individual fluid particles each having the velocity characteristics of this 
larger group. If the velocity characteristics of the individual particles of a 
fluid mass are identical, the motion can be completely described by the mo- 
tion of the mass center of that fluid mass. In Fig. 4 the dimensionless “ob- 
served velocity,” V'/V 2gho, is plotted as a function of t ¥gho/L. 

The magnitude of the error involved in assuming negligible air resistance 
can be determined by observing the same fluid mass on successive film strip 
frames. The relationships between V' and t determined at + 0.5 ft. and 
+ 2.50 ft. below the pipe were identical, indicating negligible air resistance. 

The correctness of the assumption that the jet can be visualized as a 
series of individual fluid masses is confirmed by the appearance of the jet 
(Fig. 3). In all the photographs, discrete masses of fluid appear to be forming 
in the jet. A logical explanation of why the jet can be considered to be a 
series of fluid masses is as follows. During the time the fluid is in the pipe, 
the external force of boundary shear acts on the fluid at the pipe wall. This 
boundary-shear force results in a system of internal shear forces within the 
fluid and results in a radial variation of linear momentum. The fluid par- 
ticles close to the wall are deficient in momentum. As the fluid issues from 
the pipe into the atmosphere, the external boundary-shear force becomes 
zero. However the internal shear forces are reduced to zero only if the 
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linear momentum is constant in a cross section. The internal shear forces 
persist until the jet has constant linear momentum in the cross section. Thus 
the internal shear forces are the implement for effecting the radial transfer 
of linear momentum. The result of this process is a jet with uniform linear 
momentum in a cross section which is, of course, the assumption utilized in 
determining the “observed velocity,” V'. 

However, the “observed velocity,” V', is not the mean velocity, V, of the i 
pipe flow. Because of the boundary-shear force between the fluid and the pipe 
wall the velocity distribution is not uniform within the pipe. Thus the fluid 
upon leaving the pipe has a greater linear momentum than if the velocity of 
all particles were equal to the mean velocity, V. Since the linear momentum 
is constant in a fluid mass after this fluid mass leaves the pipe, the fluid 
mass attains the velocity, V' , which is greater than V. The two velocities 
are related by the relationship, V = V'/Cm, in which Cp is the momentum 
correction coefficient. Cy) is a function only of velocity distribution and is 


defined by the relationship, C,, = (1/A) J (v/v)? dA, in which v is the point 
A 


velocity and A is the cross sectional area of the pipe. Therefore in order to 
utilize quantitatively the data of Fig. 4 in a study of pipe flow, the velocity . 
distribution had to be determined as a function of time. 

The velocity distribution, or specifically the value of Cy, in unsteady 
laminar flow was determined by a combination of theoretical analysis and 
empirical results. 

A solution of the Navier-Stokes equations was obtained in which the experi- 
mentally determined data relating piezometric head and time were introduced. . 
If the pipe flow is uniform, that is, the velocity is axial, two of the Navier- 

Stokes equations are indicative that (a) the piezometric head is constant in a | 
fluid cross section which is normal to the flow and (b) the axial piezometric- i 


pressure gradient 0 p*/@x is only a function of time. The remaining equation 
is 


Ov a*y av 1 


in which v is the point velocity and P is the mass density. The boundary con- 
ditions are that the velocity, v, is zero at the boundary and that the velocity 
gradient, 0 v/@r, is zero at the pipe centerline. The initial condition is that 
the velocity, v, is zero when time, t, is zero. 

The analogous equation with analogous boundary and initial conditions has 
been solved for temperature distribution in a cylindrical rod.(1) The solution 
for the pipe flow problem is 


t 


v = 2/p (2, g(r)e dy (2) 
n=] ° 
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The various terms are as follows: J, and Jj are Bessel functions of the first 
kind of orders zero and one, respectively; the positive numbers a, are the 
zero roots of Jo; ro is the value of the radius, r, at the pipe wall; £ is the 
dimensionless radial coordinate, r/rp; T is a dummy variable of integration; 
g(T) is - @p*/ox; and Gis 

In order to apply Eq. (2) in the analysis of the experimental results, the 
existence and location of a uniform laminar flow region must be determined 
and, then, the time variation of the piezometric-pressure gradient, 9p*/@x, 
in this region must be determined from the pressure-time data of Fig. 2. 

A physical interpretation of the pressure-time data (Fig. 2) is requisite to 
locating a laminar flow region. Since the fluid starts from rest, the flow is 
laminar in the entire pipe during the initial period after the downstream valve 
is opened. The first external indication of a change in flow regime is the 
abrupt pressure increase at the piezometers. For example, as shown in 
Fig. 2, the first of these abrupt pressure increases occurred at the location 
x/D of 189 at the time t\gho/L of 0.616. The interpretation of these pressure 
increases is that an interface between laminar flow and turbulent flow is pass- 
ing the piezometer at that instant. The reason for the rapid change is that 
the mean kinetic energy of the cross section is greater in the laminar flow 
than in the turbulent flow even though the mean velocity, V, is identical on 
both sides of the interface. The difference in these kinetic energies is due to 
the different velocity distributions. Thus at x/D of 189 the flow was laminar 
prior to tVgho/L of 0.616 and was turbulent during the remainder of the run. 

The time when turbulence passes through the outlet can also be determined 
from the pressure record. In Figs. 2 and 4, the time coincidence of the rapid 
pressure drops and the rapid drop in the “observed velocity” V' is apparent. 
Both occurrences are a manifestation that the laminar-turbulent interface is 
passing the pipe outlet at x/D of 475. Prior to t(gho/L of 1.38 the flow is 
laminar at the outlet. During this time a laminar-turbulent interface is mov- 
ing toward the outlet. The flow is turbulent upstream from the interface. As 
explained in the preceding paragraph, the piezometric head is greater in the 
turbulent flow immediately behind the interface than in the laminar flow im- 
mediately preceding the interface. The outlet of the pipe is a station of fixed 
piezometric head since the liquid jet is projected into the atmosphere. Thus 
as the laminar-turbulent interface passes the outlet, the piezometric-head 
line is rapidly lowered just upstream from the outlet. The effect is a pivoting 
of the piezometric-head line. This rapid pivoting of the piezometric-head 
line results in a piezometric-head drop at every measuring station. The 
magnitude of this drop is directly proportional to the distance from the inlet. 
This effect is clearly illustrated in Fig. 2 at the time, tVgh)/L, of 1.38. 

With this preliminary knowledge concerning the movement of turbulence, 
regions can be definitely located in which the flow is laminar. For example, 
in Fig. 2 between the time, tygh)/L, of 0.715 and 1.06 the flow is laminar be- 
tween the x/D stations 379 and 475, and prior to t (gho/L of 0.715 the flow 
is laminar between the x/D stations 284 and 475. 

The piezometric-head lines are indicative that the laminar flow regions 
are uniform. At constant values of time, t Vgh,/L, the measured values of 
piezometric head were plotted as a function of x/D. Since each plotted 
piezometric-head line was straight at a given time in this laminar flow re- 
gion, the conclusion was that the region was also a uniform flow region as 
visualized in the derivation of Eq. (2). A contributing factor in the existence 
of this uniform downstream laminar flow region is the presence of upstream 
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turbulent flow regions, since an upstream turbulent flow region has the effect 
of isolating the inlet. Thus the downstream laminar flow region is similar to 
flow in an infinitely long pipe since the boundary-layer thickness, 
piezometric-head gradient, and developing velocity profile are functions only 
of time. 

The slopes of the piezometric-head lines of the uniform laminar flow were 
determined and plotted. These slope values are shown as points on Fig. 5 for 
the run being discussed. In order to incorporate this experimental data into 
Eq. (2), an analytic function for g(7) must be selected that approximates the 
empirical results. A very simple function was selected, that is, 0 p*/@x is 
proportional to cos(™t/2t*) in which t*Vgho/L is empirically defined as shown 
on Fig. 5. The value of t*Vgh,/L is dependent upon the experimental parame- 
ters as well as the time and position of turbulence inception. The function for 
g(T) was introduced into Eq. (2) and the integration performed. The resulting 
expression for the time-developing velocity profile is 


(t/1*) 
J, (a,, 8) Leos (77/2) (t/*) + (7 20,2 2*) sin (77/2) (t/t*) -e 
Jy (a,) 1 20,” a*) 
v n= 


1 [ cos (77/2) +(17/2a,,2 a*) sin (77/2) (t/t*) 
[1 20,2 a*)"} 


All terms of Eq. (3) have been previously defined except a * which is v t*/ro2. 

Values of v/V were computed in systematic steps so that by interpolation 
of the computed results, these results could be applied to all experimental 
runs. The first nine zero-roots of the Bessel function were used in the nu- 
merical computations. The range of the other variables were as follows: 
0<B<1, 0.01< 0 *<0.05, and 0.2<t/t*< 1.0. Values of the linear momentum 
correction coefficient, Cp,, were then computed using numerical-graphical 
methods. Hence to determine the temporal variation of Cm for any experi- 
mental run only the value of t*Vgho/L was required. 

The analytical solution does not cover the complete temporal range of 
laminar flow. An empirical method was used to extrapolate the analytical 
solution. The method is based upon the “observed velocity” data. At the time 
of transition from laminar flow to turbulent flow at the pipe outlet, there is a 
sharp discontinuity in the “observed velocity,” as shown in Fig. 4. However 
there can be no discontinuity in the mean velocity. Thus the observed dis- 


continuity must be due to the difference in Cp, of the laminar flow and turbu- 
lent flow. Hence at the time of transition 


(V'/ turbulent 


(c_) 


m’ laminar = 


since (Cm)turbulent iS equal to 1.02. The value of Cy of 1.02 for turbulent 
flow is based upon the seventh-root law of velocity distribution. 
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The mean velocity, V, was then computed from the “observed velocity,” 
V'. The displacement boundary-layer thickness was also computed using the 
following relationship 


v,(F, - v2uwrdr (4) 


in which ve is the centerline velocity and 6* is the displacement boundary- 
layer thickness. The values of V/V 2gho, 5*/rg, and v,/V of Run 29 are pre- 
sented in Fig. 6 as a function of time. 


Analysis of Results 


The specific objectives of this paper pertain to the determination of the 
manner of turbulence inception and the manner in which the turbulent flow 
zones are expanded until the pipe is filled with turbulently flowing fluid. 


Downstream Movement of Turbulence 


The pressure-time records, of which Fig. 2 is an example, of all the runs 
are indicative that the turbulent flow region is moved toward the outlet in a 
systematic manner. As previously explained, the abrupt pressure increase 
at the measuring station is an indication that the laminar-turbulent interface 
is passing the piezometer, and the simultaneous pressure decrease at all 
measuring stations is an indication that the laminar-turbulent interface is 
passing the pipe outlet. A systematic movement of the turbulence is indicated 
in Fig. 2 between stations at which x/D is 284 and 379 and between stations 
379 and 475. The mean transport velocity, V*, (between the measuring sta- 
tions) of the laminar-turbulent interface was computed. The results are pre- 
sented in Table I. These results are indicative that the downstream move- 
ment of turbulence is a physical transport phenomenon in which the turbuleni 
eddies are transported downstream with a velocity equal to the centerline 
velocity of the preceding laminar flow. 


Nature of the Downstream Face of the Turbulent Zones 


A qualitative picture of the nature of the downstream face of a turbulent 
zone can be constructed from the occurrences recorded as the interface is 
moved past the pipe outlet. The appearance of a typical jet shortly after the 
interface has passed the outlet is shown on Fig. 3 (e). The simultaneous 
pressure decrease at all measuring stations as shown in Fig. 2 occurs during 
the time the interface is passing the outlet. The rapid decrease in the “ob- 
served velocity,” V', as shown in Fig. 4 also occurs when the laminar- 
turbulent interface passes the outlet. 

In order to establish the chronological order of events, the times of vari- 
ous occurrences are presented in Table II. Columns 2, 3, and 4 of Table II 
are based upon motion picture data. In columns 2 and 3 are recorded times 
of passage through this outlet of the last laminar flow and of the first turbu- 
lent flow, respectively. The fluid that passed the outlet in the interveing time 
is presumably in the interface region. Based upon this time increment the 
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computed lengths of the interface region are presented in column 4. These 
computed lengths are useful only to establish the order of magnitude of the 
length of the interface region because the time interval between the last lami- 
nar flow and the first turbulent flow is too short to be determined with pre- 
cision. The time of passage of the last laminar flow through the outlet as in- 
dicated by the pressure record is presented in column 5. The time of passage 
of the first turbulent flow through the outlet as indicated by the pressure 
record was not tabulated because the rate of pressure decrease as indicated 
on the record was of the same order of magnitude as the capability of the 
oscillograph to follow rapid changes. The interpretation was that the rapid 
pressure drop was at least as rapid as indicated and was probably more 

rapid than indicated. 

Based upon the previously determined transport velocity and the recorded 
data (Table II) at the time the laminar-turbulent interface passed through the 
pipe outlet, a possible configuration of the interface can be conceived which is 
in basic agreement with the data. Fig. 7 is a schematic drawing of the inter- 
face region. The flow is from left to right with laminar flow preceding the 
turbulent flow. The turbulence in region (a) of Fig. 7 is transported down- 
stream at a velocity, V*, equal to the centerline velocity, vc, of the laminar 
flow. The turbulence diffuses radially outward from the central region (a) in- 
to the boundary layer (b). In this advanced stage of laminar boundary-layer 
development, the boundary layer is unstable for a wide range of disturbance 
frequencies. As the turbulence is diffused into the laminar boundary layer, a 
point is reached at which the entire laminar boundary layer (c) suddenly be- 
comes turbulent. The last indication of laminar flow from the pressure re- 
cord (Table II, Column 5), most closely corresponds with the first indication 
of turbulent flow from the motion picture data (Table II, Column 3). The wall 
piezometric head, h, probably does not rapidly change (d) until the entire 
laminar boundary layer becomes turbulent (c). The deceleration of the 
velocity in the central region as the turbulence penetrates the boundary layer 
(b) is the cause of the radial piezometric head gradient, (e) to (d), necessary 
to effect the radial mass transfer involved in changing from the laminar 
velocity distribution to the turbulent velocity distribution. 


Delineation of the Turbulent Flow Regions 


The spatial and temporal location of the turbulent regions within the pipe 
can be determined from the piezometric-head lines. The details of this 
process are illustrated with the data of Run 29 (Figs. 2 and 6). The values of 
the piezometric-head ratio, h/hy, were determined from Fig. 2 at each of the 
measuring stations at constant values of the time parameter, tV gho/L. Since 
the time that a laminar-turbulent interface passed a measuring station could 
be determined from the pressure record and since this interface was moved 
with the centerline velocity of the laminar flow, the position of the downstream 
extremity of a turbulent region could be calculated at the constant values of 
the time parameter for which the piezometric-head lines were to be plotted. 
The piezometric-head changes at the laminar-turbulent interfaces were cal- 
culated as being the difference in the mean kinetic energy of the laminar and 
turbulent flow. The piezometric-head lines were then plotted as shown in 
Fig. 8 A. 

A rather complex plotting procedure was adopted in Fig. 8 in order to show 
both the temporal and spatial location of the turbulent flow regions. The 
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constant value of t y ‘ghy /L for each piezometric-head line is shown at the 
right hand margin. The neuen of the zero value of h/h, was lowered with 
increasing values of t igho/L in order to avoid a superposition of many 
piezometric-head lines. Since the plotted points at the largest values of x/D 
are for the pipe outlet, the zero value of each piezometric-head line is at 
this level. With the known positions of the downstream extremity of the tur- 
bulent regions and with the slopes of the piezometric-head lines of the lami- 
nar and turbulent flow, the upstream extremity of the turbulent regions could 
be graphically determined. The regions of turbulent flow are clearly deline- 
ated in both time and space on Fig. 8. The turbulent flow regions are below 
the dashed line and the laminar flow regions are above the dashed line. As a 
matter of incidental interest, the last laminar flow of Run 29 (Fig. 8A) at- 
tained a pipe Reynolds number, VD/r, of 5.1 (10)4. 


Inception of Turbulence 


With the turbulent flow regions located both spatially and temporally as 
shown in Fig. 8, the flow conditions at the time of origin of the turbulent spots 
can be qualitatively analyzed with reference to the stability of small disturb- 
ances.(2) Basically, the method is to introduce a two-dimensional periodic 
disturbance into the Navier-Stokes equations and to determine whether the 
disturbance is amplified or diminished as the disturbance is moved down- 
stream. The neutral curve which separates the regions of amplification and 
decay is presented on Fig. 9 for Blasius flow over a flat plate.(3) The ordi- 
nate value is a ratio of the displacement boundary-layer thickness 6* to the 
disturbance wave length A. The abscissa value of Fig. 9 is the boundary- 
layer Reynolds number R 5*. 

The salient features of turbulence inception are apparent in Fig. 8, and the 
pertinent flow parameters are presented in Table Ill. The most prominent 
characteristic is that turbulence started as spots in laminar flow. The 
downstream face of these spots was transported downstream in a systematic 
manner as shown previously. The upstream face of the turbulent spots did 
not indicate a systematic movement along the pipe; but even for the spots in 
which the upstream face moved downstream, this movement was at a lesser 
velocity than for the downstream face. The result is a joining together of the 
spots as shown clearly in Fig. 8A. The pattern of turbulent spot generation 
was approximatel Ay ~~, same for equal values of the independently controlled 
parameter, ghoD*/Ly2. That is, in Runs 29 and 24 the value of ghoD3/Lv2 
was 30(10)6, in Runs . and 23 the value was 24(10)6, in Runs 27 and 22 the 
value was 18(10)6 in Runs - and 21 the value was 12(10)6, and in Runs 25 
and 20 the value was 6(10)®. 

Since the flow throughout the pipe was laminar prior to turbulent spot 
generation, the initial transition from a laminar boundary layer to a turbulent 
boundary layer must have been aided by an instability of the laminar boundary 
layer rather than the external source of free-stream turbulence. The magni- 
tude of R§*, v.6*/y, at turbulence inception (Table II) in comparison with 
the magnitude of R§* required for boundary-layer instability (Fig. 9) is in- 
dicative that turbulence exists soon after the boundary layer is unstable. 
Direct numerical comparison between the values of R§* in Table III and 
Fig. 9 is not possible for two reasons. Firstly, the curve of neutral stability 
of Fig. 9 is for Blasius flow over a flat plate rather than for a temporally 
growing boundary layer ina circular pipe. Secondly, the experimentally 
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determined values of R§* of Table II are the values at the time of turbulence 
inception; whereas, the values of R§* of Fig. 9 are the values at which dis- 
turbances are either amplified or diminished. Since a disturbance undoubted- 
ly must undergo some amplification before the disturbance is transformed in- 
to turbulence, turbulence inception occurs at a later time and at a farther 
downstream position than the time and place of disturbance generation. 

The location of the points of turbulence inception (Fig. 8 and Table III) is 
indicative of the source of the disturbances. 

Turbulent spot number 1 occurred in all runs and in the position range of 
0< x/D< 62. The source of disturbances for turbulent spot number 1 was 
the band of sand grain roughness which was cemented to the bellmouth inlet. 
Since the source of disturbances was at the inlet, the value of x/D at the posi- 
tion of turbulence generation is evidence of the distance required to amplify 
the disturbance before turbulence is actually generated. The flow pattern is 
complex in the region just downstream from the inlet since the boundary layer 
is expanding both with increasing time and increasing distance from the inlet. 
Because the disturbance is generated in a region of small boundary-layer 
thickness, the disturbance may first diminish in amplitude due to the small 
and hence stable value of R§*. Since the boundary layer increases in thick- 
ness both with distance and time, the disturbance which was initially diminish- 
ing may pass into the region of instability or amplification. The disturbance 
is moved further downstream curing amplification. At some point the disturb- 
ance attains a sufficient amplitude to break into turbulence. The heterogene- 
ous mixture of disturbance wavelengths from the nonuniform sand grain 
roughness as well as the spatially varying laminar boundary-layer thickness 
precluded the possibility of avoiding the unstable region (Fig. 9) with the re- 
sult that turbulent spot number 1 existed in all runs. 

Turbulent spot number 2 occurred in all runs except Run 25 and in the 
position range 99 << x/D < 134. This position is indicative that the source of 
disturbances was the pipe junction at x/D of 95. In spite of the care of pre- 
paring the pipe junctions, the pipe edges could have been mismatched as much 
as 0.001 inch because of diameter variations between the ends of each pipe 
length. With the pipe junction as the source of disturbances for turbulent spot 
number 2, the distance required for amplification of the disturbance was from 
4 to 39 pipe diameters. More significantly, this distance for amplification 
varies between 1305* and 6506*. In order to determine the minimum possi- 
ble wavelength the maximum value of 277§*/) of 0.409 is utilized from Fig. 9. 
Thus the minimum wave length,A, is approximately 156*. With this minimum 
wavelength and the distance for amplification, the maximum number of dis- 
turbance waves between the disturbance source and the point of turbulence in- 
ception would be approximately 40. 

The existence of turbulent spot number 2 (Table III) in all runs except Run 
25 can be qualitatively explained by reference to the neutral stability curve. 
The dimensionless boundary-layer thickness, 6 */rg, and the boundary-layer 
Reynolds number, R§*, were comnuted for each run. In Fig. 10 the value of 
§*/To is presented as a function .. R§*, vcd*/v for the early portion of all 
runs. The boundary-layer parameters as shown in Fig. 10 are very similar 
for Runs 29 and 24, Runs 28 and 23, and so forth. This similarity of boundary- 
layer parameters is the reason for the similarity of the turbulent spot pat- 
terns which was mentioned previously. With the hypothesis that the disturb- 
ance wave length A is proportional to the roughness height at the pipe junc- 
tion, a possible neutral stability curve for turbulent spot number 2 is shown 
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on Fig. 10. This neutral stability curve is for a disturbance wavelength, A, of 
approximately 2.1 rg. This neutral stability curve is for the longest disturb- 
ance wavelength which originates at the pipe junction. Thus, the neutral 
stability curves for shorter wavelengths would be below the curve shown in 
Fig. 10. With this concept there is no conflict in that some of the points of 
turbulence inception (indicated by the circular symbols on Fig. 10) are below 
the lower branch of the curve shown. The curve of neutral stability is a line 
of zero amplification. Other lines of constant amplification are similar in 
shape to the neutral curve and resemble the contour lines of the end of a 
ridge on a topographic map. Since the boundary layer is undergoing essen- 
tially temporal growth, the boundary-layer curves of Fig. 10 can be visua- 
lized as pertaining to a single wave or disturbance. Thus a disturbance of 
Runs 20 and 25 (Fig. 10) would be subjected to lesser amplification in cross- 
ing the unstable zone than for the other runs. The disturbance of Run 20 re- 
ceived sufficient amplification to break into turbulence, whereas the disturb- 
ance of Run 25 passed into the upper stable zone before breaking into 
turbulence. Once a disturbance is in the upper stable zone, this disturbance 
will diminish without breaking into turbulence. 

Turbulent spot number 3 occurred in only Runs 24 and 29 at the x/D posi- 
tion of 200 and 213. This position is indicative that the source of the disturb- 
ances was the pipe junction at x/D of 190. This axial distance is equivalent to 
3106 * and 6006 *. The maximum number of disturbance waves between the 
source and the point of turbulence would be approximately 40. A possible 
neutral stability curve for the longest wavelength, equal to 1.5 rg, is shown in 
Fig. 10 for turbulent spot number 3. The interpretation is the same as before 
except that only in Runs 24 and 29 was a disturbance amplified sufficiently to 
break into turbulence. In all other runs a disturbance was not sufficiently 
amplified in crossing the unstable zone to break into turbulence. 


Upstream Face of a Turbulent Spot 


The experimental data are of little use in defining the nature of the up- 
stream face of the turbulent zones. Since every individual turbulent eddy is 
transported downstream, new turbulence would have to be generated even for 
the upstream face of the turbulent zone to remain in a fixed position. Since 
the position of the upstream face of the turbulent zones depends upon the 
generation of new turbulence, the probability is that this face has an erratic 
motion. If two turbulent spots are joined before the downstream face of the 
upstream spot passed a pressure transducer station, the presence of the up- 
stream spot could not be detected from the data. In fact, it is probable that 
turbulent spot number 1 of Fig. 8B and Fig. 8D is formed from two or more 
spots which are soon joined. The upstream face of the turbulent zone would 
move downstream if no new turbulence were generated. In a temporally 
growing boundary layer, a single burst of turbulence might be generated, sub- 
sequent disturbances from a fixed source of disturbance being insufficiently 
amplified to break into turbulence. Referring to Fig. 10, this occurrence ap- 
pears to be probable for turbulent spot number 2 in Run 20 and for turbulent 
spot number 3 in Runs 24 and 29. The definite downstream movement of the 


upstream face of these turbulent spots, as shown in Figs. 8A, 8B, and 8J, is 
indicative of this occurrence. 
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CONCLUSIONS 


Based upon an experimental study of unsteady flow, the following conclu- 
sions concerning transition from laminar to turbulent flow in a pipe are con- 
sidered to be significant: 


1. In an experimental investigation of laminar boundary-layer transition 
the effects of free-stream turbulence can be avoided entirely with a temporal- 
ly expanding boundary layer. 


2. Transition from laminar to turbulent flow in a pipe occurs in spots. 


3. The magnitudes of the laminar boundary layer Reynolds number at 
turbulence inception was in reasonable agreement with the theoretical values 
determined for Blasius flow over a flat plate. 


4. The downstream movement of a turbulent zone into a laminar zone is 
that of physical transport. The transport velocity is that of the centerline 
velocity of the laminar flow. In the region of this interface the free stream 
turbulence is the disturbance source for laminar boundary layer transition. 


5. The regions of laminar and turbulent flow in a pipe can be defined from 
the piezometric-head lines. 


6. From the standpoint of boundary shear, the pipe utilized in the experi- 
ment would be classified as smooth. However, from the standpoint of laminar 
boundary-layer transition, disturbance sources were found at the pipe junction 
as well as at the sand-roughened inlet. 


7. With a systematic variation of boundary-layer parameters, the presence 
or absence of boundary-layer transition at the pipe junctions was found to be 
in qualitative agreement with the laminar boundary-layer stability theory. 
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(a) T = 0.355 SEC. (d) T = 1.222 SEC. 


(e) T = 1.607 SEC. 
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(c) T = 0.511 SEC. (f) T = 1.953 SEC. 


Figure 3. Selected Film Strip Data of Run 
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Piezometric Head Lines of Runs 20-29, Inclusive. (Continued) 
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Figure 8. Piezometric Head Lines of Runs 20-29, Inclusive. (Continued) 
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Piezometric Head Lines of Runs 20-29, Inclusive. (Continued) 
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Table I 


December, 1957 


Downstream Transport Velocity of Turbulence” 


Run No. Value of x/D at 
Measuring Station 


Ratio of Transport 
Velocity to Centerline 


From To Velocity v/v, 
2la 189 284 0.97 
22 284, 380 0.99 
22a 189 380 0.96 
23 28h, 380 1.00 
23a 189 380 0.98 
23b 28h, 380 1.01 
2h 28h 380 0.97 
27 189 379 0.97 
27b 189 379 1.01 
189 379 1.00 
270 189 379 1.02 
27f 189 379 1.03 
28a 189 379 1.00 
29 28), 475 1.01 
29a 379 475 1.02 


“Data from Run Numbers 20, 25 and 26 as well as a portion of the 
data from Run Numbers 21, 27, and 28 have been omitted from this tabula- 
tion. The omitted data occurred during the time the jet was separating 
into individual fluid masses having different values of linear momentum. 
In other words, the radial transfer of linear momentum by the internal 
shear forces was incomplete and some fluid masses fell from the jet 
with a lesser momentum than the main body of the jet, The calculated 
values of the centerline velocity Vv, were questionable during this type 


of jet separation, 
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Table II 


1450-29 


Laminar - Turbulent Interface at the Time of 


Passage Through the Pipe Outlet 


Data from Pictures Data from Pressure Record 

Last Laminar First Turbulent Last Laminar 

20 1.532 1.566 9D 1.587 

21 1.447 1.493 11D 1.501 

22 1.386 1.479 23D 1.468 

23 1.433 1.462 7D 1.460 

-- -- -- 1.214 

25 1.781 1.835 17D 1.844 

26 1.531 1.553 7D 1.583 

27 1.526 1.567 13D 1.579 

28 1.530 1.567 11D 1.58) 

29 1.333 1.365 9D 1.386 


*the length of the interface region was very short with the result 
that meaningful determinations of the time of last laminar flow passage 
and of the time of first turbulent flow passage were impossible. 
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Turbulent 
Spot 
Number 


Table III 


Position of 
Turbulence 
Inception 
in Diameters 


Boundary 
Layer 
Thickness 
o*/r, 


0.100 
0.120 


0.078 
0.062 


0.070 
0.062 


0.068 
0.062 


0.060 
0.052 
0.064 


0.087 


0.078 
0.075 


0.067 
0.061 


0.060 
0.058 


0.060 
0.054 
0.076 


Flow Parameters at the Inception of Turbulence 


Boundary 


Layer 


Reynolds 


R 
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1450-30 HY 6 S| 
Number 
20 1 112 1140 
2 || 
21 1 1030 
2 99 
22 1 62 1100 ; 
2 105 710 
23 48 1020 
2k 1 37 1090 
2 100 640 
3 200 1340 
25 1 4 | 770 
26 1 4 1080 
2 110 960 
27 1 3 920 
250 
28 1 13 810 ; 
29 1 16 1030 
2 105 660 
3 213 2130 
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FLOOD FREQUENCIES DERIVED FROM RAINFALL DATA 
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(Proc. Paper 1451) 


ABSTRACT 


The feasibility of using rainfall data to obtain flood frequencies has been 
demonstrated before. Most investigations, however, utilized rainfall data and 
rainfall-runoff relations specific to the basins for which the flood-frequency 
data were desired. The present study, involving 36 basins, shows that the 
rainfall data need not be from the problem basin and that the data from one 
precipitation network can be used to synthesize streamflow data for a number 
of basins. Similarly, rainfall-runoff relations already derived for other 
basins may be easily modified to apply to the problem basin. The transpo- 
sition of rainfall data and rainfall-runoff relations not only make it possible 
to extend the streamflow record of basins not having rainfall data, but greatly 
reduces the cost of synthesizing streamflow data. 


INTRODUCTION 


Among the procedures for extending strea ped records for design purpos- 
es is that reported by Kresge and Nordenson, 1) who demonstrated that 
streamflow data, specifically annual flood peaks, could be synthesized from 
rainfall data by means of a rainfall-runoff relation and unit hydrograph. Us- 
ing the Hocking River at Athens, Ohio, area 944 sq. mi., which had a 4-station 
precipitation network and a 37-year period of concurrent precipitation and 
streamflow records, they showed that the flood-frequency characteristics of 
the basin could be reproduced from rainfall data. The procedure requires a 
great deal of work and the cost of application on a broad scale might be con- 
sidered prohibitive, especially in connection with the design of low-cost hy- 
draulic structures, such as highway culverts. The results were so encourag- 
ing, however, that in 1954-56 the U. S. Weather Bureau and the Bureau of 


Note: Discussion open until May 1, 1958. Paper 1451 is part of the copyrighted 

Journal of the Hydraulics Division of the American Society of Civil Engineers, 

Vol. 83, No. HY 6, December, 1957. 
1. Chief, Cooperative Studies Section, U. S. Weather Bureau, Washington, D.C. 
2. U.S. Weather Bureau, Washington, D. C. 
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Public Roads jointly undertook further investigations along these lines to de- 
termine if the procedure could be modified to lower the cost of application 
when used on a broad scale. 

A preliminary study based on three basins: Delaware River at Valley 
Falls, Kans. (922 sq. mi.); Meramec River near Steelville, Mo. (781 sq. mi.); 
and Ralston Creek at Iowa City, Iowa (3 sq. mi.), showed that their flood- 
frequency characteristics could be satisfactorily reproduced from rainfall 
data outside the basins. The tests were made by taking rainfall data from a 
number of networks around each basin and comparing the results of the 
frequency analysis of the 30 years of annual flood peaks at each basin with 
those of the annual flood peaks synthesized from the rainfall data for the 
corresponding networks for the same period as the streamflow record. It 
was found that a precipitation network would yield satisfactory results provid- 
ed (1) its area roughly approximated the size of the problem basin, (2) its 
rainfall-intensity —frequency characteristics were similar to those of the 


basin, and (3) the seasonal distribution of rainfall was the same as that of the 
basin. 


These preliminary studies thus led to the large-scale undertaking of synthe- 


sizing streamflow data for 36 basins (Table 1). These 36 basins, ranging in 
size from 2.2 to 521 sq. mi., are scattered in a long, narrow region having 
homogeneous physiographic characteristics* and extending from northern 


New Jersey to northern Georgia (Fig. 1). Each basin has at least 15 years of 
streamflow record. 


Procedure 


Selection of Precipitation Networks 


As it had been decided to synthesize streamflow data for no less than a 50- 
year period for each problem basin, the basic precipitation networks were 
made up of stations with at least 50 years of concurrent rainfall data. This 
restriction severely limited the number of networks available. Because of 
the relatively wide scatter of stations with 50-year rainfall records, the indi- 
vidual networks were limited to no more than 3 stations, since a larger 
number would have made the network area incompatible in size with the 
problem basins. 

The effect of the size of the network area on synthesized peak discharges 
was investigated. For this purpose the West Fork River Basin at Clarksburg, 
W. Va., 384 sq. mi., was used. This basin contains a fairly evenly distributed 
precipitation network of 8 stations, and is surrounded by other stations not 
far outside its boundaries. Fig. 2 shows the equation used to obtain an index 
of the effective area of the 1-, 2- and 3-station networks selected for this 
investigation and the effect of network area on synthesized 2-year peak dis- 
charges. The 2-year peak discharge synthesized from the storm-precipitation 
averages of the basin’s 8-station network served as the control. The curve 
of Fig. 2 indicates that for network areas of about 40 to 600 sq. mi., a differ- 
ence up to 350 sq. mi. between network and basin areas affects the results by 
no more than 5 percent. The range of network area from 0 to 40 sq. mi., 


*As indicated by the map of Problem Areas in Soil Conservation, published 
by the U. S. Soil Conservation Service in July 1950. Most of the problem 
basins were in the area designated as the Piedmont Plateau, or B-10, region. 
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Table 1 


PROBLEM BASINS 


STREAM AND GAGING STATION 


Whippany River at Morristown, N. J. 
North Branch Raritan River nr Far Hills, N. J. 
Walnut Brook nr Flemington, N. J. 
Neshanic River nr Reaville, N. J. 
Tohickon Creek nr Pipersville, Pa. 
Perkiomen Creek at Graterford, Pa. 
Chester Creek nr Chester, Pa. 

Codorus Creek at Spring Grove, Pa. 
Octoraro Creek nr Rising Sun, Md. 

Piney Run nr Sykesville, Md. 

Linganore Creek nr Frederick, Md. 

Seneca Creek at Dawsonville, Md. 
Northwest Branch Anacostia River nr Colesville, Md. 
Goose Creek nr Leesourg, Va. 

Slate River nr Arvonia, Va. 

Tye River nr Lovingston, Va. 

Appomattox River at Farmville, Va. 

Goose Creek nr Huddleston, Ve. 

Sandy River nr Danville, Va. 

Mayo River nr Price, N. C. 

Fishing Creek nr Enfield, N.C. 

Dial Creek nr Bahama, Ne C. 

Eno River at Hillsboro, N.C. 

Horsepen Creek at Battleground, N.C. 
East Fork Deep River nr High Point, N. C. 
Deep River nr Randleman, N.C. 

Fisher River nr Copeland, N.C. 

South Yadkin River nr Mocksville, N.C. 
Uwharrie River nr Eldorado, N. C. 

Brown Creek nr Polkton, N.C. 

South Tyger River nr Woodruff, S.C. 
South Tyger River nr Reidville, S.C. 
Middle Tyger River at Lyman, S.C. 
Apalachee River nr Buckhead, Ga. 

Little Tallapoosa River at Carrollton, Ga. 
Tobesofkee Creek nr Macon, Ga. 


ASCE 1451-: 
No. 
(From ee AREA 
Mic. 1) (Sg. Mi.) 
29 
26.2 
2024 
25.7 
94.4 
279 6 
61.1 
74.3 
193. 
1 11.4 
62.3 
12 101. 
13 23.23 
14 336. 
15 
16 92. 
17 306. 
18 187. 
19 113. 
20 260. 
el 
| 22 4.9 
1 23 66.5 
24 159 
| 25 14.2 
26 124, 
27 121. 
28 313. 
29 347. 
30 110. 
3A 174. > 
32 106. 
33 68.3 
34 436. 
35 89. 
36 182. 
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however, encompasses a rapid change in selectivity, suggesting that synthe- 
sized peak discharges should be adjusted when there is any appreciable differ- 
ence between basin and network area within this range. The data of Fig. 2 
also suggest that 2-station networks yield results as satisfactory as those 
obtained from 3-station networks. Use of the 2-station instead of the 3-station 
network results in a great saving of time in making storm surveys and de- 
termining time distribution of rainfall. Best results for very small basins, 
i.e., under about 15 sq. mi., were obtained from single stations. This is true, 
of course, only because closely-spaced long-record precipitation stations 
were not available. Another important consideration in the selection of pre- 
cipitation networks was their geographic proximity to groups of problem 
basins, the objective being to effect economy by applying one network to a 
group of basins. 


Storm Surveys 


After the precipitation networks had been selected, their records of daily 
precipitation were scrutinized to select the major rainstorms of each water 
year. A base value of daily precipitation was arbitrarily selected (usually 
one inch) and whenever one or more stations in the network exceeded the base 
value, precipitation was tabulated for all stations for the entire storm period, 
including days with less than the base value. If fewer than 10 or 12 storms 
were obtained in any one water year, the base was lowered and the record re- 
examined. 

The runoff for each major rainstorm was then roughly estimated by means 
of a rainfall-runoff relation (Fig. 3) as described by Kohler and Linsley. (2) 
The four or five storms producing the highest runoff in any one water year 
were then set aside for a detailed analysis of the time distribution of the rain- 
fall, because the relatively small size of the problem basins required a break- 
down of rainfall into increments no longer than 6 hours for accurate estimates 
of runoff and computations of peak discharge. 


Time Distribution of Rainfall 


The time distribution of rainfall in the several storms selected for each 
year was determined from the records of a recording-gage station in or near 
the network involved. The hourly data were scanned to determine the maxi- 
mum 6-hour rainfall increment for a storm, and the entire storm was divided 
into 6-hour increments synchronizing with the maximum. The average storm 
rainfall for the network was distributed into 6-hour increments in accordance 
with the indicated distribution for the recording-gage station. Six-hour incre- 
ments of runoff were then obtained from a rainfall-runoff relation. For basins 
smaller than about 15 sa. mi., the peak discharges generally result from 
short-duration, high-intensity rainfall. Consequently, 1-hour increments of 
rainfall were used for such basins, the 1-hour runoff cs being esti- 
mated by means of a special rainfall-runoff relation. (3 


Rainfall-Runoff Relations 


Rainfall-runoff relations were available for very few of the problem basins. 
Development of these relations is a time-consuming procedure. In order to 
keep the cost of the study to a minimum, available rainfall-runoff relations 
from other basins were utilized. In order to determine which of several 


ASCE PAULHUS - MILLER 1451-5 


relations was best suited to a particular problem basin, the runoff for 25 or 
more storms over the basin was computed from the hydrographs. The run- 
off for the same storms was also estimated by means of the rainfall-runoff 
relation being tested, and plotted against the corresponding observed runoff. 
The relation yielding values nearest to a 45° line was the most representative 
of the rainfall-runoff characteristics of the problem basin. The mean curve 
for the points yielded by this most representative relation was used to adjust 
all runoff values obtained from that particular rainfall-runoff relation for use 
on the problem basin involved. Fig. 4 presents such a mean curve indicating 
so high a degree of representativity for the rainfall-runoff relation that no 
adjustment was considered necessary. 


Unit Hydrographs 


Unit hydrographs of 1-hour or 6-hour duration were derived from discharge 
data obtained from the original records of the U. S. Geological Survey for 
converting runoff increments into peak discharge. 4) It is generally recog- 
nized that one unit hydrograph for a basin may net properly reproduce the 
stream hydrograph for all magnitudes of floods. 9) Different unit hydrographs 
are often required for low and high runoff, and sometimes for intermediate 
values. No matter what criteria are established for the choice of the proper 
unit hydrograph, situations will exist where almost equal values of runoff will 
require different unit hydrographs, probably resulting in appreciable differ- 
ences in the peak discharges obtained. In order to avoid such discontinuities, 
only one unit hydrograph was constructed for each basin, and a curve was 
then developed to adjust the peak discharges yielded by the unit hydrograph. 

The unit-hydrograph adjustment curve for a particular basin was derived 
by first obtaining from observed flood hydrographs the total storm runoff for 
a good range of peak discharges. The next step was the breakdown of storm 
runoff into 1- or 6-hour incremental runoff on the basis of the chronological 
distribution of the rainfall producing the runoff. The peak discharge was 
then computed by applying a unit hydrograph like that of Fig. 5 to the runoff 
increments. The computed peak discharge for each storm was then plotted 
against the observed peak discharge for the storm. The mean curve (Fig. 6) 
was used to adjust the synthesized peak discharges estimated by the unit 
hydrograph, 

In addition to eliminating the need for several unit hydrographs for one 
basin, the peak-discharge adjustment curve also allowed for base flow. Base 
flow was not included in synthesizing the peak flow from the unit hydrograph, 
but it was included in the observed values of peak discharge used to derive 
the adjustment curve. Consequently, a part of the adjustment provided by 
application of the curve was for base flow. 


Results 


Comparative Test 


After selection of the most representative rainfall-runoff relation and unit 
hydrograph and their adjustment curves for a particular basin, peak discharg- 
es corresponding to the rainstorms for the period of discharge measurements 
were estimated. The estimated or synthesized annual peak discharges, i.e., 
the maximum for each water year, were arrayed in decreasing order of 
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magnitude and then compared with the observed, similarly arrayed. If the 
average difference between the two was less than 20 percent, the record of 
synthesized peak discharges was considered to be representative of the 
streamflow record. Annual peak discharges were then synthesized for the re- 
mainder of the precipitation record. 


Synthesized Annual Peak Discharges for Codorus Creek, Pa. 


Table 2 shows such a comparison between synthesized and observed annual 
peak discharges (annual floods) for Codorus Creek Basin above Spring Grove, 
Pa., area 74.3 sq. mi. Codorus Creek has a streamflow record starting with 
water year 1933, and annual flood peaks were available through water year 
1954. Peak discharges were therefore synthesized for the corresponding 
period, 

The peak discharges for Codorus Creek were synthesized from major rain- 
storms observed by a 3-station network (Ephrata, Lancaster and Lebanon, 
Pa.) about 45 miles northeast of the basin. The time distribution of the rain- 
fall in the selected storms subsequent to 1940 was determined from either or 
both of the recording gages (Lancaster, Lebanon) and prior to 1940 from the 
recording gage at Harrisburg, about 30 miles west of the network. The run- 
off from these storms was estimated by means of a rainfall-runoff relation 
previously developed for the Anacostia River Basin, about 50 miles south of 
the problem basin. This rainfall-runoff relation was so well adapted to the 
problem basin that no adjustment was considered necessary (see Fig. 4). 

The 6-hour unit hydrograph (Fig. 5) used to estimate the peak discharge 
from the computed runoff was based on five storms of approximately 6 hours’ 
duration in the problem basin. An adjustment curve was required (Fig. 6). 

The first three columns of Table 2 show poor agreement on an annual basis 
between observed and synthesized annual floods. This is not surprising as 
it would indeed be strange if floods of the same approximate frequency oc- 
curred simultaneously over the problem basin and over the network 45 miles 
away. Tests based on two pairs of adjoining basins showed poor agreement 
between observed annual peak discharges for adjoining basins on a year-to- 
year basis. 

The last two columns of Table 2, showing observed and synthesized annual 
floods in decreasing order of magnitude, suggest that the two may well repre- 
sent samples of the same population. The synthesized annual floods are 
generally well within the 20-percent limit of acceptability. A frequency analy- 
sis of the observed and synthesized annual floods for the 22-year period by 
the Gumbel method yielded the results shown in Table 3. 

The frequency data of Table 3 might be criticized by some statisticians 
because they were based on all annual floods for the 1933-54 water years. 
Referring to Table 2, it is seen that the 1933 observed and synthesized annual 
floods were both extremely high, being about twice the second highest. Some 
would argue that the 1933 annual floods are not of the same population as the 
lesser floods and should therefore not be included in the frequency analysis. 
Table 4 shows the results of the frequency analysis of observed and synthe- 
sized annual floods for water years 1934-54. This table shows even better 
agreement between observed and synthesized annual floods than when the 
1933 data were included. In any event, it appears that the record of synthe- 
sized annual floods compares reasonably well with that of observed annual 
floods for Codorus Creek Basin. 


— 
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Table 2 


ANNUAL FLOODS FOR CODORUS CREEK AT SPRING GROVE, PA. 
(74.3 Sq. Mi.) 


CHRONOLOGICAL LISTING : STATISTICAL ARRAY 
Water Observed Synthesized : Order Observed Synthesized 
Year (cfs) (cfs) (cfs) (cfs) 
1954 1910 2890 : 1 11200 11600 
1953 2180 22k0 : 2 6070 5370 
1952 1750 2340 : 3 4880 4900 
1951 3200 4900 y 4180 4800 
1950 1610 2590 5 3200 3910 
1949 1730 1640 : 6 3180 3670 
1948 1810 1930 : 7 3070 2890 
1947 1430 1320 : 8 2710 2810 
1946 3070 4800 : 9 2490 2730 
1945 2490 2810 ¢ 2380 2590 
1944 2380 3910 ; & 2210 2360 
1943 4180 1750 ‘.. 22 2180 2340 
1942 2710 5370 33 1910 22k0 
1941 870 1340 — 1810 1970 
1940 4880 1970 : 2S 1750 1940 
1939 1620 1940 ; 8 1730 1930 
1938 3180 2730 7) RT 1620 1780 
1937 2210 2360 2 1610 1750 
1936 1380 1780 1430 1640 
1935 1110 1330 : = 1380 1340 
1934 6070 3670 = 1110 1330 


1933 11200 11600 oS 870 1320 
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The suitability of the above precipitation network, rainfall-runoff relation 
and unit hydrograph for synthesizing peak discharges for Codorus Creek 
Basin having been established, the next step was to synthesize annual floods 


for the remainder of the precipitation record. The results are given in Table 
5. 


The results of a frequency analysis of the 55-year period of synthesized 
annual floods of Table 5 are given in Table 6. Also listed are the results ob- 
tained with the 1933 flood eliminated. Both sets of data suggest that the 1933 
observed flood peak of 11,200 cfs is an unusual event with a return period 
probably much greater than 100 years. 


Synthesized Annual Floods for the Other 35 Project Basins 


Synthesized annual floods for the other 35 project basins (Table 1) provid- 
ed comparable degrees of agreement with the observed for the respective 
basins. Altogether, only 13 precipitation networks and 17 rainfall-runoff 
relations were used in synthesizing annual flood records for all 36 basins. 


CONC LUSIONS 


This study substantiates the conclusion that reasonably reliable flood- 
frequency data may be synthesized from rainfall records. Moreover, the 
precipitation stations need not be in the basin for which the synthesized flood 
peaks are required. The fact that the rainfall record from one network and 
one rainfall-runoff relation can apparently be used to synthesize flood- 
frequency data for a number of basins makes the procedure economically 
practicable. 

One point that must not be overlooked, however, is that there were no ade- 
quate precipitation networks in the project basins. This meant that there 
were no entirely suitable controls for determining whether the synthesizing 
procedure was yielding reliable values of peak discharge storm by storm. It 
is obvious that tests involving even but a few years of streamflow data with 
concurrent rainfall observations from an adequate network in the basin proper 
are to be preferred to the tests employed in these investigations. 

In spite of this shortcoming, the results of the studies made thus far 
strongly suggest that the synthesizing of peak discharges from rainfall is 
practicable and economically feasible. 
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Table 3 


FOR 1933-1954 


Observed 
(cfs) 


2530 
4850 
6380 
8320 
9750 

11180 


Table 


FOR 1934-54 


Observed 
(cfs) 


2280 
3610 
4500 
5610 
6440 
7260 


FREQUENCY ANALYSIS OF OBSERVED AND SYNTHESIZED ANNUAL FLOODS 


Synthesized 
(cfs) 


2730 
5040 
6570 
8500 
9940 
11360 


FREQUENCY ANALYSIS OF OBSERVED AND SYNTHESIZED ANNUAL FLOODS 


Synthesized 
(cfs) 


2470 
3730 
4560 
5610 
6390 
7160 
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Table 5 


SYNTHESIZED ANNUAL FLOODS FOR CODORUS CREEK AT SPRING GROVE, PA. 
FOR 1900-1954 


CHRONOLOGICAL LISTING STATISTICAL ARRAY 


Water 
Year 


Annual 
Flood 


Order 
No, 


Annual 
Flood 


1954 2890 : 1 11600 

1953 2240 : 2 5420 1908 
1952 2340 : 3 5370 1942 
1951 4900 4 4900 1951 
1950 2590 5 4800 1946 
1949 1640 6 4270 1902 
1948 1930 7 4100 1913 
1947 1320 8 3910 1944 
1946 4800 : 9 3670 1934 
1945 2810 : 10 3420 1916 
1944 3910 : 11 3040 1915 
1943 1750 : 12 2950 1927 
1942 5370 3 13 2890 1954 
1941 1340 : 14 2850 1903 
1940 1970 : 15 2810 1945 
1939 1940 : 16 2730 1938 
1938 2730 : 17 2660 1926 
1937 2360 : 18 2590 1950 
1936 1780 : 19 2360 1937 
1935 1330 : 20 2340 1952 
1934 3670 : 22k0 1953 
1933 11600 : 22 2210 1906 
1932 2200 : 23 2200 1932 
1931 1390 : ok 2200 1924 
1930 2030 : 25 2140 1919 
1929 1380 : 26 2120 1907 
1928 1530 : 27 2030 1930 
1927 2950 : 28 1970 1940 
1926 2660 : 29 1960 1900 
1925 1380 : 30 1950 1917 
1924 2200 : 31 1940 1939 
1923 1200 : 32 1930 1948 
1922 1400 : 33 1820 1910 
1921 1220 : 34 1810 1905 


1920 


(Continued next page) 
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7 12 : 35 1780 1936 
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Table 5 - Con. 


SYNTHESIZED ANNUAL FLOODS FOR CODORUS CREEK AT SPRING GROVE, PA. 
FOR 1900-1954 


CHRONOLOGICAL LISTING : STATISTICAL ARRAY 

Water Annual Order Annual Water 
Year Flood $ No. Flood Year 

1919 2140 : 36 1750 1943 
1918 1520 : 37 1730 1911 
1917 1950 : 38 1640 1949 
1916 3420 : 39 1600 1904 
1915 3040 : 4o 1590 1914 
1914 1590 : 41 1530 1928 
1913 4100 : ke 1520 1918 
1912 1250 : 43 1470 1901 
1911 1730 : 4h 1400 1922 
1910 1820 $ 45 1390 1931 
1909 1270 : 46 1380 1929 
1908 5420 : 47 1380 1925 

a 1907 2120 : 48 1340 1941 
1906 2210 : kg 1330 1935 
1905 1810 : 50 1320 1947 
1904 1600 : 51 1270 1909 
1903 2850 : 52 1260 1920 
1902 4270 : 53 1250 1912 
1901 1470 : 54 1220 1921 
1900 1960 : 55 1200 1923 

Table 6 


FREQUENCY ANALYSES OF SYNTHESIZED ANNUAL FLOODS 


Return period 1900-54 1900-32 


(years) 1934-54 
2 22h0 2160 
5 3840 3220 
10 4890 3920 
25 6230 4810 
50 7220 5470 


100 8200 6130 
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PRECIPITATION NETWORKS 


Plainfield-Elizabeth, N. J. 
Ephrata-Lebanon-Lancaster, Pa. Ns 
Emmitsburg, Md.-Gettysburg, Pa, 
Baltimore-Woodstock, Md. 
Harper's Ferry, W. Va. -Lincoln- 
Mt. Weather, Va. 

Columbia-New Canton, Va. 3 
Henderson-Louisburg, N. C. 


Rae 


“Abaitimo Brevard-Hendersonville, N. C, 
Gainesville-Maysville, Ga. 


Talbotton-Woodbury, Ga. 
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(Based on streamflow measurements 
near Riverdale and Hyattsville, Md) 
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RAINFALL-RUNOFF RELATION FOR ANACOSTIA RIVER 
FIGURE 3 
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COMPUTED VS OBSERVED STORM RUNOFF 


CODORUS CREEK AT SPRING GROVE, Pa. (74.3 Sq. Miles) 
FIGURE 4 
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SIX-HOUR UNIT HYDROGRAPH 
CODORUS CREEK AT SPRING GROVE, PA., (74.3 Sq. Miles) 
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COMPUTED VS OBSERVED PEAK DISCHARGE 
CODORUS CREEK AT SPRING GROVE, PA., (74.3 Sq. Miles) 
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COMPUTED PEAK DISCHARGE (1,000 cfs) 
COMPUTED VS OBSERVED PEAK DISCHARGE 
CODORUS CREEK AT SPRING GROVE, PA. (74.3 Sq. Miles) 
FIGURE 6 
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FLOW CHARACTERISTICS ON THE OGEE SPILLWAY 
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(Proc. Paper 1452) 


SYNOPSIS 


Although the ogee spillway has been the subject of more research than per- 
haps any other hydraulic structure, there are still deficiencies in our 
knowledge of its functioning. There is more to be learned about the process 
of energy loss on the spillway face and about the effects of submergence of 
the spillway jet. A general formula is here suggested for calculation of the 
spillway energy loss, and a theoretical derivation is presented defining the 
limit between the true hydraulic jump and the drowned jump. A criterion for 
estimating bottom velocities along a level spillway apron is considered. 


INTRODUCTION 


The determination of energy losses in flow on the face of an ogee spillway 
and the solution of the energy dissipation problem on the spillway apron often 
have been approached by means of hydraulic models. The reliability of such 
tests is universally accepted. The hydraulic laboratory plays a most im- 
portant role in the design of these structures. The wealth of data on spillway 
energy loss produced by the many laboratories offers a good source of ma- 
terial for use in the development of general mathematical criteria. This 
paper attempts to bring this potentiality into focus. 

Experience gained thus far in design of the overflow spillway indicates 
that this hydraulic structure has certain elements which lend themselves to 
standardization. In shaping the ogee crest, for example, it has become wide- 
ly accepted practice to conform to the lower nappe of a stream issuing over 
a sharp-edged weir. Uniformity has been obtainable in this case because of 


Note: Discussion open until May 1, 1958. Paper 1452 is part of the copyrighted 
Journal of the Hydraulics Division of the American Society of Civil Engineers, 
Vol. 83, No. HY 6, December, 1957. 


*Supervising Hydr. Engr., Dept. of Water Resources, State of California, 
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the easily predictable nature of the phenomenon of overflow. Probably due to 
the relative difficulties of measurement, the prediction of flow characteristics 
on the spillway face and apron has not been simplified to this extent. How- 
ever, by correlation of the results of experimentation by our hydraulic labora- 
tories, such simplification should be possible. 


Energy Loss on Spillway Face 


A simple and exact criterion for determining the energy loss which will 
occur on the face of an overflow spillway would be a valuable design tool. 
Various mathematical procedures are now employed for this purpose. The 
application of the Manning formula has many advocates and is believed by 
them to yield satisfactory results. Solution by this means, though, is not di- 
rect and there are differences of opinion as to how the formula can best be 
applied. For example, in using the Manning equation, a common practice has 
been to determine the friction loss that would occur for flow at the depth with- 
out air entrainment, but to allow for an assumed lower actual loss by reducing 
the value of the Manning coefficient of roughness. Yet, in the words of one 
reviewer of this paper, it is not even universally agreed that entrainment of 
air in the water has the effect of reducing the friction loss. 

Flow characteristics on the apron below the dam will depend, among other 
factors, upon the energy that remains in the jet when it enters the tail water. 
The height of the spillway, the head on the crest, the shape of the crest, and 
the position of the tail water will all affect the amount of energy thus involved. 
The theoretical energy content of the jet will obviously vary directly with the 
vertical distance of travel. The amount of water involved, and therefore the 
total energy rate, will depend primarily upon the crest head and the crest 
form. 

It is reasonable to assume that the energy loss on the spillway face is de- 
termined by a simple natural law. The factors involved are few, and their 
influences are relatively uncomplicated. The principal components of the 
formula would be the rate of discharge and the height of fall. The roughness 
and slope of the spillway face also affect the energy loss, and would have to 
be taken into account. The influence of these factors could probably be in- 
corporated into a coefficient in the formula. 

In the Madden Dam spillway prototype tests, (1) Pitot-tube measurements 
of jet velocities were made for a range of crest heads. The loss in energy 
head on the spillway face was taken as the difference in elevation between 
the reservoir water surface and the energy grade line calculated on the basis 
of the Pitot-tube readings. In order to appraise the accuracy of these ve- 
locity tests, the method of construction of the Pitot tubes should be noted. 

As described by Mr. Randolph, each installation consisted of a plate steel 
pedestal embedded in the concrete of the spillway apron, and projecting from 
14 to 20 inches above the floor. The Pitot tube itself was a length of 1-inch 
pipe protruding 4 inches from the nose of the pedestal at a height between 3 
inches and 18 inches from the floor. Although readings on these Pitot tubes 
were thoroughly and carefully taken, it can not be assured that the resulting 
accuracy is comparable with that obtained with laboratory instruments. Also, 
since the tubes were located close to the floor, they gave indications of the 
maximum velocity rather than the mean. Despite these apparent limitations, 
the writer considers the Madden Dam velocity tests to be reliable enough to 
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serve in a preliminary general analysis of basic relationships. 

The head loss between the spillway crest and any point on the spillway 
face, for a given concrete roughness, should be a function of the height of 
fall and the rate of discharge over the crest. It is suggested, therefore, that 


the free-flow loss on the spillway face may be defined by the general ex- 
pression: 


” (P_+ H}® 
(1) 


wherein K is a coefficient, P is the height of crest above the downstream 
apron, H is the total head on the crest, q is the discharge per unit length of 
crest, and “a” and “b” are constant exponents. Fig. 1 shows the relationship 
of the basic dimensions. 

Application of this general formula to the prototype observations on the 
spillway of the Madden Dam shows encouraging possibilities. Using values 
of K = 1/8, a = 3/2, and b = 1/3, an acceptably close correlation can be ob- 
tained. However, the face of the Madden spillway has a roughness somewhat 
greater than average. For a typical concrete spillway of common roughness 
and comparable in other respects to Madden Dam, it is believed that the 
following equation will give a reasonably close measure of the head loss. 


_ 0.10(P + 
h, = 5 (2) 


Substitution of CH3/2 for q in Equation (2) leads to 


0.10(P + 
hy, (3) 


in which C is the coefficient of discharge over the crest at the given head. 
Mr. Randolph in his paper attempted to derive an expression for the 
friction loss on the ogee spillway face. The following two formulas, convert- 


ed to the notation used herein, were developed by him for the Madden Dam 
spillway: 


h, = (P + H) 4) 


_ 0,00934 P(P + H) 


he = 
L 0.6 (5) 
23 


Equation (4) is dimensionally correct but Equation (5) is not. The deficien- 
cies of these expressions, either individually or comparatively, are demon- 
strated by the substitution of the Madden Dam spillway height of 172 feet and 
a crest head of one foot. The result in each case is an energy loss greater 
than the actual head available. It is believed that Equation (3) will be found 
to be a more generally applicable conclusion. The implication is not intend- 
ed, however, that Equation (3) itself is a final refinement. It is merely an 
approximation independently drawn from the basic data presented by Mr. 
Randolph in his highly commendable work. 

Mr. Randolph, in his paper on the Madden Dam tests, stated that he was 
suggesting the installation of Pitot tubes on the aprons of other dams designed 


| 
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by the U. S. Bureau of Reclamation. If data from such Pitot-tube measure- 
ments are available for a full range of heights of prototype dams, a formula 
probably of a form similar to Equation (1) can likely be derived which will 
remove the uncertainties from this design problem. The use of data obtained 


from high prototype spillways will enable the extension of the criteria to in- 
clude the effect of air insufflation. 


Type of Flow on Spillway Apron 


In order to provide full assurance against structural failure at the toe of 
an overflow dam, the desirable flow characteristics on the spillway apron 
must be known, along with the means of maintaining these characteristics. 
Formation of a true hydraulic jump may be desired, but not always feasible 
to attain. When the tailwater depth exceeds a certain limit, the jump can be 
expected to become submerged and the efficiency of energy dissipation ac- 
cordingly reduced. Of importance to the designer, therefore, is the determin- 
ation of the limiting condition at which the jump will lose its true character- 
istics. 

The criterion of demarcation between the true hydraulic jump and the 
drowned jump can be established theoretically. This is demonstrated in the 
analysis which follows. 

Using the term H to denote the total energy head on the crest (including 
velocity head) and D, as the depth of flow at the base of the overfall, the ve- 
locity at the base of the overfall is 


wherein hy is the total head loss from the crest to the apron, and g is the ac- 
celeration of gravity. 


The depth of flow at the base of the fall can be expressed as 
ae /2e(P + H- D, - hy) (7) 
The hydraulic jump will become drowned at a tailwater depth greater than 
the sequent depth Dg which would occur at the downstream end of a true jump. . 


The limit between the true jump and the submerged jump can therefore be 
defined by this sequent depth. The theory of conservation of momentum gives 


2 
Dp = Ph -1 


(8) 


By substitution of terms derived above, the following dimensionless ratio 
is obtained: 


q 
| 
v, = /2e(P + H- Dy - (6) 
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Dy} + H- - 


Do 


P x (9) 


By plotting the ratios H/P versus D9/P, obtained through application of 
Equations (7) and (9), a family of limiting curves for various values of the 
coefficient of discharge can be developed. Fig. 2 serves as an illustration. 
For convenience in determining the depth at the end of the overfall, Fig. 3 
may be used instead of Equation (7). 

In planning for the design of the Boulder Canyon Project, the U. S. Bureau 
of Reclamation conducted an extensive hydraulic model testing program, much 
of which was devoted to a study of performance of overflow spillways. This 
program has provided experimental data relating to the limit between the 
true jump and the drowned jump. 2 Fig. 4 was derived from these obser- 
vations by the Bureau of Reclamation. Comparison of Figs. 2 and 4 indicates 
that the limiting depth observed in the laboratory was somewhat higher than 
that obtained by the theoretical formulas. This difference is noted even when 
the loss on the spillway face is neglected in applying Equations (7) and (9). It 
is suggested that perhaps Fig. 2 defines the limit at which the jump first tends 
to move upstream, and that the limiting curve of Fig. 4 marks the tailwater 
depth at which the drowning process is complete. 


Bottom Velocities on Spillway Apron 


Of primary concern to the designer is the distance traveled by the high- 
velocity stream under the tail water. On this will depend the extent of the 
works necessary for protection of the streambed against erosion. It is im- 
portant that a criterion be established to locate the point at which the velocity 
of the jet has been reduced to a safe magnitude. - 

In the testing program for the spillway of Madden Dam, a set of measure- 
ments was taken of bottom velocities along a horizontal apron below a model 
overflow dam. An analysis has been made of the data reported from these 
tests. The relations drawn in Fig. 5 show the most promise in interpreting 
the experimental results. The graph is dependent upon limited data and there- 
fore should not be applied to general use. It is presented here to illustrate 
the possibilities of such analysis. The curves in Fig. 5 apply to a level apron 
without baffle piers or deflector sill. The variables used in the graph are 
defined as follows: The term “x” is the horizontal distance from the inter- 
section of the downstream spillway slope with the level apron (neglecting the 
curve of the bucket); Vo is the velocity of flow at the point where the stream 
on the spillway face passes the elevation of the tail water at the downstream 
end of the apron; d, is the depth of flow, measured perpendicular to the spill- 
way face, at the same point as Vj; hyo is the velocity head corresponding to 
Vy; D is the depth of tail water at the downstream end of the apron; and V is 
the bottom velocity of the jet on the apron. 


CONC LUSION 


This presentation has been exploratory in nature. The primary factors 
affecting the characteristics of spillway flow have been considered, along with 
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suggested means of relating them. Basic relations developed here are de- 
pendent upon a limited experimental foundation and are therefore intended to 
serve only for comparison with other model or prototype observations. It is 
expected that the discussions which are invited will produce further records 
from observed spillway flows that may enable the formulation of more wide- 
ly usable criteria. 

Particular emphasis has been placed herein on the determination of ve- 
locities of flow on the spillway face and apron, and on the phenomena influ- 
enced by these velocities. It is believed that a simple natural law governs 
flow on the face of an ogee spillway, although the exact relationship has not 
yet been written. The formula for energy loss here offered for comparison 
may help in the writing of a dependable rule that can be applied generally to 
overflow spillways of any shape or size. 

It is hoped that this paper will provoke discussions which will contribute 
knowledge not yet available to the profession at large. If the skeleton of cri- 
teria appearing here is given sustenance drawn from the experience of others, 
the purpose of this work will be served. 
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SYNOPSIS 


A comprehensive solution is proposed for the discharge characteristics of 
rectangular, thin-plate weirs. The solution is based on a simple equation of 
discharge and experimentally derived coefficients which account for the influ- 
ence of the fluid properties and the physical characteristics of the weir and 
the weir channel. The effects of viscosity and surface tension are related to 
an increase in the effective head and an increase or decrease in the effective 
notch width. Thus, the combined effects of the fluid properties are accounted 
for with adjustment coefficients which are applied to measured values of the 
head and width. Consequently, the coefficient of discharge is defined as a 
function of the width-contraction ratio and the head-weir-height ratio only. 

Values of the coefficients required to describe the discharge of water over 
rectangular, thin-plate weirs are determined from original tests covering a 
wide range of the significant geometric ratios. Additional data needed to con- 
firm the proposed method of analysis are obtained from the published results 
of experiments by other investigators. The proposed discharge equation is 
compared with some of the most widely used formulas. 


INTRODUCTION 


The rectangular, thin-plate weir is defined in this report as a basic weir 
form which includes rectangular-notch weirs as well as full-width or 
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“suppressed” weirs. A typical rectangular-notch weir in a rectangular 
channel is shown in Fig. 1. 


Fig. .— THe Basic THin-PLate NotcH WEIR 


For many years hydraulicians have sought a comprehensive equation that 
would describe the discharge characteristics of rectangular weirs for a full 
range of fluid, flow, and geometric variables. Because weir discharge in- 
volves free, curvilinear flow and the combined influences of several fluid 
properties, it is not subject to complete mathematical description. Most in- 
vestigators have resolved their difficulties by restricting their research to 
the flow of water over full-width weirs of considerable height and width. Un- 
fortunately, however, different workers continue to disagree regarding the 
physical limitations as well as the discharge characteristics of even this so- 
called standard measuring weir. Thus, there still exists a variety of weir 
formulas, each based on a different set of experimental data, which disagree 
so much that the weir is questioned as an accurate measuring instrument. 

Most published weir formulas are related to a quasi-rational analysis in 
which the weir is described as a limiting example of the two-dimensional 
orifice. The traditional form of the basic equation of discharge is obtained 
by integrating an approximate velocity equation over an approximate area. 
The result is a complex equation which must be modified to make it practic- 


able and adjusted to make it agree with a particular set of experimental data. 
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The fact that a large number of different empirical formulas have been evolved 
by this procedure is evidence that the “theoretical” equation is inadequate. 

The obvious reason for the inadequacy is the fact that the basic assumptions 
are inaccurate. 

Because the flow pattern for the weir is not yei subject to mathematical de- 
scription, the most direct solution for the discharge function involves a com- 
bination of dimensional analysis and experiment. Such an approach to the 
problem is used in this report. 


Analysis of the Discharge Function 


Description of the Weir 


Fig. 1 shows a typical example of the basic weir form which is the subject 
of this investigation—a symmetrical, sharp-edge, horizontal-crest, rectangular 
notch in a smooth, vertical, thin plate located in a smooth, horizontal, rec- 
tangular channel. The nappe is fully ventilated and unsubmerged. The channel 
upstream from the weir, described hereafter as the approach channel, is as- 
sumed to be of sufficient length to develop the normal velocity distribution for 
all discharges. This ideal condition is seldom realized in practice. In fact, 
the use of short channels with various arrangements of baffles to control the 
velocity distribution is a major cause of disagreements in the results of ex- 
periments made by different investigators. 


Dimensional Analysis 


The geometry of the flow pattern shown in Fig. 1 is described by the width 
of the notch (b), the width of the approach channel (B), the height of the weir 
crest above the bottom of the channel (P), and the piezometric head (h) re- 
ferred to the level of the crest and measured in the uniform flow section up- 
stream from the weir. The fluid properties involved are the specific weight 
the density (P), the viscosity (u), and the surface tension Only one 
independent flow characteristic is involved, and it can be represented by either 
the discharge (Q) or the head (h). However, h is already involved as an inde- 
pendent geometric variable, and Q is conveniently selected as the dependent 
variable. Thus, a complete statement of the discharge function will include 
both h and Q, as in the expression 


Q £,(b, B, P, h, (1) 


From Eq. 1, dimensionless ratios which describe the discharge function 
can be formed as follows: 


(2) 
in which g = 7 //, the acceleration due to gravity. The dependent ratio in Eq. 
2 is a coefficient of discharge. The first three independent ratios describe 
the geometry of the flow boundaries, and the last two symbols represent the 
Reynolds and Weber numbers. 
In American engineering practice the acceleration due to gravity is includ- 
ed in the definition of the coefficient of discharge. Thus, for practical pur- 
poses, a convenient definition of the coefficient (C) is 
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(3) 


from which 


Q = c (4) 


The dimensions of C in Eqs. 3 and 4 are the dimensions of Vg. Because of 
its obvious simplicity, and despite its lack of dimensional purity, Eq. 4 is used 
as the basic discharge equation in this report. 

The discharge function represented by Eq. 3 cannot be evaluated by analyti- 
cal procedures. Accordingly, the relative influence of each of the independent 
ratios must be evaluated by experiment. 


Significance of the Geometric Ratios in Eq. 3. 


One of the most significant ratios in Eq. 3 is b/B. It is a measure of the 
channel-width-contraction characteristic of the weir. In combination with h/P, 
it is also an area-contraction ratio. The influence of b/B is similar to that 
of the corresponding width or diameter ratios which are used to describe the 
geometry of orifices. A better analogy is the width ratio which is used to de- 
scribe open-channel constrictions. 1) Strangely, in view of the great volume 
of recorded research on weirs, the relative influence of b/B has received little 
attention. The weir which has been investigated most extensively in the labo- 
ratory is the full-width weir, for which b/B = 1.0. 

The b/h ratio can be described as a shape parameter. The independent 
influence of this ratio is believed to be negligible over the full practical range 
of the other variables. An earlier investigation 2) at the Georgia Institute of 
Technology supports this conclusion. The fact that most of the published re- 
sults of research on full-width weirs ignore the b/h ratio indicates either that 
it was not recognized as an independent ratio or that its influence was not evi- 
dent from the experimental data. A few recorded efforts to incorporate the 
b/h ratio in discharge formulas are believed to be based on misinterpretations 
of influences related separately to the magnitudes of b and h. 

The h/P ratio is a measure of the depth-contraction characteristic of the 
weir. It is complementary to b/B as an area-contraction ratio. Its influence 
is said to be represented by a velocity-of-approach term in several published 
formulas. It has also been described as a Froude number for the flow in the 
approach channel. Actually, for the boundary conditions specified, the flow in 
the approach channel is always tranquil, and flow with minimum specific ener- 
gy always occurs in the vicinity of the weir crest. In other words, the Froude 
number is not a significant independent ratio. (The Froude number would be 
significant if, for example, the approach channel were sloped or if the weir 
were located at the foot of a spillway or sluice.) Thus, to describe the h/P 


ratio as anything but a geometric parameter is to disguise its fundamental 
significance in the discharge function. 


Significance of R and W in Eq. 3 

Very little is known about the separate influences of viscosity and surface 
tension which are represented by the Reynolds and Weber numbers in Eq. 3. 
It is known, however, that the relative influence of the two fluid properties in- 
creases as the head on the weir and the size of the weir decrease. 


The Reynolds number (R) is a measure of the relative influence of viscosi- 
ty. It is usually expressed as 
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in which V is a typical velocity, L is a conveniently evaluated and physically 
significant length, and y is the kinematic viscosity of the fluid. For large, 
full-width weirs the most significant length is the head. For small, narrow 
weirs, however, the width (and possibly the height) of the weir, as well as the 
head, are independently significant. In general, therefore, because the Rey- 
nolds number must be related to both the head and the width of the weir, there 
are at least two forms of R. Furthermore, because velocities in the vicinity 
of the weir crest are proportional to the square root of the head, Vh can be 
substituted for V in Eq. 5, whence Vh(h)/v and Vh(b)/v are independent ratios 
which have the significance of Reynolds numbers. Thus, for a given fluid (V 
constant), the total influence of viscosity on a given weir can be expressed in 
terms of the absolute magnitudes of h and b, or h alone. 

The Weber number (W) is a measure of the relative influence of surface 
tension. It is usually expressed as 


Vv JL 
(6) 


in which the surface tension (a) and the density () are essentially constant 
for a given liquid. Because V is proportional to Vh, and because either h or b, 
or both, can be significant lengths (L), the ratios Vh(Vh)/Vo/o and Vh(Vb)/ 
Vo/p are independent ratios which have the significance of Weber numbers. 
Thus, there are at least two independent forms of the Weber number as well 

as the Reynolds number. Furthermore, the effects of both viscosity and 
surface tension for a given weir and liquid are related to the absolute magni- 
tudes of h and b. Therefore, it is impossible to distinguish the separate ef- 
fects of the two fluid properties from experiments with a single liquid. 

Because the Reynolds number for any flow pattern is inversely proportional 
to a typical viscous-shear force, the relative influence of viscosity decreases 
as R increases. Similarly, because the Weber number is inversely proportion- 
al to a typical surface-tension force, the relative influence of surface tension 
decreases as W increases. It follows from the foregoing definitions of R and 
W in terms of h and b that the relative influence of the combined fluid proper- 
ties diminishes as either h or b becomes larger. Thus, for large heads on 
large weirs, the influence of viscosity and surface tension is negligible. This 
conclusion is substantiated, for example, by the observation that the coefficient 
of discharge for large weirs is essentially constant for all heads above a 
certain minimum. In other words, C is independent of the fluid-property para- 
meters for large values of h and b. Here the distinction between “small” and 
“large” must be based on systematic experimental investigations of the inde- 
pendent relationships between C and h and between C and b. Many investigators 
have attempted to determine the relationship between the absolute head and 
the discharge coefficient for full-width weirs. The writers know of no system- 
atic effort to evaluate the independent effects related to the absolute width of 
the weir. 

If attention is henceforth restricted to water at ordinary temperatures, it 
follows from the preceding discussion that the quantities h and b can be used 
effectively to replace both R and W in Eq. 3. Furthermore, if the b/h ratio 
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is discarded on the basis of experimental evidence that it is insignificant, the 
discharge function represented by Eq. 3 can be changed to 


c= Bh, b) (7) 


which must be evaluated by experiment. 


Evaluation of Eq. 7 

The function represented by Eq. 7 involves four independent variables, a 
condition which precludes the presentation of experimental data in simple 
graphical form. However, if the relatively minor effects represented by h 
and b are described by an approximate, independent, empirical equation, the 
coefficient of discharge (C) can be expressed as a function of b/B and h/P 
alone. In other words, if the parameters representing the influence of viscosi- 
ty and surface tension are effectively “removed” from Eq. 7, the relationship 
between the coefficient of discharge and the principal geometric ratios can be 
shown by a family of plane curves. 

The effects of viscosity and surface tension cannot be described by exact 
physical equations. Nevertheless, any procedure whereby these effects are 
eliminated from the discharge function must be based on a general understand- 
ing of the manner in which the fluid properties influence the flow pattern. 


The Influence of Viscosity 

A flow pattern which is determined by boundary conditions alone is de- 
scribed as a potential-motion pattern. The influence of viscosity in a real 
fluid motion can be illustrated by comparison with its potential-motion counter- 
part. 

One of the effects of viscosity which can be ignored in many accelerated 
fluid motions is the energy loss which results from internal shear. The total 
loss of energy between the section of head measurement and the crest of a 
weir, for example, is generally insignificant. Therefore, the principal effects 
of viscosity on weir iiow are those which are associated with flow-pattern 
modifications due to boundary resistance and separation. 

A separation zone occupied by a large eddy occurs in the corner between 
the weir plate and the bottom of the approach channel, as shown in Fig. 2(a). 
The effect of this occurrence on a low weir, in comparison with the potential 
pattern for the same weir, is similar to the effect of sloping the weir plate 
downstream. That is, compared with its potential-motion counterpart, the 
vertical trajectory of the nappe is lower and the coefficient of discharge, being 
virtually a coefficient of contraction, is larger. The comparable effect on the 
discharge, therefore, is the same as an increase in head. 

A separation zone similar to that shown in Fig. 2(a) also occurs in the 
corner between the sides of a notch-weir plate and the walls of the approach 
channel. The effect of this occurrence, in comparison with potential motion, 
is to reduce the width contraction of the weir nappe. Its influence on the dis- 
charge, therefore, is the same as an increase in the width of the notch. 

Viscous shear causes the flow to be retarded in the vicinity of the bounda- 
ries. This condition can be described with respect to potential motion by 
means of the discharge-displacement, boundary-layer thickness, 6 * , shown 
schematically in Fig. 2. The cross-hatched areas in the figure represent the 
virtual displacement of the potential motion due to boundary resistance. On 
the surface of the weir plate, near the vertical edges of the notch as well as 


=| 


ASCE KINDSVATER - CARTER 1453-7 


Potential 
Real velocity 


(b) 


| 


— 


Fic. 2—INFLUENCE OF ViscosiTy AND SurFACE TENSION (SCHEMATIC) 


the horizontal crest, the effect of the boundary layer is similar to the effect of 
rounding the upstream corners of the notch. Compared with potential flow 
through a sharp-edged notch, therefore, the coefficient of contraction for the 
real fluid motion is larger. The effect on the discharge, on the other hand, is 
similar to the effect of increasing both h and b. 

A boundary layer is also formed on the bottom and side walls of the ap- 
proach channel. The effect of this occurrence is particularly significant for 
large values of h/P and b/B. For example, the flow through full-width weirs 
(b/B = 1.0) is characterized by displacement boundary layers on both sides of 
the channel at the weir section. The effect of this occurrence on the discharge 
is similar to the effect of a reduction in the width of the weir (b). For values 
of b/B less than 1.0 the effect of the side-wall boundary layers is to increase 
the effective value of b/B by decreasing the effective channel width (B). Simi- 
larly, the influence of the boundary layer on the bottom of the approach channel 
resembles the influence of a decrease in P. The effect of this occurrence on 


the discharge function, therefore, is similar to the effect of an increase in 
h/P. 


A 
A 
| | 
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The influence of viscosity has been associated with the occurrence of 
boundary layers and separation zones. It should be emphasized, however, that 
the thickness of the boundary layers and the size of the separation zones are 
not directly proportional to the magnitude of the head or the size of the weir. 
(For this reason, the effects of viscosity are often described as “scale” ef- 
fects.) Thus, for example, considering the small absolute magnitude of 6 *, 
the relative effect of each of the several boundary layers depends on the abso- 
lute magnitudes of h, b, and P. This observation is consistent with the previ- 
ous conclusion that the relative effects of viscosity diminish as the size of the 
weir and the head on the weir increase. 


The Influence of Surface Tension 

Lindquist, 3) discussing the influence of surface tension on thin-plate, full- 
width weirs, described two significant phenomena. He observed, first, that 
the nappe clings to the top edge of the crest as shown in Fig. 2. Thus, com- 
pared with its potential-motion counterpart, a real fluid flow is characterized 
by a lower trajectory and a larger coefficient of contraction because of this 
surface-tension phenomenon. The effect of the occurrence is similar to the 
effect of rounding the upstream edge of the plate. Therefore, the effect on the 
discharge is the same as an increase in head. 

The second phenomenon described by Lindquist is the result of surface 
tension in both the upper and lower nappe surfaces. The effect of this oc- 
currence is a resultant force, F,, shown schematically in Fig. 2(b), which 
acts in the direction of the center of curvature of the nappe. Because the radii 
of curvature of the nappe surfaces decrease with decreasing head, and because 
the surface-tension force varies inversely with the radius of curvature, the 
resultant force F, increases as the head decreases. Furthermore, because 
the resultant surface-tension force has a dominant downward component, it 
has the same effect on the discharge as an increase in head. 

The surface-tension phenomena described for full-width weirs are equally 
effective on notch weirs. In addition, similar phenomena occur at the vertical 
sides of the notch as well as the crest. Thus, the total influence of surface 
tension on the discharge through rectangular-notch, thin-plate weirs can be 
compared with an increase in the width of the notch as well as an increase in 
the head on the crest. It should be emphasized, however, that the relative 
influence of surface tension depends on the absolute magnitudes of h and b. 
Therefore, like the influence of viscosity, the influence of surface tension di- 
minishes as the size of the weir and the head increase 


A Comprehensive Discharge Equation 

By comparing the flow of a real fluid with the potential flow for the same 
boundary conditions, the influence of the fluid properties has been related to 
an increase in the effective head, an increase or decrease in the effective 
width of the weir, and a decrease in the effective height of the weir. It has 
been suggested also that the coefficient of discharge can be expressed in terms 
of the principal geometric variables alone if the effects of the fluid properties 
are described by approximate, empirical equations. It is proposed, therefore, 
that the effective values of h and b be defined as 


= h+ ky 


(8) 


bo = b+ ky (9) 
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in which the quantities kj, and kp represent the combined effects of the several 
phenomena attributed to viscosity and surface tension. It is also proposed 
that a coefficient of discharge which is independent of the size of the weir or 
the magnitude of the head be defined as 


be hed/2 (10) 


It is assumed that kj, and kp as well as C, can be defined by experiment. Thus, 
from Eq. 10, the equation 


3/2 
6, (11) 
is proposed as a comprehensive equation of discharge for all rectangular, thin- 
plate weirs. The adequacy of this equation, and, in fact, the adequacy of the 
effective-head and effective-width concepts, must be examined in the light of 
experimental data. 


Sources of Verification Data 


Experiments of Bazin 

The classic experiments of Henri Bazin'4; 5) on full-width, thin-plate weirs 
were made in 1886 and 1887 in a flume located on the bank of the Canal de 
Bourgogne, near Dijon, France. Bazin’s setup included a sluice-controlled 
forebay and a concrete flume 6.56 ft wide and approximately 700 ft long. His 
weirs consisted of a crest plate mounted on a timber bulkhead. The nappes 
were fully ventilated on the downstream side. The crest pieces were described 
as being iron, carefully straightened, 0.276 in. thick, sharp-edged, and with 
top surface perpendicular to the plane of the weir. Several different weirs 
were actually used in making the tests, but they were described as being 
“perfectly comparable with each other.” All of the weirs were of the full- 
width type (b/B = 1.0). In order to cover a full range of values of h/P, Bazin 
made ten series of tests. The scope of his work is shown in Table 1. 

Series 1, 2, and 3 of the Bazin experiments were described as reference 
tests, because they alone involved volumetric measurements of discharge. 

For all other tests the weirs were calibrated by comparison with the reference 
weirs. Head measurements were made with hook gages in stilling wells con- 
nected to large orifices in the walls of the flume 16.4 ft upstream from the 
weirs. 

In general, Bazin’s instrumentation was good and his technique was meticu- 
lous. His results differ considerably from those of most other investigators, 
however, and this is usually attributed to the conclusion that, for different 
reasons, the crest pieces were not effectively thin and sharp-edged. (6) 


Experiments of Schoder and Turner 

One of the most extensive investigations of full-width weirs reported in the 
literature is that performed under the direction of. E. W. Schoder and K. B. 
Turner at Cornell University from 1911 to 1920, (7) Nearly 2500 tests were 
made in this investigation, but many of these involved velocity distributions 
which were not uniform and crest plates which were not sharp-edged. The 
two series of tests which have been selected for use in this report are believed 


— 
— 


1453-10 HY 6 December, 1957 


TABLE 1. 


- SCOPE OF THE BAZIN SXPERINENTS 


Series P in feet b in feet h in feet 
l 3.72 6.56 0.19 to 1.01 
2 3.72 3.26 0.19 to 1.34 
3 3.30 1.64 0.19 to 1.78 
4 2.47 6.54 0.48 to 1.28 
5 2.47 6.55 0.49 to 1.44 
6 1.64 6.54 0.31 to 1.22 
7 1.64 6.55 0.34 to 1.40 
& 1.16 6.53 0.29 to 0.94 
9 1.14 6.55 0.30 to 1.34 
10 0.79 6.55 0.30 to 1.34 
to have been made with “normal” velocity distributions. The scope of the - 


tests used (Series E and F) is shown in Table 2. 


TABLE 2. - SCOPE OF THE SCHODER-TURNER SERIES E AND F TESTS (b/B = 1.0). 


Series P in feet b in feet h in feet 
E 0.50 4.22 0.029 to 1.41 
F 1.00 4.22 0.029 to 2.21 


The Cornell tests described in Table 2 were made in a wooden flume 4.22 
ft.wide and approximately 30 ft.long. Because of the short length of the flume, cs 
baffles were provided at the inlet to control the velocity distribution. A mov- 
able bottom was used to obtain different weir heights. Discharges were 
measured volumetrically. The crest piece for this group of tests was made of 
rolled brass, 1/4 in. thick. It was described as having a straight, square, and 
sharp upstream top corner, with a machined 3/32-in. wide top surface and a 
beveled downstream face. The crest piece was mounted on a timber bulkhead. 
Head measurements were made with a float gage connected to piezometers 


11.74 — from the weir. Complete details are given in the paper 
cited. 
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Experiments of the U. S. Bureau of Reclamation 

Experiments on full-width weirs were performed by the U. S. Bureau of 
Reclamation in connection with hydraulic investigations for the Boulder Canyon 
Project. 8) A distinctive feature of the equipment used for the tests was a 
knife-edge, stainless-steel crest piece. The purpose of this innovation was 
to prevent the clinging-nappe phenomenon described previously in this report. 
The crest piece was mounted on a steel plate located in a sheet-metal-lined 
wooden flume 2 ft.wide and effectively 18 ft-long. A movable floor was used 
to obtain different values of weir height. Baffles located at the upstream end 
were used to produce a uniform velocity distribution in the flume. Discharges 
were measured with venturi meters in the supply line. Heads were measured 
with a hook gage in a stilling well connected to a piezometer in the movable 
floor 9.9 ft.upstream from the weir. The scope of the USBR tests is shown in 
Table 3. 


TART 
al 


SLE 3. - SCOPE OF THE USBR EXPERIMENTS (b/B = 1.0). 


Series P in feet b in feet h in feet 
1 5.00 2.007 0.14 to 1.26 
2 4.60 2.007 0.46 to 1.25 
3 3.50 2.007 0.14 to 1.22 
4 2.50 2.007 0.14 to 1.25 
5 2.00 2.007 0.24 to 1.19 
6 1.50 2.007 0.28 to 1.20 
7 1.01 2.007 0.25 to 1.18 
8 0.57 to 0.65 2.007 0.44 to 1.11 
9 0.32 to 0.41 2.007 0.40 to 1.04 
10 0.13 to 0.14 2.007 0.35 to 0.89 
ll 0.08 to 0.09 2.007 0.34 to 0.71 


New Data on Rectangular Weirs 

Original experiments on rectangular weirs were performed in 1955 in the 
Hydraulics Laboratory of the Georgia Institute of Technology at Atlanta. For 
the Georgia Tech tests the weirs were located at the downstream end of a 
steel flume 30 in. deep, 10 ft. wide, and 25 ft.long. The floor of the flume was 
level and smooth. Its width was varied by means of movable walls faced with 
1/4-in, aluminum plates. For tests on full-width weirs the wall plates were 
made to project through the notch and past the crest a distance of 6 in. Baffles 
in the forebay and curved guide walls at the upstream ends of the movable 
walls were adjusted to produce a uniform velocity distribution. Surface floats 
were used at the entrance to damp waves in the flume. 
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The basic weir bulkhead was made of 3/8-in. aluminum plate. The width 
and height of the weir notch was varied by means of lap-joined plate sections 
made of the same material. The notch edges, including the crest piece, were 
made of 1/8-in. stainless-steel plates which were fastened to the aluminum 
plates with the upstream surfaces flush. The edge pieces were machined on 
a flat-bed planer to be smooth and flat on the upstream side and beveled on 
the downstream side. The top surface of the notch edges was 1/16-in. wide 
and perpendicular to the upstream face of the weir. The upstream corners of 
the edge pieces were sharp. 

Water was supplied to the flume from a constant-head system. Discharges 
were measured in an electrically operated weighing tank at the upstream end 
of the flume. Heads were measured with a hook gage in a stilling well which 
was connected to piezometers located in both movable walls 2 in. above the 
floor and 5 ft.upstream from the weir. 

The purpose of the Georgia Tech tests was to obtain the additional experi- 
mental data necessary to describe the discharge function over the full, practi- 
cal range of the variables shown in Eq. 7. The scope of the 249 tests which 
were made at Georgia Tech is shown in Table 4. 


TABLE 4. - SCOPE OF THE GEORGIA TECH TESTS. 


P in feet b in feet b/B h in feet 
1.44 2.68 1.00 0.16 to 0,53 
1.44 1.80 0.90, 0.80, 0.60, 0.40, 0.20 0.14 to 0.58 
0.56 2.68 1.00 0.10 to 0.66 
0.56 1.80 0.90, 0.80, 0.60, 0.40, 0.20 0.10 to 0.72 
0.30 2.68 1.00 0.17 to 0.68 
0.30 1.80 0.90, 0.80, 0.60, 0.40, 0.20 0.10 to 0.76 
0.30 1.20 0.13 0.11 to 0.66 
0.30 0.577 0.06 0.28 to 0.71 
0.30 0.281 0.03 0.27 to 0.64 
0.30 0.118 0.01 0.08 to 0.65 
0.30 0.400 1.00, 0.80, 0.50 0.15 to 0.72 
0.30 0.200 1.00, 0.80, 0.50 0.10 to 0.62 


0.30 0.100 1.00, 0.80, 0.50 0.10 to 0.60 
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Verification of the Proposed Solution 


Recapitulation 

The coefficient C, in the proposed discharge equation (Eq. 11) is defined in 
Eq. 10 in terms of the effective head, he, and the effective notch width, be. It 
is also defined as a function of the geometric ratios h/P and b/B. Unlike C 
(Eqs. 4 and 7), therefore, C,. is expected to be independent of the fluid-property 
effects represented by the magnitudes of h and b in Eq. 7. 

It is assumed that C, can be determined from laboratory data if the 
measured values of h and b are adjusted as indicated in Eqs. 8 and 9. In Eq. 
8 the effective head is defined as the measured head plus a quantity kp. Simi- 
larly, in Eq. 9 the effective width is defined as the measured width plus a 
quantity kp. It is reasonable to conclude that the proposed discharge equation 
will be practicable and convenient to use only if k, and kp are constants or 
functions of not more than one independent variable. Thus, verification of the 


proposed solution is critically related to the experimental evaluation of the k 
quantities. 


Evaluation of Kh 


The quantities k, and kj, are evaluated from experimental data by a trial 
computation procedure. For a series of tests in which h is the principal vari- 
able, kj, is determined by successive approximations as the quantity which will 
“remove” the correlation between the coefficient of discharge and the head. 
The procedure is illustrated by its application to the results of the Schoder- 
Turner tests. 

Fig. 3(a) shows C as a function of h/P for two series of tests made by 
Schoder and Turner on full-width weirs (Table 2), The scatter in values of C, 
which is most apparent at lower values of h/P, is systematically related to 
the magnitude of h. If a small quantity (kj) is added to h, the value of the coef- 
ficient is reduced by an amount which is inversely related to the value of h. 
Values of the adjusted coefficient (Ce) from a trial computation, using Ky = 
0.003 ft, are shown in Fig. 3(b). A second trial, with kj, = 0.004 ft, is shown 
in Fig. 3(c). From Figs. 3(a) and 3(c) it is apparent that a constant value of 
kp = 0.004 ft. effectively eliminates the correlation between the coefficient of 
discharge and the head for Series E and F in the Schoder and Turner tests. 
From Fig. 3(c), furthermore, Ce is conveniently related to h/P by a straight- 
line equation. 

The Schoder and Turner Series E and F were selected for two reasons. 
First, they included some tests with very low values of h. For these tests, 
therefore, the effects of viscosity and surface tension which are related to the 
magnitude of h are relatively large. Second, the width of the weir used for 
the tests was large enough (4.22 ft) that the influence of kp and Ce could be 
assumed to be negligible. Proof of this conclusion is based on results of 
Georgia Tech tests which are described subsequently (See Fig. 12). 

Fig. 4(a) shows C as a function of h/P for the Bazin tests on full-width 
weirs (Table 1). Fig. 4(b) shows the same data in a plot of C, as a function of 
h/P. The value of kp, used for computing C, for these tests was 0.012 ft. This 
value was selected on the basis of trial computations as the value required to 
eliminate the correlation (evident in Fig. 4(a)) between the coefficient of dis- 
charge and the head. The effect of the width-adjustment factor, k), was as- 
sumed to be negligible for all the Bazin tests. This assumption is believed 
to have had a small effect on the computed values of Ce for only one group of 
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(a) C versus h/P. 
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tests, Series 3, for which b was 1.64 ft. From Fig. 4(b) it is apparent, again, 
that C, is a simple, linear function of h/P. 


40 


Legend 
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Fic.5-U S. BUREAU OF RECLAMATION TESTS (b/B=1.0) 


Fig. 5 shows C, as a function of h/P for the USBR tests on full-width weirs 
(Table 3). For these tests, k, was determined to be 0.003 ft. The scatter in 
the points for larger values of h/P is believed to be indicative of experimental 
errors resulting from observational difficulties related to low values of P. 
Symbols are used in Fig. 5 to indicate different values of P. The width of the 
flume was approximately 2 ft. The effect of k, was assumed to be negligible. 

Only full-width weirs (b/B = 1.0) were involved in the tests made by Schoder 
and Turner, Bazin, and the USBR. Therefore, one of the principal purposes 
of the Georgia Tech tests (Table 4) was to determine the validity of the pro- 
posed solution as applied to notch weirs. Tests made specifically to evaluate 
ky included a complete series for each of the following values of b/B: 0.20, 
0.40, 0.60, 0.80, 0.90, and 1.00. For the tests on notch weirs, the width (b) 
was 1.80 ft. For the tests on full-width weirs, b was 2.68 ft. Thus, for all of 
the kj tests the relative influence of effects related to b were negligible. Both 
h and P were varied in each series of tests in order to cover a full range of 
values of h/P. Values of h varied from 0.10 to 0.76 ft, and values of P varied 
from 0.30 to 1.44 ft. Sufficient tests were made to define the coefficient of 
discharge from h/P = 0 to h/P = 2.4. Fig. 6 shows the results of a typical 
series (b/B = 0.80). 

The successive approximations procedure was used to evaluate kj, for each 
of the series of tests described above. Remarkably, the value of k, required 
to satisfy the effective-head criteria was the same for all values of b/B and 
all values of h/P included in the Georgia Tech tests. In other words, for a 
full practical range of the significant variables a constant value of kj}, (0.003ft) 
is adequate to compensate for the fluid-property effects related to the head. 
This conclusion was subsequently verified by successfully applying the same 
value of kj, to additional tests covering values of b/B from 0.01 to 1.00 and 
values of b as small as 0.10 ft. Thus, the effective-head concept is substanti- 
ated as an important part of the proposed solution. 
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Fic. 6.—Georcia Tests (b/B=0.80) 


Differences in values of kp, and in the Ce curves determined from Schoder 
and Turner, Bazin, USBR, and Georgia Tech tests on full-width weirs are be- 
lieved to have been caused by differences in test equipment. The larger value 
of k, which characterizes the Bazin tests, for example, is believed to be an- 
other indication that Bazin’s weir plate was not effectively thin and sharp. 
This conclusion is supported by the reasoning that the effect of crest-rounding 
is similar to some of the effects of surface tension and viscosity. The results 
of the USBR tests are doubtless influenced by the fact that the crest piece of 
the weir was knife-edged. The Schoder and Turner tests may have been influ- 
enced by the fact that their flume was very short and entrance conditions were 
poor for large values of P. Most of Bazin’s tests, on ilic other hand, were 
made in flumes long enough to develop a normal, open-:}iannel velocity distri- 
bution, Thus, the results of each investigation are characterized by the equip- 
ment used in making the experiments. 

A constant value of kj, was found to be adequate for all weirs tested in the 
Georgia Tech investigation. Different values of ky were determined for full- 
width weirs tested by other investigators. Therefore, although it is indicated 
that the same value of kj, can be used for all values of b/B, the value of kj, to 


be used is also dependent on the physical characteristics of the weir and weir 
channel, 
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Evaluation of kp 

The weirs investigated by Bazin, Schoder and Turner, and the Bureau of 
Reclamation were too wide to reveal the fluid-property effects related to 
small values of b. Therefore, the results of these tests could not be used to 
evaluated kp: Some tests on weirs of relatively small width were included in 
the Georgia Tech experiments. It was the purpose of these tests to evaluate 
kp and to investigate the validity of the effective-width concept. 

In Fig. 6 the total displacement of values of Ce with respect to correspond- 
ing values of C is a measure of the combined effects of k, and kj. However, 
only a few points at low values of h/P represent heads so small that the adjust- 
ment caused by k, was appreciable. The nearly uniform displacement at high- 
er values of h/P, therefore, is the result of using a value of kp = 0.014 ft in 
the computation of C,. This value was determined from trial computations as 
the value of ky) which would eliminate the correlation between b and C, in the 
results from a special series of tests with b/B = 0.80 (not shown in Fig. 6). 
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Fic. 7.—Georcia Tests (b/B 0.13) 


Fig. 7 shows the results of special tests made to evaluate kj, for another 
value of b/B. Actually, several values of b/B are represented on the figure, 
but the results of the tests indicated that kp is essentially constant for values 
of b/B less than 0.2. In Fig. 7, as in Fig. 6, the effect of the kp adjustment is 
shown by the displacement of values of Ce with respect to corresponding 
values of C. The value of kp, used in computing Ce for this example was 0.008 
ft. 


Fig. 8 shows the results of all tests made to evaluate kp. The figure 


3.6 
Note: ky =+ 0.003ft; kp =+0.008f 
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indicates that kj, is a function of b/B, that it increases with increasing values 
of b/B until it reaches a maximum positive value of 0.014 ft at b/B = 0.80, and 
that it decreases rapidly thereafter, reaching a minimum value of -0.003 ft at 
b/B = 1.0. The authors acknowledge some uncertainty regarding the generality 
of the curve in Fig. 8. It is believed that the value of k, obtained from the 
Georgia Tech tests may have been influenced, to an unknown extent, by the 
physical characteristics of the equipment used. Tests made especially to 
evaluate kj involved very small values of b and correspondingly small values 
of B and Q. Thus, unavoidable errors of measurement may account for com- 
paratively large relative errors in the computed results. Nevertheless, the 
adequacy of the effective-width concept was demonstrated by the successful 
elimination of b as an independent variable in the experimentally determined 
relationship between Cg, h/P, and b/B. 

A remarkable example of the experimental evidence which substantiates 
both the effective-head and the effective-width concepts is shown in Fig. 9. In 
this figure a single, straight line satisfactorily represents the relationships 
between Ce and h/P for full-width weirs involving values of b from 0.10 to 
2.68 ft and values of h from 0.10 to 0.72 ft. 


Evaluation of 


The results of the tests from all sources indicate that Ce, the “effective” 
coefficient of discharge, is a function of b/B and h/P alone. Thus, in confir- 
mation of Eq. 10, Ce is independent of the influences of viscosity and surface 
tension. Furthermore, the results show that the relationship between C, and 
h/P for all values of b/B is conveniently defined by a family of straight lines. 

Fig. 10 shows a summary of the results of the Georgia Tech tests. The 
information on this figure, in combination with values of kp, and kp (Fig. 8) de- 
termined from the same tests, is the information needed to evaluate the 
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proposed discharge equation. Thus, Eq. 11, complemented by the effective- 
head and effective-width concepts, is confirmed as a comprehensive equation 
of discharge for all rectangular, thin-plate weirs. 


Summary 


Experimental Results 

The proposed discharge equation (Eq. 11) and the concepts of effective head 
and effective width (Eqs. 8 and 9) have been substantiated by experiments from 
various sources covering a wide range of test conditions. It has been demon- 
strated that C, is a linear function of h/P for each value of b/B, that kj is a 
constant for all values of b/B and h/P, and that kj is a function of b/B. 

From the results of tests made at the Georgia Institute of Technology, 
values of Cg, ky, and kp have been determined for water at ordinary tempera- 
tures. Fig. 10 shows Cg for a full range of values of b/B, with h/P from 0 to 
2.0. From the information contained in this figure, simple equations for C, 
can be written if desired. The constant value of kj determined from the 
Georgia Tech tests is 0.003 ft. Values of kp are shown in Fig. 8. 

Comparisons of values of C, and kp obtained from tests made by different 
investigators of full-width weirs confirm the general conclusion that weir dis- 
charge is critically influenced by the physical characteristics of the experi- 
mental setup. Especially critical are the weir-notch edges and the channel 
properties which control the velocity distribution and turbulence upstream 
from the weir. 

Neither superior equipment nor superior technique is claimed for the 
Georgia Tech tests. However, it is believed that the over-all quality of these 
data is as good as that of any other data available. Furthermore, this is the 
only investigation known to cover a full, practical range of the significant vari- 
ables. Therefore, as a reasonably accurate solution for the comprehensive 
discharge characteristics of rectangular weirs with effectively sharp, thin 
crests, the authors recommend the use of the proposed discharge equation and 
the coefficients derived from the Georgia Tech experiments. 


Effect of Neglecting kp and kp 

It has been established that the relative influence of viscosity and surface 
tension diminishes as both the head and the notch width becomes larger. In 
the proposed discharge equation the fluid-property effects are represented by 
the head- and width-adjustment factors, kp and kp. Thus, the relative influ- 
ence of ky decreases as the head increases in magnitude, and the relative 
influence of k, decreases as the width increases in magnitude. It follows that 
for “large” heads on “large” weirs the influence of the adjustment factors is 
negligible. The experimental results can be used to evaluate the influence of 
kp and kp in relation to the magnitudes of h and b. 

A convenient approximation for the discharge can be obtained from Eq. 11, 
using Cg from Fig. 10, with the unadjusted (measured) values of h and b, The 
relative accuracy of this approximation is a significant measure of the influ- 
ence of k), and kp. It is also a measure of the influence of the combined effects 
of viscosity and surface tension. Thus, if the absolute error (AQ) due to 
neglecting kp and kp is 


= 3/2 C.bh3/2 (12) 
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the relative error is 


Q h 3/2 

Q (13) 

The influence of k; alone can be demonstrated by assuming that the notch width 
is very large, whence b= be in Eq. 13. Then, from Eqs. 8 and 13, 


3/2 


which is shown in Fig. 11 for kp = 0.003 ft. Similarly, the separate influence 
of k, can be demonstrated by assuming that h = he. Thus, from Eqs. 9 and 
13, 


(14b) 


which is shown in Fig. 12 for values of k, from Fig. 8. Figs. 11 and 12 are 
recommended for use in connection with Eq. 11 as a means of determining 


when (depending on the accuracy desired) the influence of the fluid properties 
may be neglected. 


Upper Limit of h/P 

Boss has shown that critical depth will occur in the uniform flow up- 
stream from a weir with h/P greater than about 5. For values of h/P greater 
than 5, therefore, the weir is not a control, and the ordinary weir discharge 
equations are not applicable. A lower, practical limit on h/P is imposed by 
observational difficulties and measurement errors. For precise measure - 
ments the writers recommend that h/P be limited to values less than 2.0. 


Influence of Weir Height 

The weir height, P, was acknowledged as a possible third significant length 
parameter in the discharge function. There is disagreement in the literature, 
however, regarding experimental evidence of a systematic correlation between 
P and the coefficient of discharge. Bazin 9) warned against the use of his dis- 
charge formulas for a “very low weir. . . which should always be avoided.” 
Rehbock'‘10 acknowledged “some disturbance, which is not yet clearly under- 
stood (and because of which) the coefficient C is not exactly the same for any 
two sharp-crested weirs of different absolute heights, P, even if the ratios, 
h/P, are identical.” Weirs as low as 0.08 ft were used in the USBR tests in 
order to achieve the higher values of h/P. The results of these tests show a 
definite correlation between Ce and P for values of h/P greater than 2 (not in- 
cluded on Fig. 5). The Georgia Tech tests, however, failed to show a similar 
correlation for values of P as low as 0.30 ft and values of h/P as high as 2.4. 

It is generally acknowledged that observational difficulties and errors in 
measurement may account for appreciable errors in the results of experiments 
involving small values of P and large values of h/P. For this reason alone, 
weirs lower than 0.3 ft are not recommended for precise measurement of dis- 
charge. Nevertheless, although questions regarding the influence of weir 
heights less than 0.3 ft may be academic, additional experiments are needed 
to resolve the contradictions to be found in the literature. 
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Comparison with Other Experiments and Formulas 


Some Classic Formulas 

One equation is frequently described as the basic equation of discharge for 
all weirs. It is founded on an assumed analogy between the weir and the two- 
dimensional orifice. In its derivation an approximate velocity equation is inte- 
grated over the approximate area limits of the nappe in the plane of the weir. 
The effects of viscosity and surface tension are ignored. The result, in the 
form usually ascribed to Weisbach, is 


(15) 


in which C, is a coefficient of contraction and V, is the average velocity in the 
approach channel. 

Some writers have erroneously implied that Eq. 15 is a comprehensive and 
theoretically correct equation of discharge. It is claimed, for example, that 
the velocity-head terms account for the influence of approach-channel geome- 
try. From this viewpoint h/P is eliminated as an independent variable, and 
Ce; is a function of b/B alone. These conclusions are not substantiated by ex- 
periment, 

The inadequacy of Eq. 15 can be attributed partly to the fact that the effects 
of viscosity and surface tension are ignored in its derivation. However, ex- 
periments show conclusively that these effects are appreciable only when h or 
b, or both, are relatively small. Eq. 15 is deficient largely because its deri- 
vation was based on several erroneous assumptions. 

Unjustifiable alterations of Eq. 15 have been made for the purpose of im- 
proving agreement between the equation and experimental data. Thus, it is 
quite common to omit the second velocity-head term in the brackets, and, 
subsequently, to attribute great significance to the “total head” term which 


remains. The result is an equation such as that proposed by Hamilton Smith, 
Jr., in 1884,(12) 


(16) 


in which both Cg, and ky are coefficients determined from experiment. Smith 
described the coefficient kj as a number which represents “the portion of the 
kinetic energy [which] produces a useful effect at the opening of discharge.” 
He recognized that k; is a variable which depends on the velocity distribution 
in the approach channel. Nevertheless, he recommended a constant value of 
1.33 for all full-width weirs and 1.40 for all “fully contracted” weirs. He 
described C, as a discharge coefficient, and he concluded that it was a function 
of the head, the notch width, and the width-contraction ratio (b/B). Actually, 
both k; and C, are functions also of the ratio h/P, and they are critically de- 
pendent on the physical characteristics of the weir and the approach channel. 
Therefore, Eq. 16 is little more than a basis for the experimental definition 
of the discharge function. It is obviously not the most convenient formula for 
this purpose, nor is it the comprehensive solution which is needed. 

Hamilton Smith’s exhaustive study and detailed exposition of the influence 
of the velocity of approach was subsequently ignored or misinterpreted by 
others. Some writers have inferred that his kj coefficient can be associated 


— V12.3/2 
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with the kinetic-energy coefficient (a1). Accordingly, they have perpetuated 
the misconception that the coefficient of discharge is a constant when the 
quantity within the brackets (Eq. 16) is evaluated as the true total energy head 
in the approach channel [h + (V,2/2g)]. 

The coefficient @ is defined as the ratio of the true average velocity head 
in a cross section to the velocity head computed on the basis of the average 
velocity in the section. It is a dimensionless measure of the degree of velocity 
nonuniformity, but it does not account ane for the shape of the velocity- 
distribution curve. Smith, Schoder and Turner, ( ) and others who have studied 
the influence of the velocity distribution in the approach channel have shown 
conclusively that the pattern of nonuniformity is more critical than the degree 
of nonuniformity. Thus, experiments performed in a large variety of weir 
channels have failed to reveal a satisfactory correlation between a, and kj. 

Most of the defects of Eq. 16 are common also to a large number of 
equations in the form suggested by Weisbach in 1844, 


h 
3 2elk, + k3(5 =) Ipn?/2 (17) 


in which Cy, Kg and kg are experimentally determined coefficients which ac- 
count for the influence of the fluid properties, the velocity distribution, and 
the physical characteristics of the weir and approach channel. The principal 
advantage of Eq. 17 is that the velocity-head term in Eq. 16 is replaced by a 
term which involves the h/P ratio. The effect of this modification is to avoid 
a successive-approximations solution for discharge. Equations of this gener- 
al variety include some of the most widely used formulas for full-width weirs 
(b/B = 1.0). Among these are the Bazin formula, (5) 


h 3/2 


2 
(18) 


Q = (3.25 + pli + 0.55¢ 


the Rehbock formula of 1912, (13) 


Q = (3.24 + 0.43 ph3/2 (19) 
and the S.I.A. (Swiss Society of Engineers and Architects) formula,(14) 
0.011 3/2 
Q = (3.29 + + bh (20) 


The complex expressions which are the coefficients of the quantity bh3/2 in 
each of these equations are equivalent to C in the simple equation of discharge 
(Eq. 4). The numerical quantities are experimentally derived coefficients 
which account for the influence of the fluid properties as well as the weir and 
channel characteristics. Specifically, a term involving and having the di- 
mensions of h is now recognized as a means of compensating for the combined 
effects of viscosity and surface tension. Because it is a dimensional quantity, 
the h term restricts the application of Eqs. 18, 19, and 20 to the fluid used 

for the experiments (water at ordinary temperatures). However, it is apparent 
that the influence of the head term is negligible except when h is relatively 
small. 


~- — formula for full-width weirs was proposed by Rehbock in 
1929(15); 


| 
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Q = (3.22 + 0.44 8) pn?/? (21) 


This was the first formula known to involve the effective head (hg = h + 0.004 
ft). The quantity 0.004 ft, corresponding to k, in Eq. 8, was apparently derived 
by Rehbock from an analysis of his own experiments, 

Rehbock did not associate the h term in Eq. 19 or the effective-head term 
in Eq. 21 with the influence of the fluid properties. He described both terms 
as a means of accounting for a disturbance related to the magnitude of P. He 
hesitated to endorse Prandtl’s suggestion that they were related to the effect 
of capillarity, (10) 

Various restrictions were placed on Eqs. 18, 19, 20, and 21 by their auth- 
ors. Bazin(9) specified that “the arrangement of our standard weir [must be] 
exactly reproduced.” Rehbock(15) required that the head be large enough to 
ensure a free nappe, but he claimed a high degree of accuracy for both of his 
formulas when applied to full-width weirs of all sizes. The S.IL.A., in its 
Standards for Water Measurement, 14) restricts the application of Eq. 20 to 
values of P not less than 1 ft, values of h between 0.08 and 2.6 ft, and values 
of h/P not greater than 1.0. 


One of the new formulas developed for rectangular-notch weirs (b/B < 1.0) 
is that proposed by the S.1.A.:(14) 


0.063 - 0.053(2)4 
B h + 0.005 


22 
(1 + 0.5(2)* pn3/2 (22) 
B P+h 


The application of this formula is restricted to values of P not less than 1 ft, 
values of h between 0.08(B/b) and 2.6 ft, values of h/P not greater than 1.0, 
and values of b/B not less than 0.3. Eq. 22 is widely quoted in European 
technical literature, but correspondence and library research failed to reveal 
the source of the data on which it was based. The formula is virtually unknown 
in this country. 

In American practice, approximate formulas for notch weirs are adapted 


from formulas for full-width weirs Mm means of the contracted-width equation 
ascribed to James B. Francis:(16, 17 


bo =b-O0.2h 


(23) 


In this equation b, is the width of the contracted stream, b is the gross width 
of the notch, and the quantity 9.2 is a coefficient recommended by Francis for 
“fully contracted” rectanguJ.r-notch weirs. Thus, from Eqs. 4 and 23, the 
Francis discharge formula for notch weirs is 


Q = C, (b - 0.2 bh) (24) 


in which C, is the coefficient of discharge in any of the formulas for full-width 
weirs (b/B = 1.0). 


Hamilton Smith, Jr., after reviewing the experiments of Francis, Fteley 


| 
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and Stearns, and others, described the requirements for “full contraction” as 
follows:(17) The weir height (P) should be not less than twice the head (h) or 
not less than 1 ft; the width from the side of the approach channel to the end 
of the notch should be not less than twice the head or twice the width of the 
notch (b), Thus, Smith’s recommendations included the requirement that h/P 
be not greater than 0.5 and that b/B be not greater than 0.2. Francis(17) add- 
ed the requirement that b/h be not less than 3.0. It is recognized now that the 
concept of “full contraction” is inconsistent with the true significance of b/B 
as an independent variable. 


Comparison of Experiments and Formulas for Full-Width Weirs 

The experiments of Schoder and Turner, Bazin, and the Bureau of Recla- 
mation, as well as the Georgia Tech experiments, were used to confirm the 
method of analysis and the discharge equation proposed in this report. Conse- 
quently, values of Cp and ky, were determined, and these values can be used 
to compare the results of the tests made by the different investigators. Com- 
parable values of k, could not be determined because only the Georgia Tech 
experiments included tests on weirs narrow enough to reveal the fluid-property 
effects associated with small values of b. 

From the results of the Georgia Tech tests on full-width weirs, kj, = 0.003 
ft, k, = -0.003 ft, and C, is given by the simple equation, 


= h 
Ce = 3.22 + 0.40 > (25) 


From the Schoder and Turner tests (Series E and F) in Fig. 3(c), kp, = 0.004 
ft, and 


h 
Ce = 3.21 + 0.45 5 (26) 


The results of the Bazin tests on full-width weirs are shown in Fig. 4. The 
value of k}, required for the Bazin tests is considerably larger than that de- 
termined from the other experiments. From Fig. 4(a), kj, = 0.012 ft, and 


Cc 


h 
= 3.25 + 0.44 


(27) 


The results of the USBR tests (Fig. 5), in comparison with the others, are 
characterized by smaller values of C, at the larger values of h/P. For the 
USBR tests, kp = 0.003 ft, and 


= 3.25 + 0.35 (28) 


Rehbock’s 1929 formula is also based on the effective-head concept. 
Furthermore, his equation for C, is in the linear form of Eqs. 25 to 28, in- 


clusive. Thus, for comparison with those equations, Rehbock’s 1929 formula 
gives 


h 
C. = 3.22 + 0.44 = 
(29) 


with kj}, = 0.004 ft. Rehbock did not recognize the existence of k,,, although 


| 
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some writers have indicated that his experiments were made in a very narrow 
flume. Therefore, it is possible that Eq. 29 would be modified if kp were taken 
into account. This conclusion could not be verified because Rehbock’s experi- 
mental data were not available. 


4.2 


Value of Ce 


4.0 
3.8 
Georgia Tech tests (Eq.25) 
USBR tests (Eq.28) 
3.4 
= 
10) 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8 2.0 
Value of h/P 


Fic. 13.- COMPARISON OF EXPERIMENTAL DaTaA (b/B=1.0) 


Eqs. 25 to 29, inclusive, are compared in Fig. 13. The Schoder and Turner, 
Bazin, and Rehbock curves agree reasonably well. The curve determined from 
the Georgia Tech tests lies below these curves and above the USBR curve. It 
was suggested previously that the differences in the results of the experiments 
shown in Fig. 13 were caused by differences in the test equipment. Thus, the 
disti.ctive character of the C, curve from the USBR tests can be associated 
with the fact that the crest piece for these tests was knife-edged. Similarly, 
the larger value of kj, determined from Bazin’s tests can be related to the 
conclusion that his crest was not truly thin and sharp. Such details as these, 
frequently unknown or uncontrollable, are believed to be the cause of the vari- 
ance shown in Fig. 13. Nevertheless, it is noteworthy that the equations for 
C, derived from various sources are all simple, linear functions of the h/P 
ratio. 

Fig. 14 shows a comparison of the results of the Georgia Tech tests with 
some of the most widely used formulas for full-width weirs. The comparison 
is made for two values of h and a full range of values of h/P. The Bazin 
formula (Eq. 18), the Rehbock 1912 formula (Eq. 19), and the S.I.A. formula 
(Eq. 20) are in the non-linear form suggested by Weisbach. For these formu- 
las C was computed as Q/bh3/2 from Eqs. 18, 19, and 20. For the Georgia 


Tech tests and the Rehbock 1929 formula, C was computed from the relation- 
ship 


Sebelh + ky)3/2 + 3/2 


— Ca 
pho/2 pho/2 h ) (30) 


|| 
Rehbock (1929) formula (Eq.29) Jf 
| 
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in which b is assumed to be so large that the influence of k, can be neglected, 
whence b=be. Values of Ce in Eq. 30 were computed from Eq. 25 for the 
Georgia Tech tests and from Eq. 29 for the Rehbock 1929 formula. 

Fig. 14(a) shows a comparison of the different formulas for h = 0.2 ft. For 
this condition the influence of the fluid properties and effects of differences in 
the experimental equipment are relatively large. Fig. 14(b) shows a compari- 
son of the formulas for h equal to 1.5 ft. For this condition the effects related 
to small values of h are negligible (See Fig. 11). The curves in Fig. 14(b) 
show reasonably good agreement. 

In general, the comparisons of experiments and formulas shown in Figs. 13 
and 14 cast doubt on the adequacy of some of the widely used formulas and on 
the reproducibility of the test conditions from which they were derived. In 
particular, the relatively complex Weisbach-type formulas are not justified 
as comprehensive discharge equations. This conclusion is emphatically demon- 
strated in Fig. 4. From the family of curves shown in Fig. 4(a), it is apparent 
that the Bazin formula compares unfavorably with the results of the Bazin 
tests for all but the intermediate values of h included in the experiments. In 
Fig. 4(b), however, a single, straight-line curve fits the data reasonably well 
for all values of h and h/P, Thus, the simple, linear equation for Cg, comple- 
mented by the effective-head concept, is substantiated as an improved method 
of defining the discharge function. 


Comparison of Formulas for Notch-Weirs 

Fig. 10, based on the Georgia Tech tests, shows the relationship between 
Ce, i/P, and b/B for notch weirs. A comparable family of curves can be 
plotted from the S.I.A. formula for notch weirs (Eq. 22) if the term involving 
the absolute head is omitted. For relatively large values of h (and b), the 
influence of the h term is negligible, and values of C from Eq. 22 can be com- 
pared with values of C, from Fig. 10. The comparison is shown in Fig. 15. 
The maximum difference between corresponding values of Cg shown in Fig. 15 
is about 4 percent. Howev TS within the limits of applicability specified by 
the S.A. (b/B 2 0.3, h/P , the maximum difference between the two fami- 
lies of curves is about 2 Arcisose| It is emphasized that Fig. 15 does not in- 
clude a comparison of the different methods of compensating for the fluid- 
property effects. However, it is believed that the advantages of the effective- 
head concept were clearly demonstrated in the previous comparisons of formu- 
las and experiments for full-width weirs. 

The results of the Georgia Tech tests on notch weirs (Fig. 10) indicate that 
Ce is a function of b/B and h/P. In the Francis contracted-width formula 
(Eq. 23), however, a correction for “end contractions” is assumed to be a 
function of h alone. Because of this fundamental difference, the Francis dis- 
charge formula (Eq. 24) cannot be compared with the Georgia Tech tests or 
the S.I.A. formula on a graph such as Fig. 15. Nevertheless, it is possible to 
show that the Francis formula is not substantiated by the curves shown in Fig. 
10. 

Eq. 23 can be written in the form 


(31) 


in which, according to Francis, K is a constant equal to 0.2 for “fully contract- 
ed” rectangular-notch weirs. From Eqs. 11, 24, and 31, 


4. 
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Q = Cebghe?/? = Co - 


in which C, is the coefficient of discharge for notch weirs and C,, is the corre- 
sponding coefficient for full-width weirs (b/B = 1.0). If the fluid-property ef- 
fects related to small heads and widths are neglected, h = he and b = be. 
Thus, from Eq. 32, 
K = b 
(33) 


Fig. 16 shows K as a function of h/P and b/h. Values of K shown in this figure 
were computed from Eq. 33, using values of C, and Cg, from Fig. 10. It is 
apparent that the constant value of K recommended by Francis is not substanti- 
ated, even within the limits specified by Smith and Francis for “fully contract- 
ed” weirs. 


Value of K 


— K=0.2 (Francis) 
c 


Curves based on 
Georgia Tech Tests 


5 10 15 20 
Value of b/h 


Fic. |6- ComPARISON WiTH FRANCIS FoRMULA FOR NoTcH WEIRS 


The comparison shown in Fig. 15 indicates that the S.I.A. formula for notch 
weirs agrees reasonably well with the results of the Georgia Tech tests. How- 
ever, the complexity of the formula is a deterrent to its use, and the limits of 
applicability specified by the S.I.A. prevent its being a comprehensive solution 
for a full, practical variety of notch weirs. Furthermore, the effective-head 


| 
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and effective-width concepts which are complementary to the simple C, 
equations derived from the Georgia Tech tests are believed to be superior to 


any other method of accounting for the influence of viscosity and surface 
tension. 


CONC LUSIONS 


The flow pattern for rectangular, thin-plate weirs is not subject to com- 
plete mathematical analysis. Consequently, an analytical solution for the dis- 
charge characteristics has not been developed. A comprehensive solution 
based on dimensional analysis and experiment is presented in this report. It 
provides a simple, direct solution for the discharge and a convenient method 
of compensating for the influence of viscosity and surface tension. It is sub- 
stantiated by experimental data from various sources. 

Widely used formulas adapted from the classic Weisbach equation are not 
justified as comprehensive solutions. In fact, they do not satisfy the need for 
a simple equation which can be used to correlate experimental data covering 
a full range of the significant variables. Unjustified emphasis has been placed 
on the significance of the orifice analogy in the derivation of weir equations. 
Furthermore, there is no logical basis for the assumption that formulas involv- 
ing the true total-energy head are independent of the h/P ratio. The Francis 
contracted-width formula is fundamentally inadequate as a means of adapting 
formulas for full-width weirs to notch weirs. 

Weir discharge is critically influenced by the physical characteristics of 
the weir and weir channel. It is especially dependent on features which control 
the velocity distribution in the approach channel. For this reason the results 
of experiments made by different, capable investigators do not agree, and 
formulas based on a particular set of data reflect the individual characteristics 
of those data. For this reason, too, a truly reproducible, standard measuring 
weir and a precise, universal discharge formula are impractical. 
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SYNOPSIS 


Air Binding is a potentially serious cause of excess head loss in water 
pipelines flowing under less than atmospheric pressure. This paper dis- 
cusses the phenomenon of air binding and methods of control, with emphasis 
on the accumulation of air in sloping pipes. The specific installations 
described are circulating water systems of steam-electric stations, but the 
principles discussed are applicable to any water pumping station with similar 
profile. 


INTRODUCTION 


Air binding in a pipeline is the trapping of air within the water passage in 
a manner which prevents the pipe from flowing full. In poorly laid out or 
poorly vented pipelines air binding may be the source of head losses which 
can seriously reduce the pipeline capacity. 

The accumulation of air is particularly severe in a pipeline flowing as a 
siphon under less than atmospheric pressure. Such a condition may occur, 
for example, where the pipe passes over a hill. 

Important water pumping installations which normally utilize a siphon are 
the circulating water systems of steam-electric generating stations. These 
pumping plants may move as much as 150,000 gallons per minute through the 
condenser of each large generating unit. Large-diameter circulating water 
pipelines, several hundred feet long, are often designed to flow under a 
vacuum between minus 20 and minus 30 feet of water. The phenomenon of air 
binding is illustrated in this paper by two examples of circulating water sys- 
tems which were seriously obstructed by air. Excess head losses due to 
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entrapped air were 15 feet and 19 feet for systems designed originally to 
pump against normal pressures of 43 and 36 feet. 

It will become apparent by study of these examples that the air binding 
problem may be applicable to any other pipeline with similar profile. 


Basic Elements of the Air Binding Problem 


Air accumulation in a pumping system may originate from several 
sources: air exists in the empty pipeline before it is filled, air is entrained 
in the water carried into the system and, at sections under negative pressure, 
it can leak in at joints and fittings. This air will be present in varying con- 
centrations from very small, widely dispersed bubbles to large pockets. 

The air is carried along the pipeline by the moving water. However, where 
the pipe slopes downward from a high point, the bubbles tend to collect along 
the top of the pipe. When the buoyancy of the air accumulation becomes 
greater than the water force pushing it down the pipe the pipeline becomes 
air bound. 

The most severe head losses encountered in circulating water systems 
tlowing under vacuum are the result of large air pockets lying along down- 
ward slopes and creating a condition similar to open channel flow. The ener- 
gy gradient thus roughly parallels the actual slope of the pipeline instead of 
representing only the friction losses in a pipe flowing full. Figure lisa 
diagrammatic view of the phenomenon of airbound slopes. The air pocket 
terminates in a hydraulic jump which absorbs further energy. 

Generally speaking, for a given water velocity and pipe slope air binding 
difficulties increase with the pipe size. For example, an 8-foot-per-second 
velocity may be sufficient to carry all accumulating air out of a 12-inch pipe 
on a given slope, but is not sufficient to move such air through an 84-inch 
pipe on the same slope. For large-size pipe there is very little published in- 
formation relating the factors of air entrainment, air release, water velocity, 
pipe size, pipe slope and natural air removal. 

The author has been concerned particularly with conduits of such a size 
that water velocities up to approximately 10 feet per second have not been 
sufficient to move air down a slope. Pipes of five feet or more in diameter 
fall into this class. The discussion is, of course, applicable to smaller con- 
duit where water velocities are low. 

Experience indicates that where a large diameter pipe normally flows un- 
der vacuum there will be air binding in any non vented section which slopes 
downstream from the high point in the system. Accumulating air must be 
continuously removed or the efficiency of the system will be impaired. Ex- 
ceptions to this general rule have been at slopes immediately following sharp 
horizontal bends; the turbulence induced by the horizontal bend may break up 
large air pockets into more easily transported smaller bubbles. 


Self-Priming Vertical Drops 


When a change in grade is made with a vertical or nearly vertical drop the 
possibilities of air binding are much reduced. Such a vertical drop is gener- 
ally “self-priming.” The turbulence created by the abrupt change in water 
direction breaks up and disperses large pockets of air. The smaller bubbles 
are then carried away by the moving water. See Figure 2. 
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Vertical Versus Sloping Pipelines 


It is apparent from the comments above that air evacuation equipment 
must be used unless the change in grade is made with a vertical drop. This 
equipment, discussed briefly below, is expensive and requires constant main- 
tenance. However, vertical drops have disadvantages also, and in many cases 
may be more costly to install and operate than a continuously evacuated slop- 
ing pipeline. 

First, the turbulence which entrains and carries away the air at a vertical 
bend will also create a permanent loss of head which raises the cost of pump- 
ing water. This loss for one vertical drop with two reverse 90° bends may 
amount to 1 foot of head and represent an annual excess pump operating cost 
of several hundred dollars and a capitalized cost of several thousand dollars. 
One company solved a potential air binding problem in one of their older 
plants by traversing a downgrade with a series of three vertical drops, a total 
of six 90° bends. At a large plant the cost of head lost in such an arrange- 
ment would be considerably more than the cost for the installation and main- 
tenance of an air removal system on a sloping pipeline. 

A second serious objection to the vertical drop is the thrust produced at 
90° bends. Massive thrust blocks must be provided to anchor such a bend 
securely in the event of water hammer surges. At one large power station 
the thrust block required at a vertical drop in twin 84-inch conduits contained 
over 400 cubic yards of concrete and required complicated formwork. A new 
circulating water system at this same station has a sloping conduit equipped 
with an air removal system at a fraction of the cost of the thrust block. 


Methods of Air Removal 


Air taps at high points are standard practice on pipelines, but it is appar- 
ent from the previous discussion and a knowledge of laboratory observation 
that the air carried by moving water will be found along the downstream 
slope. Air removal systems on the conduits discussed below were not effec- 
tive when they were connected to the high point taps alone. Accordingly, a 
system of taps along the downstream slope was provided. These taps were 
connected to a common header which in turn was connected to the air removal 
system. 

Where a vertical drop is used a single air removal connection should be 
placed at the high point at the outside edge of the bend. Generally, this tap 
will not be required during operation unless water velocities are very low. 

Where an airbound pipe flows under vacuum some type of air evacuator 
will be required to remove the air accumulations. Such an installation may 
consist of one or more vacuum pumps, a control system, power supply and, 
for water-sealed pumps, a continuous supply of clean water (or a closed- 
circuit sealing system using some other fluid). The author’s firm has recent- 
ly installed a water jet evacuation system which is considerably less complex 
and costly than vacuum pumps. Operating water for the jets is supplied from 
the power station elevated water tank, but could be supplied by booster pump 
taking suction from the pipeline itself. 

The primary consideration in sizing air removal equipment is to provide 
enough capacity to maintain the line free of air during operation. A second- 
ary consideration is the rate of initial priming of the system. 
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The determination of capacity requirements is a difficult matter. Uncer- 
tainty is caused by several factors. A first unknown is the quantity of air 
which may leak in through joints and fittings in long conduits under vacuum. 
Also, experimental data for these large-size pipes is not sufficient to permit 
determination of the amount of entering air which will be carried out natural- 
ly for given velocities, slopes and pipe sizes. 

It may be desirable to install one fairly large vacuum pump for initial 
priming, then a smaller pump in parallel to maintain the prime. The smaller 
one can be selected and installed after actual operating experience with the 
large vacuum pump has indicated the air binding characteristics of the pipe- 
line. 

The details of two typical circulating water systems are described below. 
Both of these are older stations at which air binding problems had to be cor- 
rected. For all new stations the original design should provide for air re- 
moval if the profile indicates a potential problem. 


Station A - 132,000-KW Power Plant Installation on a Large River 


Station “A” has a 1,750-foot condenser discharge conduit, all of which 
should normally operate under a partial vacuum. The majority of the pipe is 
78-inch diameter, the design velocity is 7.3 feet per second. The profile of 
this conduit is shown on Figure 3. The head recovery in this discharge 
siphon line was very poor and caused considerable concern because of the 
loss in water pumping capacity. 

The vacuum measured at the condenser discharge should be minus 22.9 
feet of water for river elevation 718. The actual measurement at the north 
condenser was minus 7.9 feet, an unexplained loss of 15 feet. It was deter- 
mined that with river water levels at 719 or above air accumulated along the 
57:1 sloping section shown on Figure 3, giving a head loss at water level 719 
of about 7.5 feet. This represented the approximate difference in elevation 
between the top and bottom of the slope, and indicated that the pressure 
gradient was parallel with the pipe slope as could be expected for open chan- 
nel flow down a similar incline. At water levels under 719 the air pocket also 
extended along the entire top of the horizontal section and into the 45° sloping 
condenser discharge pipe, thereby reducing the siphon an additional 7.5 feet. 

When the river level was at 709, flow through the pipeline was recorded as 
follows: 


Actual Flow 


Design Flow in Pipeline Total Loss 

With Proper Venting No Venting of Capacity 

Two Pumps 70,000 gpm 58,000 gpm 12,000 gpm 
Three Pumps 89,000 gpm 74,000 gpm 15,000 gpm 


The points of air evacuation had been installed in accordance with original 
design, but proved to be inefficient. A 95-cfm, 15-inch Hg vacuum pump was 
located at Point A where the discharge line passed the intake structure. It 
worked continuously and was of insufficient capacity. Two-inch air-powered 
ejectors at the condensers were also ineffective. 

As a result of investigation, a 200-cfm vacuum pump was installed at 
Point B on Figure 3. Additional taps were also installed on either side of 
Point B, and on the upper end of the 45° condenser piping. These all take 
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suction from the same pump. The 95-cftm pump at the intake structure was 
also replaced with a 200-cfm pump. 
This arrangement was effective in eliminating the air binding problem. 


Station “B” - 166,000-KW Power Plant Installation on a Large River 


The circulating water system of this station has three pronounced high 
points, the usual one at the condensers, and two where the intake and dis- 
charge conduits pass over a flood protection levee. The principal airbound 
section is 66-inch diameter, design velocity is 10.7 feet per second. An ap- 
proximate profile is shown on Figure 4. 

The pressure at the crest of the discharge levee crossing was about minus 
3.5 feet without venting. The design pressure was minus 23 feet, a difference 
of 19.5 feet representing excess lost head. A vacuum pump taking suction 
from a Single tap at the crest had been able to maintain a maximum vacuum 
of slightly less than 15 feet, even though the gage glass at the air receiving 
tank showed full prime at this point. Taps were therefore drilled at three 
additional points on the downstream slope of the pipe and connected to the 
vacuum pump. Maximum design siphon of minus 23 feet was obtained within 
less than one hour and could be maintained indefinitely with a 6-cfm pump. 
The air had been trapped on the slope beyond the levee crest as could be ex- 
pected from previous discussion in this paper. 

With the increased vacuum provided by the air evacuation system four- 
pump flow in the system was increased from 115,000 gpm to 134,000 gpm. 
This increased flow was so great that the pump impellers had to be rede- 
signed to prevent excessive impeller pitting due to cavitation. 

The amount of air to be removed from this system was very small since 
little leakage of air could take place. The only section under vacuum was 
that section of welded steel pipe over the levee. 

When only one or two of the four circulating water pumps were operating 
at low water conditions a slight vacuum also occurred at the intake levee 
crest and a Similar air-binding condition arose. Vacuum pump connections 
were run to this point also. A completely separate parallel circulating water 
system has recently been installed for Units 3 and 4 of this power plant. A 
bank of vacuum pumps with interconnected controls and suction lines draws 
air from both circulating water systems at the intake and discharge levees. 
The maximum length of air suction line is about 600 feet. 


CONCLUSION 


These two examples clearly indicate the extent of capacity loss which may 
be experienced in pipelines flowing downhill under a vacuum. This phenome- 
non probably occurs much more frequently than records show. In many 
steam plant circulating water systems overdesigned pumps have made up for 
lack of hydraulic efficiency. Air binding may well have gone unnoticed at 
Station B if the extent of head loss had been a little less severe. At Station A, 
on the other hand, air binding threatened to reduce the efficiency of the entire 
power station as summer weather brought a greater demand for cooling 
water. 

There is a particular need for experimental air-binding data applying to 
large diameter pipelines. Observations at Station “B” for several different 
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pumping conditions have indicated that the pipe may be self-evacuating at 
velocities slightly above the design. Empirical formulas based on enough 
carefully evaluated field data may well permit the confident design of a self- 
evacuating sloping pipeline flowing under vacuum. 
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FLOW THROUGH CIRCULAR WEIRS 


J. C. Stevens,! M. ASCE 
(Proc. Paper 1455) 


SYNOPSIS 


There are certain advantages inherent in the use of the circular weir. How- 
ever it has not come into general usage in this country. It is more extensively 
used in Europe where there has been developed a rigid, accurate formula by 
which the flow can be determined. 

There is presented herein a formula for the theoretical flow through sharp 
crested circular weirs which involves the use of elliptic functions. Most of 
the European treatises on such weirs are not generally available now. Tables 
appended give (1) a summary of available experimental data with discharge 
coefficients, (2) the flow through weirs of various diameters based on these 
coefficients and (3) values of the elliptic functions involved by which the flow 
through a circular weir of any diameter may be found. 


INTRODUCTION 


Weirs with circular crests have the advantage that the crest can be turned 
and beveled with precision in a lathe. Moreover the weir crest does not have 
to be leveled, hence the point of zero flow is readily determined. 

When wells are tested by pumping the pump discharge is usually measured 
by an orifice fixed to the end of the discharge pipe. The flow is determined 
by the orifice formula in which the down-stream head is constant—that of one 
atmosphere. When the pipe is not running full the orifice becomes a circular 
weir and it is often important to know the flow under such conditions. 

The problem of circular weir discharge first thrust itself upon the author’s 
attention in connection with a rather novel fishway where the discharge from 
it was to be carried in a pipe to a power canal a short distance away. This 
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flow had to be measured. The circular orifice in the end of the pipe would be 
running partly full most of the time, and therefore became a circular weir. 


Experimental Data 


A few laboratory experiments are available on the flow through circular 
and semi-circular weirs. A summary of all such experiments the author could 
find appears in Appendix A. Doubtless others were made in European countries, 
but are not available to the author. There is also given in Appendix B a table 
of flows through such weirs from 0.25 to 5.0 feet in diameter. The theoretical 
flow was first determined by the rigid formula presented herein which was | 
then modified by an average discharge coefficient disclosed by the above labo- | 
ratory experiments. Appendix C presents a table of the elliptic function in- 
volved in the formula by which the theoretical flow for a weir of any diameter 
can easily be tabulated. 
The following laboratory experiments to determine the flow through circular 
and semi-circular weirs have been found in available technical papers. 


Experiments by H. V. Cone(1) 

Experiments on the flow through two circular and two semi-circular weirs 
are reported. The results are given graphically by 4 curves. No attempt was 
made to develop a flow formula for them. 


Experiments by F. V. Greve(2) 
Extensive experiments on the flow through circular, parabolic, and triangu- 
lar weirs were conducted by the late F, V. Greve at the Engineering Experi- 
ment Station of Purdue University, Lafayette, Indiana. ; 7 
The differential formula he developed for the flow through such weirs, he 
states, could not be solved but the flow through semi-circular weirs was found 
by mechanical integration of his differential equations for several diameters. 
The flows in Greve’s bulletin are all given in pounds per second but they + 
have been converted to cubic feet per second herein. 


Experiments by Eldon R. Dodge(3) 

These were made to determine the flow through circular orifices in steel 
plates fixed to the outlet end of sewer pipes intended for storm drain culverts. 
The orifice plates were attached to the pipe with various eccentricities, that 
is with the centers of some orifices concentric and others not concentric with 
the centers of the pipes. Only the flows for concentric orifices are tabulated 
herein. Flows were determined under heads with orifices flowing partly full 
and also under surcharged heads wherein the discharge changes from weir 
flow to orifice flow. 

Mr. Dodge presents a trigonometrical series formula, aided by curves, by 
which the theoretical flows through his weirs were approximately found. He 
also presents curves of his discharge coefficients. 7 

It is doubtful whether any of the above authors and experimenters were 
aware of the work done in Europe on flow through circular weirs. It is prac- 
tically certain that the first one in this country to make recorded tests of such 
weir flow was H. V. Cone who made his experiments in 1916 at Fort Collins, 
Colorado. Dr. Eng A. Staus’ first publication of which the writer has knowledge 
appeared in 1926, while G. H. Gulliver made flow tests in England on small 
circular weirs as early as 1909. (See References) 
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Experiments by G. H. Gulliver) 


These were made on a small circular weir 2-1/2 inches diameter. The 
differential formula developed is 


h 
q 12.05 [D-(H-hj] dh 


where D = diameter of the circular weir-inches; H = head above the lowest 
point of the weir-inches, q = discharge - gpm; h = distance to elementary area 
above the circular weir invert. The formula was mechanically integrated and 
the results are given by curves. Flow coefficients found were between .625 
and .630. The curves were on too small a scale to be read with accuracy and 
his results are not included in this paper. 


Experiments by H. Th. Thijsse(5) 

Circular weirs have been in use for many years in his laboratory to 
measure the flow through hydraulic models of structures and other devices. 
They were adopted because of their ease of construction but more particularly 


because they did not have to be leveled and the point of zero is so readily de- 
termined. 


The formula given is 


Q=m VgR = kH2 


Where Q, flow in liters per second 

g, the metric gravity constant 

R, radius of Circular Weir - centimeters 

H, head on the Weir invert - centimeters 

m, discharge coefficient 
The value of m was found by calibration in place. Curves were prepared for 
values of k vs H2 from the calibration for each size of weir. These curves 
however have no general validity if the conditions of use depart materially 
from those under which they are obtained. 


These data are given in metric units but they have been converted to English 
units in the tables herein. 


The Elliptic Formula 


The articles by Dr. Eng Staus do not give the transformation to the Elliptic 
form which is rather involved. This however was developed entirely inde- 
pendently by the author with substantial assistance from Dr. W. E. Milne, re- 
tired, formerly head of the mathematics department of Oregon State College. 
This transformation is given in Appendix D. 

Referring to Fig. 1 it is obvious that 


2 (5) - 5 y ) y) (1) 


The head on the elementary area 1 dy is H - y and the elementary theoretical 
flow becomes 


dQ 2 \ D-y) y (H-y) dy 


This is true whether H is greater or less than the radius of the circle. Then 
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Q, 2 \few Ke H - y) dy (2) 


This can be developed into a series and ae term by term obtain- 
ing as was done by E. R. Dodge but not in a trigometric form 


Q, 2 yég 2°) 214, 218, 
Where z = H/D 


This, however, is an infinite series and laborous to deal with. 
The elliptic formula for the theoretical flow through circular weirs is 


Q oD? [ec ke 4 E - (2-«©) (1-k*) x | (3) 


Application of the Elliptic Formula 


In equation (3) the modulus k depends solely on the head and diameter of 
the weir. In brief the formula may be written 
Q = BA 
where 


B = 4/15 V2gD9/2 2.14p9/2 


and Ais the quantity in[ ] 

In evaluating (3) for various diameters a table of k2 vs A was prepared (see 
Appendix C). This table gives A to 4 decimal places for values of D/H to 3 
decimal places. If we multiply the values of A by 2.14 we obtain the theoretical 
flow for a 1-foot diameter weir since 


B = 4/15 V2g D9/2 = 2.14 


The value of the factor B for any other diameter is obtained by 
B = 2.14D9/2 


Discharge Coefficients 


The tabulation of experimental data in Appendix A gives the discharge co- 
efficients obtained by 


c = Q/Q 


These values are plotted in Fig. 2. While there is some evidence of a non- 
linear trend it is not definable in the evidence at hand. No coefficients have 
been discarded. The average of each experimenter’s coefficients for each weir 
diameter is also indicated in Fig. 2. From this it appears that an average co- 
efficient of 59 per cent is justified for the experimental data presented. 

The coefficients from the Thijsse experiments indicate that the flows were 
influenced by velocity of approach. However this would not affect the validity 


| 
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of his experiments since his weirs were all calibrated in place and the flows 
observed include this factor. The Dodge experiments also appear to show 
this same effect since the orifices were in the end of drain pipes. 

The experiments by H. V. Cone are all near the 59 per cent line except 
those for the 3-foot weir for which there is no evident explanation unless his 
zero point was in error. 

The Greve experiments are the most consistent of all and are believed to 
be entirely trustworthy. 

It is hoped that some laboratories in this country will make more experi- 


ments on the flow through circular sharp-crested weirs and thereby obtain 
more discharge coefficients. 


Notation 


In this paper the letter symbols for Hydraulics, approved by the American 
Standards Association January 1942 has been used as far as practicable. 
The foot-seconds system has been used throughout. See also Fig. 1. 


diameter of sharp crested weir opening 


head on weir above the invert measured upstream of the surface curve 
of approach 


theoretical flow through weir 

observed flow 

coefficient of discharge = Q/Q; 

acceleration of gravity 

length of chord of elementary area 

distance above circle invert to an elementary area 

modulus of the elliptic integrals = VH/D 

the complete elliptic integral of the first kind 

the complete elliptic integral of the second kind 

Value of the elliptic function involved such that Qt = BA and Q = BCA 
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APPENDIX A 


Flow through Circular Weirs and discharge coefficients determined by 


laboratory experiments. 


See text for description of each author’s experiments. 


E. R. Dodge 


H/D 


F. W. Greve 


H/D 


-262 
-314 
-320 
-412 
-546 
.610 
-688 
-782 
-880 


Q 


-071 
.101 
-105 
.168 
-282 
-344 
-424 
.925 


H. V. Cone 


.159 
.334 
-437 
-945 
-655 
.758 


E. R. Dodge 


-289 
-437 
-916 
-633 
-750 
-814 
-920 
-933 


-085 
.189 
-256 
-493 
.963 
-680 
-694 
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D Q Q c D H Q 
-211 .084 .398 .018 .012 .666 00 .131 -043 .606 
-100 .474 .025 .016 .640 -060 .594 
108 .512 .029 .019 .655 .160 .061 .580 
-132  .626 .042 .026 .619 .206 .098 .583 
.677 .048 .031 .646 -212 -104 .588 
150 6.71106 «6.052 «836.081 -164 .580 
-152  .720 .053 .032 .604 .305 -200 .580 
-169 .802 .063 .036 .572 .344 244 .575 
.810 .064 .039 .619 .391 -304 .580 
-177 .068 .043 .633 -440 -633  .368 .580 
-189 .896 .075 .046 .614 .585 
-193 .915 .078 .050 .640 
-202  .958 .082 .051 .623 
-211 1.000 .088 .055 -625 
-624 
.50 .20 -40 -10 -628 
.50 -540 
F. W. Greve .30 -60 
.35 -70 
.525 .046 .027 .587 -40 .80 .30 -550 
-133° .533 6584 -45 -90 .38 .580 
.548 .050 .029 .580 =1.00 -44 -980 
-146 .576 .055 .032 .582 .585 
-240 .960 .127 .074 .583 
-244 .976 .129 .076 .990 
.585 .50 .145 -054 .635 
.219 -113 .598 
.259 -154 
E. R. Dodge -317 -218 .592 
.376 -297 =.603 
-336 .128 .380 .054 .034 .628 -407 -340 .605 
.069 .040 .580 -461 -421 .619 
-161 .478 .085 .050 .588 -467 -421 
-168 .500 .089 .055 .618 .607 
.652 .144 .084 .583 
.682 .155 .097 .625 
-235 .700 .162 .096 .593 
.780 .195 .114 .584 
-294 .875 .233 .136 .585 
-299 .888 .239 .139 .583 
-302  .900 .244 .151 .619 
-321 .955 .265 .156 
-602 
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-666 


F. W. Greve 
H/D Qt 
.234 =.117 
-244 =.127 
-298 .187 
.326 
§=.296 
-451 .408 
.003 .091 
-656 .799 
-748 .998 
.832 1.186 
-924 1.392 
1.000 1.550 
294 .185 
-411 .461 
-416 .471 
-475 .603 
-484 .624 
.269 .269 
.332 .402 
-348 .440 
-418 .620 
.869 
-090 1.162 
1.183 
-685 1.508 
.785 1.894 
-794 1.929 
-803 1.965 
.902 2.354 
.906 2.376 
1.000 2.750 


th, 
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So 


F. W. Greve 


H/D 


Q 


.251 
.379 
-754 
1.029 
1.461 
1.491 
2.079 
2.659 
3.530 
4.280 


1.463 
2.128 
2.827 
3.701 
4.639 
6.300 


-247 
-953 
-870 
1.300 
1.684 
2.168 
2.700 
3.788 
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| 
D H Q ike D H | a Q c 
.068 .581 1.00 -205 .205 .146 .580 
.163 .074 .582 208. -220 .580 
-198 .108 .579 -314 .576 
217 .364 .364 435 .577 
aan -430 .430 -994 .578 
.300 286 6.509 .920 .520 836 .572 
.368 .336 220 .570 
.437 -454 .569 .634 1.190 .573 
.498 .068 .569 1.520 .572 
.554 75 .570 .873  .873 2.000 .566 
.615 .798 61.000 2.520 .588 
.665 .922 .594 .576 
.578 
i H. V. Cone 
|_| .191 107 6.579 
.223 -145 .580 1.00 .628 
-268 .570 -899 .52 
.356 245 .572 1.360 .80 .588 
.363 -356 .570 .6 1.887 1.10 .583 
.584 2.466 1.41 .572 
8 3.077 1.76 
9 3.698 2.15 .981 
.833 .189 112 «86.578 1 1.0 4.280 2.50 
.224 .158 .585 
.276 
.290 204 .578 
348 .398 8.577 F. W. Greve 
-417 495 .570 
.491 .659 .568 .598 
.495 -671 .567 316 .271 .633 .595 
| 864 .573 -390 .334 .944 .585 
.654 076 .569 -492 .422 .594 
.661 093 .567 -603  .517 .610 
.669 114.568 -708 .606 
345 .571 .708 .586 
| 352 .569 .947 .812 .985 
| .833 570.574 1.166 1.000 -602 
574 578 
ae 


D 


1.181 


H 


-236 


J. Th, Thijsse 
H/D 


-200 
.212 
.254 
-296 
-338 
-381 


-466 

-907 

-950 

2, 
-635 3.163 
-676 3.526 
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H/D 


-143 
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| .210 .578 1.50 -214 .345 .207 .600 
-406 .246 .605 -226 .151 .384 .230 .599 
.30 -576 .350 .607 -248 .165 .456 .272 .596 
-40 -995 .613 .616 -309 .206 .700 .417 .596 
-45 1.246 .767 .616 .232 .879 .526 .598 
.50 -944 .620 -388 .259 1.086 .643 .592 
1.155 .638 .283 1.289 .755 .586 
-60 1.365 .645 =.314 1.574 .916 .582 | 
.65 1.590 .649 1.779 1.035 .582 | 
.70 1.860 .664 048 .365 2.091 1.205 .576 | 
2.110 .666 -990 .393 2.401 1.383 .576 | 
80 2.410 .680 642.428 2.815 1.620 .576 
-630 -674 .449 3.077 1.780 .579 | 
-720  .480 3.481 2.003 .575 | 
-920 .614 5.422 3.070 .567 
F, W. Greve 1.050 .700 6.804 3.800 .559 h | 
1.196 .798 8.456 4.780 .565 
1.333 .280 .210 .540 .322 .597 1.406 .938 10.830 6.250 
.242 .709 .417 .588 .982 
415 .311 1.149 .675 .588 
.372 1.611 .935 .580 
-977 .433 2.139 1.244 .580 1.748 .159 .091 .204 .128 .627 
982 .437 2.176 1.260 .580 .118 .345 .214 .620 
 .527 3.069 1.797 .586 204 .145 .520 .316 .608 
.610 3.986 2.296 .575 .184 .825 .482 .584 
-967 .726 5.380 3.055 .568 -388 .222 1.183 .710 .600 
1.124 .844 6.873 3.900 .567 435 .249 1.476 .875 .593 
1.146 .860 7.077 3.995 .565 487 .279 1.837 1.080 .589 
1.332 1.000 8.780 5.070 .978 928 .302 2.141 1.250 .584 ie 
.580 .345 2.756 1.540 .560 
-647  .370 3.143 1.830 .582 
-703  .402 3.670 2.140 .583 
H. V. Cone -750 .429 4.142 2.390 .577 
801 .458 4.677 2.710 .579 
1.50 .20 .298 =.18 .604 -853  .488 5.257 3.050 .580 
.30 .660 .38 .991 
.35 233 511 .575 
-40 -2607 1.152 .69 -598 
-50 333 1.758 1.03 
-60 -400 2.480 1.45 .584 
-70 -467 3.308 1.90 .975 
.80 -933 4.211 2.44 -980 
588 
| 


ASCE 


2.25 


_ 


. Thijsse 


Q 


-406 

-507 

-722 
1.263 
1.935 
2.749 
3.663 
4.707 
9.835 
7.068 


. Greve 


-293 

-483 

-640 
1.014 
1.328 
1.846 
2.195 
2.733 
3.071 
3.544 
4.286 
4.828 
5.971 
7.134 
7.520 


1.128 
1.534 
2.441 
3.451 


. 186 
300 
-392 
-615 
-792 
1.097 
1.305 
1.560 
1.808 
2.070 
2.504 
2.810 
3.485 
4.223 
4.430 


.450 

-690 

-840 
1.50 
2.02 


STEVENS 


-670 


.630 


-635 


2.495 


.409 
-453 
-490 
-926 
-063 
.985 
-620 
.665 
.699 
-723 


J. Th. Thijsse 
H/D Qt 


095 .473 
-106 .590 
.127 = .846 
.169 1.484 
-212 2.294 
-254 3.248 
-296 4.362 
.339 5.646 
.382 7.066 
-424 8.583 
-466 10.23 
.908 11.98 
.990 13.85 
-994 15.87 
-635 17.85 
-695 20.86 


F. W. Greve 


125 .942 
-164 1.607 
-181 1.943 
2.265 
.211 2.608 
-226 2.978 
-235 3.209 
-248 3.560 
-267 4.105 
-280 4,498 
.290 4,814 


H. V. Cone 


-100 .953 
-133 1.686 
-167 2.640 
-200 3.730 
-233 5.006 
-267 6.506 
-300 8.142 
.333 9.932 


.283 
.354 
-507 
.894 
1.380 
1.975 
2.680 
3.470 
4.360 
5.430 
6.580 
7.900 
9.330 
10.90 
12.60 
15.20 


.973 

-966 
1.175 
1.359 
1.550 
1.765 
1.895 
2.118 
2.424 
2.648 
2.823 
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.467 


7 
J. Th 
| 1.770 .223  .126 246 .606 2.36 .223 .600 
| «(141 306 .604 .600 
| 30.169 437.605 .30 .600 
| 40.226 .768  .608 .40 .603 
50 1.185 .613 .603 
60  .339 1.700 .618 .60 .608 
70.395 2.330 .636 .70 
80.452 3.060  .650 .80 615 
| .565 4.920 .695 .633 
| .643 
.659 
| .675 
.685 
.706 
| 2.00 .185 .092 .728 
| 237.118 .622 .636 
| 273.136 .612 
345.172 .607 
| 396 .198 
470.235 
515.257 .608 
567.288 .601 
613 .306 .589 .604 
.661  .330 .600 
730.365 585 594 
873.436 .583 .590 
.960 .480 .594 
.988  .494 .589 .590 
.595 .589 
-566 
.594 
8.280) «.124 .630 
351.156 611 | 
413.183 | 
525 .238 .614 
 .279 3.00 48 .504 
.592 .40 81.481 
.50 1.22 .463 
.60 1.72  .462 
.70 2.30 .459 
.80 2.95 .454 
.90 3.67 .452 
1.00 4.44 
| 
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APPENDIX B 


Flow Through Circular Weirs as Calculated from the Following Formula 


Q = 1.262 D?/2 2(1 - k2 - k4) E - (2 - k2)(1 - k2) K 


where 

Q is the flow in ft3/sec 

D, the diameter of the Circular Weir - feet 

H, the head on the Weir invert - feet 

K2, H/D the head in terms of diameter 

E and K taken from tables of complete elliptic integrals 

b=4 C Vag D5/2 where C = .59 

15 
based on the experiments in Appendix A 
Values of b 
/ 
D b D 2 b 
25 .0312 .039 2.00 5.656 7.14 
50 .1768 223 2.50 9.882 12.47 
15 -4871 615 3.00 15.588 19.67 
1.00 1.000 1.262 4.00 32.000 40.38 


| 
1.50 2.756 3.48 5.00 55.900 70.55 


ASCE 


D=.60 D=.75 D=1.00 


STEVENS 


Q 


“oO 
308 
anon 


© 
wo 


OS 


-950 


o- 


690 


D=2.00 D= 


t 
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2.50 D = 3.00 


| D = .25 
| 
H Q H Q Q Q | a Q Q Q . 
| .001 .006 .008 .009 .054 .063 
J .04 .003 10 .025 .031 .036 .204 
.06 .006 15 .054 .069 .081 .514 .563 
.08 .011 094 .141 .908 .995 
10.017 25.142 .184 1.227 1.395 1.558 
.12 .023 .30 .197 .259 .308 1.744 1.981 2.201 
| 14.031 .35 .257 .346 .413 2.336 2.667 2.954 
.039 .321  .441 .530 3.000 3.445 3.840 
18 45 .385 .544 .661 3.739 4.303 4.305 
.20 50 .446 .655 .802 4.535 5.241 5,862 
.22 .066 55 .769 5.390 6.259 7.038 
.24 .60 .888 1.113 6.294 7.352 8.267 
.25 65 1.008 1,281 7.245 8.511 9.583 
.70 1.124 1.454 8.222 9.725 11.03 
15 1.235 1.633 9.237 11.00 12.50 
.80 1.815 10.26 12.32 14.04 + 
.85 2.000 10.98 13.68 15.69 
.90 2.181 12.33 15.06 17.35 
.95 2.357 13.33 16.50 19.06 
1.00 2.524 _ 14.28 17.93 20.87 
1.05 007 19.38 22.66 
1.10 330 20.82 24.49 
1.15 671 22.25 26.42 
1.20 001 23.63 28.29 
1.25 335 24.94 30.17 
1.30 676 32.10 
1.35 6.010 33.99 
| 1.40 6.333 35.82 
nies 1.45 6.657 37.65 
1.50 6.956 39.34 
| 
| 
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D=4.00 D=5.0 D=4.00 D = 5.00 


Q Q 


.081 .092 
-295 .325 
-646 .720 
1.155 1.284 
1.809 2.018 
2.592 2.907 
3.497 3.951 
4.518 9.136 
5.669 6.448 
6.941 7.895 
8.342 9.482 
9.865 11.21 
11.49 13.08 
13.22 15.09 
15.04 17.24 
16.97 19.49 
19.00 21.86 
21.25 24.35 
23.36 26.94 
25.66 29.65 
28.04 32.47 
30.49 35.41 
33.02 38.46 
35.61 41.60 
38.26 44.83 
40.99 48.15 
43.73 51.54 
46.52 95.22 
49.38 58.58 
52.25 62.21 


CO OO 


1. 
2. 
2. 
2. 
2. 
2. 
2. 
2. 
2 

2 

2. 
3 


Note: Although Q is given for the complete circle of the weir, the flow begins 
to increase at about 70% of the diameter beyond that given by the 
formula, thus, partaking of a transition between weir flow and orifice 
flow as the weir becomes surcharged. : 


H Q Q 
55.16 65.91 
58.07 69.68 
60.97 73.51 
63.92 77.39 
66.83 81.27 
69.78 85.22 
72.60 89.24 
75.83 93.34 
78.14 97.28 
80.76 101.5 
105.5 
109.6 
113.7 
117.8 
121.9 
125.9 
129.8 
133.7 
137.4 
141.1 


STEVENS 


APPENDIX C 


Values of A = 2(1 - k2 - k4) E - (2 - k2)(1 - k2) K 


where k2 = H/D; E and K to be taken from tables 
of complete elliptic integrals. 


Note: for k2 values the 3d decimal is at the top. 


| 
ASCE a 1455-17 
| 
| 
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-0004 
-0013 
-0027 
-0046 
-0071 


-0102 
-0139 
.0182 
-0231 
-0286 


-0346 
-0412 
.0483 
-0560 
-0642 


-0728 
-0819 
-0914 
.1014 
-1119 


.1229 
.1344 
1464 
-1589 
.1719 


-1854 
-1994 
-2139 
-2289 
-2443 


-2601 
-2763 
-2929 
-3099 
.3273 


-3451 
3633 
-3819 
-4009 
-4203 


-4401 
-4603 
-4809 
-5019 
-5233 


-9451 
-9672 
.5896 
-6123 
.6354 


ouwo ule 


| 
HY 6 
k2 1 2 3 4 5 6 7 8 9 Diff 
for k2 
- ,001 
0 oO 0000 .0001 .0001 .0002 .0002 .0002 .0003 .0003 .0004 0 
.01 .0005 .0006 .0007 .0008 .0009 .0010 .0011 .0011 .0012 0 
.02 0015 .0016 .0018 .0019 .0020 .0022 .0023 .0025 .0026 1 
.03 .0028 .0030 .0032 .0034 .0035  .0037 .0039 .0041 .0043 1 
.04 0048 .0050 .0053 .0055 .0058  .0060 .0063 .0065 .0068 2 
.05 .0077 .0080 .0083 .0087  .0090 .0093 .0096 .0099 3 
.06 0106 .0109 .0113 .0117 .0120 .0124 .0128 .0132 .0135 3 
.07 0144 .0148 .0152 .0156 .0160 .0165 .0169 .0173 .0178 4 
.08 0187 .0192 .0197 .0202 .0206 0211 .0216 .0221 .0226 4. 
.09 10236 .0242 .0248 .0253 .0258  .0264 .0270 .0275 .0281 ‘5. 
.0298 .0304 .0310 .0316  .0322 .0328 .0334 .0340 6. 
0353 .0359 .0366 .0372 .0379  .0386 .0392 .0399 .0405 6. 
0419 .0462 .0433 .0440 .0448  .0455 .0462 .0469 .0476 
0491 .0498 .0506 .0514 .0522 .0529 .0537 .0545 .0552 7. 
14 0568 .0576 .0585 .0593 .0601  .0609 .0617 .0626 .0634_ 8. 
115 0651 .0659 .0668 .0676 .0685  .0694 .0702 .0711 .0719 8. 
.16 .0746 .0755 .0764 .0773  .0783 .0792 .0801 .0810 9. 
17 .0838 .0847 .0857 .0866  .0876 .0886 .0895 .0904_ 9. 
.18 0924 .0934 .0944 .0954 .0964 0974 .0984 .0994 .1004 10. 
.19 1024 .1035 .1046 .1056 .1066 1077 .1088 .1098 .1108 10 
.20 1130 .1141 .1152 11163 .1174 1185 .1207 .1218 11. 
.21 1240 .1252 .1264 .1275 .1286 .1310 .1321 .1332 11.5 
.22 1356 .1368 .1380 .1392 .1404 1416 .1428 .1440 .1452 12.0 
.23 .1489 .1514 .1526 .1564 .1576 12.5 
.24 .1615 .1628 .1641 .1654 .1680 .1693 .1706 13.0 
.1746 .1773 .1786 .1814 .1827 .1840 13.5 
.26 .1882 .1869 .1910 .1924 .1938 .1952 .1966 .1980 14.0 
.27 .2008 .2023 .2038 .2052 .2066  .2081 .2095 .2110 .2124 14.5 
.28 .2169 .2184 .2199 .2214 .2244 .2259 .2274 15.0 
.29 -2304 .2320 .2335 .2351 .2366 2381 .2397 .2412 .2428 15.4 
.30 2459 .2475 .2490 .2506 .2522 .2538 .2554 .2569 .2585 15.8 
.2633 .2650 .2666 .2682 .2698 .2714 .2731 .2747 16.2 
.32 .2780 .2796 .2813 .2829 .2846  .2863 .2879 .2896 .2912 16.6 
.33 2946 .2963 .2980 .2997 .3014 .3031 .3048 .3065 .3082 17.0 
.34 3116 .3134 .3151 .3169 .3186  .3203 .3221 .3238 .3256 17.4 
.3291  .3309  .3326 .3344 .3362 .3398 .3415 .3432 17.8 
.36 3469 .3487 .3506 .3524 .3542 3560 .3578 .3597 .3615 18.2 
.37 3652 .3670 .3689 .3707 .3725 .3745 .3763 .3782 .3800 18.6 
.38 3838 .3857  .3876 .3895 .3914 .3933 .3952 .3971 .3990 19.0 
.39 4028 .4048 .4067 .4087 .4106 4125 .4145 .4164 .4184 19.4 
.40 4223 .4243 .4262 .4282 .4302 .4342 .4361 .4381 19.8 
41 4421 .4441 .4462 .4482 .4502 .4542 .4563 .4583 20.2 
.42 4624 .4644 .4665 .4685 .4706 .4747 .4768 .4788 20.6 
.43 4830 .4851 .4872 .4893 .4914 4935 .4956 .4977 .4998 21.0 
.44 5040 .5062 .5083 .5105 .5126 5147 .5169 .5190 .5212 21.4 
45 5255 .5277 .5298 .5320 .5342 .5364 .5386 .5407 .5429 21.8 
.46 5473 .5495 .5517  .5539 .5561 5584 .5606 .5628 .5650 22.1 
5694 .5717 .5739 .5762 .5784 .5806 .5829 .5851 .5874 22.4 
.48 5919 .5964 .5987 .6010 .6032 .6055 .6078 .6100 22.7 
.49 -6146 .6169 .6192 .6215 .6238  .6262 .6285 .6308 .6331 23.1 
.50 6377 .6401 .6424 .6448 .6471 6494 .6518 .6541 .6565 23.4 


ASCE 
k2 0 

.51 .6588 
.52 .6825 
.53 .7064 
.7306 
.7551 
.7799 
.57 .8050 
58 .8304 
.59 .8560 
.60 .8818 
.9079 
.62 .9342 
.63 .9608 
.64 
.65 1.0147 
.66 1.0420 
.67 1.0694 
.68 1.0969 
.69 1.1246 
1.1524 
1.1804 
.72 1.2085 
.73 1.2368 
.74 1.2653 
.75 1.2939 
.76 1.3226 
17 1.3514 
.78 1.3802 
.79 1.4091 
.80 1.4380 
.81 1.4670 
.82 1.4960 
83 1.5250 
.84 1.5540 
.85 1.5830 
.86 1.6120 
.87 1.6410 
.88 1.6699 
.89 1.6988 
.90 1.7276 
91 1.7561 
.92 1.7844 
.93 1.8125 
.94 1.8403 
.95 1.8678 
.96 1.8950 
.97 1.9219 
.98 1.9484 
.99 1.9744 
1.00 2.0000 


_ 
. 


-6659 
.6897 
-7137 
-7380 
-7625 


-7874 
-6126 
-8381 
.8637 
.8896 


.9158 
.9422 
-9688 
-9957 
1.0229 


1.0502 
1.0777 
1.1052 
1.1329 
1.1608 


1.1888 
1.2170 
1.2453 
1.2739 
1.3025 


1.3312 
1.3600 
1.3889 
1.4178 
1.4467 
1.4757 
1.5047 
1.5337 
1.5627 
1.5917 
1.6207 
1.6497 
1.6786 
1.7074 
1.7362 


1.7650 
1.7928 
1.8208 
1.8485 
1.8760 


1.9031 
1.9298 
1.9562 
1.9821 


STEVENS 
4 5 

-6683 .6706 
-6921 .6944 
7161 .7185 
7404 .7428 
7650 .7675 
7899 .7924 
8152 .8177 
8406 .8432 
-8663 .8689 
8922 .8948 
9184 .9210 
.9448 .9475 
.9715 .9742 
.9984 1.0012 
1.0256 1.0284 
1.0530 1.0557 
1.0804 1.0832 
1.1080 1.1108 
1.1357 1.1385 
1.1636 1.1664 
1.1916 1.1944 
1.2198 1.2226 
1.2482 1.2510 
1.2767 1.2796 
1.3054 1.3082 
1.3341 1.3370 
1.3629 1.3658 
1.3918 1.3946 
1.4207 1.4236 
1.4496 1.4525 
1.4786 1.4815 
1.5076 1.5105 
1.5366 1.5395 
1.5656 1.5685 
1.5946 1.5975 
1.6236 1.6265 
1.6526 1.6554 
1.6815 1.6843 
1.7103 1.7132 
1.7390 1.7418 
1.7674 1.7702 
1.7956 1.7984 
1.8236 1.8264 
1.8513 1.8540 
1.8787 1.8814 
1.9058 1.9084 
1.9325 1.9352 
1.9588 1.9614 
1.9846 1.9872 


. 
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Diff 
2 3 6 7 8 9 for k2 
= .001 
.6612 .6635 .6730 .6754 .6778 .6801 23.7 
.6849 .6873 6968 .6992 .7016 .7040 23.9 
.7088 .7112 .7209 .7233 .7258 .7282 24.2 
.7330 .7355 .7453 .7478 .7502 24.5 
.7576 .7700 .7725 .7749 .7774 24.8 
.7824 .7849 .7950 .7975 .8000 .8025 25.1 
.8075  .8101 .8202 .8228 .8253 .8279 25.4 
.8330  .8355 8458 .8483 .8509 .8534 25.6 
.8586 .8612 .8715 .8741 .8766 .8792 25.8 
.8844 .8975 .9001 .9027 .9053 26.1 
.9105 .9132 .9237 .9289 .9316 26.3 
.9369 .9502 .9528 .9555 .9581 26.6 
| .9635 .9662 9769 .9796 .9822 .9849 26.8 
9903 .9930 1.0066 1.0093 1.0120 27.1 
| Hoi74 1.0202 0311 1.0338 1.0365 1.0393 27.3 
0447 1.0475 584 1.0612 1.0639 1.0667 27.4 
0721 1.0749 D859 1.0886 1.0914 1.0941 27.5 
0997 1.1024 1135 1.1163 1.1191 1.1218 27.7 
| 1274 1.1302 1413 1.1441 1.1468 1.1496 27.8 
1552 1.1580 1692 1.1720 1.1748 1.1776 28.0 
: 832 1.1860 1973 1.2001 1.2029 1.2057 28.1 
1.2142 2255 1.2283 1.2311 1.2340 28.3 
2396 1.2425 2539 1.2568 1.2596 1.2624 28.5 
2682 1.2710 2825 1.2853 1.2882 1.2910 28.6 
| 2968 1.2996 B111 1.3140 1.3169 1.3197 28.7 
| B255 1.3284 1.3399 1.3428 1.3456 1.3485 28.8 . 
| 8543 1.3572 1.3716 1.3744 1.3773 28.8 
8831 1.3860 1.4004 1.4033 1.4062 28.9 
4120 1.4149 264 1.4293 1.4322 1.4351 28.9 
4409 1.4438 1.4583 1.4612 1.4641 29.0 
1.4699 1.4728 1.4873 1.4902 1.4931 29.0 
| 1.4989 1.5018 5134 1.5163 1.5192 1.5221 29.0 7 
1.5279 1.5308 D424 1.5453 1.5482 1.5511 29.0 
1.5569 1.5598 $714 1.5743 1.5772 1.5801 29.0 
1.5858 1.5880 5004 1.6033 1.6062 1.6091 29.0 
1.6149 1.6178 5294 1.6323 1.6352 1.6381 29.0 
1.6439 1.6468 5583 1.6612 1.6641 1.6670 28.9 
1.6728 1.6757 5872 1.6901 1.6930 1.6959 28.9 
1.7017 1.7046 161 1.7190 1.7218 1.7247 28.8 
1.7304 1.7333 447 1.7475 1.7504 1.7532 28.5 
1.7589 1.7618 731 1.7759 1.7787 1.7816 28.3 
. 1.7872 1.7900 8013 1.8041 1.8069 1.8097 28.1 
1.8153 1.8181 8292 1.8320 1.8347 1.8375 27.8 
1.8430 1.8458 8568 1.8596 1.8623 1.8650 27.5 
1.8705 1.8732 1.8841 1.8868 1.8896 1.8923 27.2 
1.8977 1.9004 9111 1.9138 1.9165 1.9192 26.9 
1.9246 1.9272 9378 1.9404 1.9431 1.9458 26.5 
1.9510 1.9536 9640 1.9667 1.9692 1.9718 26.0 
. 1.9770 1.9795 9900 1.9923 1.9949 1.9974 25.6 
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APPENDIX D 


Transformation of the basic hydraulic equation (2) to the elliptic equation 
(3) in the text; (20) in this appendix. 


ASCE STEVENS 1455-21 


Transformation to Elliptical Integrals 


Equation (2) may be transformed into a form of complete elliptic integral 
and a solution found to almost any degree of precision desired. Tables are 
available of these elliptic functions to 8 or more decimal places. 

Following are a few fundamental relationships that are essential to this 
transformation. The elliptic integral of the first kind is 


ite <u: f {2 sin 


ke sin @ 
= 

¢ is the amplitude and k the modulus 751.1 
%=amu 751.2 
sin = snu=x 751.3 
cos $ = cnu= V1 - x2 751.4 
dnu = V1 - k2x? 751.5 
am -u=-amu 753.1 am 0=0 754.1 
sn - u = -snu 753.2 sn 0=0 754.2 
cn -u=cnu 753.3 ca = 1 754.3 
dn -u=dnu 753.4 dn 0=1 754.4 
sn2u+cn2u=1 755.1 
dn2 u +k? sn2u = 1 755.2 
dn2u - k2cn2u = 1 - k? 751.7 and 755.3 
d_snx=cnxdnx 768.1 


dx 
The elliptic integral of the second kind is 


dx [x sin? | 


yi-x Dwight 771 


If the integration is taken between the limits 0 and 7/2 we obtain complete 
elliptic integrals thus 


| 1-k 2sineg 


Dwight 773.1 


1455-22 HY 6 
and 


ag 
Let 
y = Hsn2 x = Dk? sn? x 
H = Dk2 


The first derivative of (3) is 


dy = 2Hsnxd snx 


dx 
and from (4) and Dwight 768.1 
dy = 2D k* snx cnx dn x dx 
from(3) and 755.1 
H - y = H(1 - sn2x) = D k@cn2x 
from (4) and 755.2 


D - y = D(1 - k@sn2x) = D dn2x 
From (3), (4), (6), (7) and (8) 


\(D-y)y (H-y) * ay 
|Ddn*x Dk? 


} 
[4D5x8sn ‘x cn x dn? x dx 


2 2 42 
sn“x D k* | (2Dk*sn x cn x dn x) dx 


December, 1957 


Dwight 774.1 


(3A) 


(4) 


(5) 


(8) 


(9) 
(3) may be written y =H sn2 (sin-! ¢). Hence wheny = 0, sin"! ¢ =x =0 


also when y = H, from Dwight 751.3 since y/H = 1 
x 


u = K the complete elliptic integral of the first kind. 
Adding the limits of integration and substituting (9), (2) becomes 


K 
> 2 
Q, => 4V2g p°/2 k4 


But by Dwight 755.1 and 755.2 


en2x =j- sn2x and dn2x -l- k2sn2x 
substituting in (10) 


4 Vag xt 
where 


(10) 


(11) 


__ 
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is | sn*x dx - (1 - dx + k 


u 


rk 


2 


sn dx 


% (12) 
Reverting to Pierce) we have 


(m+1) sn™ x dx = (m+ 2) (1+ k*) “sn ™+2x ax 


+4 
-(m + 3)k ! sn™* *x dx + sn™+ 1, cnxdnx, Pierce 567 


‘which for m = 2, and adding the limits of integration becomes 
7K 
3 dx = 4(1 +k?) sn4dx - 5k2 sn®x dx +sn3x cn x dn x 
“9 
since by Dwight 751.3 and 754.2 sn 0 - 0 and cn K =x = 0 theintegrated terms 
(the last one) of Pierce 567 disappears and there remains 


K 
k© sn®°x dx = 4(l+k ) j sin x dx -3 a x dx (13) 


o 


also for m = 0 in Pierce 567 
-K 
2 
dx = 2(1 + k*) sn x dx - 3k~ | sn*x dy 
“o “o “0 (14) 
Since the left hand term = K we obtain by Pierce 564 


| on”x (K- E) 
k 


(15) 
since 
E(o, k) = 0 and E(am k, H) = E 
substitute (15) in (14) 
K 
3k? | sn*x dx = 2(1+k2) (K-E)-K (16) 
or 
rk 
sn4dx = 2(1+k*) (K-E)-K 47) 
° 3k4 3k 
substitute (15) in (17) and (13) 
pK 
5k? | snox dx = 4(1 +k?) 2(1+k2) (K-E)-K | 
0 
= 8(1+k2)* (K - E) - 4(1+k2) K-3 (K-E) 
34 3 ke 
whence 
k? | sn6x dx =| 8 (142)? -3 | [K-E| -4 1+«?k 
0 (15 ké 15 ke (18) 
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substitute (15), (17) and (8) in (12) and obtain 


2 
2(1+k*)* (K-E) + 1+k2 K 
3 3k? 
+( 8 -3 )(K-E) - 4 1+k* 


15 5k 15 ké 
collecting E - terms 


E  (- + 10+ - 8 - 16k? 8k? + 
15k? 


2(1-k*+k4)E 


collecting K - terms 


K (15k* 10 - 20k? - 10k + 5k2 + 5k? 4 84 16k2 


15 k4 
2 - 3k“ + 
15 k4 
combining E - and K - terms 
3 [a1 - - (2 - - K] (19) 
15k4 
substitute (19) in (11) and obtain finally 
a = J2zg - k2 + - (2 - k*) - k*) K] 


(20) 


a rigid formula for the theoretical flow through a sharp edged circular weir 
with complete contractions and no appreciable velocity of approach. 
In (20) k2 = H/D the head on the invert of the circle in terms of the diameter. 
E and K may be taken from tables (see references). 
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PROPORTIONAL WEIRS FOR SEDIMENTATION TANKS* 


Closure by J. C. Stevens 


a. G. Stevens, ! M. ASCE.—The author is grateful indeed to Prof. Guido 
Di Ricco of the University of Rome, for his contribution to the subject of 
proportional weirs. 

He is quite correct in assuming that the author was entirely unaware of 
the very extensive work done by him in this field. This must be charged to 
the handicap of the language barrier that separates us. To illustrate: when 
Prof. Ricco’s discussion appeared any attempt to obtain copies of his papers 
only met with frustration. From his footnotes, they could not be found in the 
Engineering Societies Library in New York nor any other library where in- 
quiries were made. They were finally obtained through ASCE for whom Prof. 
Ricco had photostat copies made of some of the papers cited, which were 
passed on to the writer. Two translations of the principal paper were ob- 
tained from two local Italian sources, but both had to be retranslated by the 
author into the King’s English no doubt because neither translator know hy- 
draulic terms — more frustrations. 

The credit for inventing the proportional weir must undoubtedly be given 
to the late Dr. Oscar Van Pelt Stout, M. ASCE, while he was Professor of 
Civil Engineering at the University of Nebraska. He designed, tested and put 
to practical use such a weir in 1896 and presented a paper on it at the first 
meeting of the Nebraska Engineering Society held in Lincoln in 1897. 

These facts were submitted to the author by Prof. Steponas Kolupaila of 
Notre Dame University, who also submitted a bibliography on the subject 
which, with his permission, is being appended hereto. It is always best to 
give full credit where credit is due. 

The writer was a student under Prof. Stout. He faintly remembers being 
told about the Stout proportional weir while in college, but did not know or 
had forgotten that its mathematical details had been made public some 8 
years prior. During the intervening 50 years it was easy to forget until Prof. 
Kolupaila’s most excellent system of references brought the facts cogently to 
mind, for which grateful acknowledgment is hereby made. 

There is a typographical error in the footnote to Prof. Ricco’s discussion 
as printed in No. HY 1, Feb. 1957 p 1177-41. The date of his paper in 
L’Energia Elettrica should be October 1936 instead of 1956. 

Beginning with the author’s basic equation(7) Prof. Ricco presents his 
generalized formula(2) in which the lateral width of the weir (x) for any depth 
(y) is evolved in his equation(4) in terms of the gamma function. 


a. Proc. Paper 1015, June, 1956, by J. C. Stevens. 
1. Cons. Hydro. Engr., Stevens & Thompson Engrs., Portland, Ore. 
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Had the author known of this he might have saved himself much of the 
laborious series developments with which his paper is adorned. 

The use of a series of chords instead of a purely curvilinear profile for 
the weir is also an innovation that may well be adopted in this country. 
Compare Figs 4 and 5 with the author’s Fig 2. It appears that by this expe- 
dient the proportional weir can be designed to produce almost any desired 


pattern of outflow from a tank or reservoir and might, with profit, be applied 
to some cases of river regulation. 


December, 1957 


It is hoped that this paper and the discussion will stimulate a more exten- 
sive use of the proportional weir than exists at present in this country. It has 
many advantages and could well be used more extensively. 
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THE PROBLEM OF RESERVOIR CAPACITY FOR LONG-TERM STORAGE® 


Closure by Aly Fathy and Aly Salem Shukry 


ALY FATHY! and ALY SALEM SHUKRY“.—The authors find that the dis- 
cussions of their paper have brought forth no consequential information on the 
subject under study with the exception of Mr. Hanna’s exposition on the ques- 
tion of sub-annual storage. That exposition fills a gap in the inquiry that was 
passed over rather lightly. 

In Mr. Mostafa’s discussion an attempt was made to demonstrate the in- 
adequacy of a normal capacity of 45 milld. m3, as arrived at through applica- 
tion of the authors’ method to the projected long-term reservoir on the Nile 
at Aswan. In that attempt, however, sub-annual effects were introduced. In 
the light of Mr. Hanna’s observations it will be seen that when these effects 
are taken into consideration the extra capacity provided for flood relief should 
be included in the working capacity. According to Mr. Mostafa’s own calcula- 
tions the total capacity of 75 milliards (45 + 30) would be more than sufficient. 

The discussion by Messrs. Simaika and Boulos is devoted mainly to 
demonstrating the unreliability of conditions derived from past observations 2g 
as a guide in forecasting the future behaviour of a natural event. In a field of 
enquiry where probability considerations play the leading part such an attempt 
is really superfluous. As far as hydrological phenomena are concerned, the 
principal criteria to which reference may be made are: 


1. The Mean. 

2. The Standard Deviation. 

3. The Maximum Deviation of a single observation from the Mean. 
4 


. The Maximum Deviations of the arithmetic mean for groups of observa- 
tions of varying number. 


Of these criteria, the Standard Deviation is the most erratic and the least 
amenable to provisional adjustment. Thus a method of treatment which de- 
pends primarily on that parameter would be the least reliable. 

In Dr. Hurst’s discussion a slip has been made in the attempt to recon- 
struct the expression of the ideal range (Eq. b). As the length of the critical 
period is supposed = N/2, the derived expression should be multiplied by V2 . 
That would raise the numerical factor to the same value as given by the 
authors in Eq. 15. 

Dr. Hurst considers it faulty to adopt the mean value N/2 as a standard for 
the length of the critical period. It would certainly be more legitimate to do 
so with regard to a parameter that fluctuates about a definite mean than to 


a. Proc. Paper 1082, October, 1956, by Aly Fathy and Aly Salem Shukry. 
1. Part-time Prof., Faculty of Eng., Alexandria Univ., Alexandria, Egypt. 
2. Asst. Prof., Faculty of Eng., Alexandria Univ., Alexandria, Egypt. 
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one that varies progressively such as the range R. According to Dr. Hurst’s 
“Century Storage” theory a standard value is assigned to R corresponding to 
a working period of 100 years. It is here that disregard of the variability of 
R may be questioned. 

Incidentally, it may be noted that by assigning a standard value of 1/2 to 
the relative length of the critical period the resultant expression of R is 
brought to the same form (though not the same value) as arrived at by Dr. 
Hurst himself through a very elaborate mathematical analysis. Mr. Hanna’s 
observations on that point will be found most illuminating. 

The authors wish to stress the fact, however, that in the practical applica- 
tion of their method of treatment no reference at all is made to the ideal 
range which is an indeterminate quantity. The capacity S needed to guarantee 
a draft less than the mean (as would always be the case) is derived straight- 
forward from the deviation-curve. The tying up of S with R, as proposed by 
Dr. Hurst, greatly complicates the problem. Both of his expressions of S 
(Eq. 8 in the original paper and Eq. 27 in Mr. Hanna’s discussion) will be 
found to lead to absurd results when the difference between the guaranteed 
draft and the mean is large. According to (8), S becomes zero when (M - B) = 
1.04 o . That quantity is usually much in excess of the maximum deviation of 
one observation or of the mean for a small group. For the Nile at Aswan, 

M = 93 and o = 19.8. Thus S becomes zero when B = 72.4. That figure is 
30.4 milliards greater than the observed minimum supply of one year, which 
is 42. Evidently, the storage provided should cover the deficit in one year at 
least. 

The logarithmic relation (27) yields a value of S = 4.9 for a guaranteed 
draft of 60 millds. /ann. Actually, the storage needed to cover the deficit in 
the worst year would be 18 milliards. 

The above examples bear sufficient evidence of the need for a totally dif- 
ferent line of approach. It is the authors’ belief that the line they have sug- 
gested is the only alternative. 


Corrections: Fig. 7 of the original paper did not contain the construction 
referred to in the text. The correct figure is reproduced herewith. 


DISC USSION 
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Fig. 7 - DERIVATION OF THE RANGE 
FROM THE DEVIATION-CURVE 
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TECHNICAL PROBLEMS OF FLOOD INSURANCE* 


Closure by H. Alden Foster 


H. ALDEN FOSTER,* M. ASCE.—The discussion by Mr. Alexander em- 
phasizes some of the difficulties encountered in preparing probability esti- 
mates for floods. From a practical viewpoint, it is impossible to prepare an 
absolutely reliable probability curve from a limited flood record. The record 
itself may be represented by a “standard” flood-frequency curve; but the cor- 
responding curve for any future “sample” of equal length may be appreciably 
different within so-called “tolerance limits” or “confidence limits.” 

It has been pointed out by numerous writers that the probability method is 
not suitable for the determination of the flood magnitude for design of the 
spillway of an important dam. As pointed out by the writer, however, in the 
present paper and on other occasions,(1) probability methods have a proper 
and useful application in certain fields, particularly in estimating average an- 
nual flood benefits for a flood-control project or in the design of cofferdams 
for sites subject to floods, where precision of the estimate is not of primary 
importance. For satisfactory determination of flood-insurance rates, some 
kind of probability estimate is essential. In this case, however, the possible 
range of the probability estimate between the confidence limits may add ap- 
preciably to the risk assumed by the insuring agency. 

Mr. Alexander points out the differences of opinion among mathematicians 
as to the proper method of estimating the confidence limits for a particular 
record which is essentially a statistical problem. A fellow-countryman of 
Mr. Alexander(2) has discussed methods for determining such confidence 
limits for the log-normal and Pearson Type III distributions. 

The writer agrees with Mr. Alexander that the method of estimating floods 
by use of the “maximum possible precipitation” requires clarification, par- 
ticularly as to its relation to probability estimates. It is understood that a 
paper by Mr. Alexander discussing this question, together with the interesting 
subject of the relation of probability estimates to “factors of safety” for dams, 
will soon be published by the Institution of Civil Engineers of Australia. 

Professor Kolupaila’s remarks about the writer’s 1924 paper are greatly 
appreciated. The problem of “sampling errors” and confidence limits was 
not well understood by most engineers and hydrologists 30 years ago. An ap- 
preciation of the importance of these problems would doubtless have pre- 


vented unsuitable application of the theoretical probability methods in hydro- 
logical studies. 


a. Proc. Paper 1165, February, 1957, by H. Alden Foster. 


*Prin. Associate, Parsons, Brinckerhoff, Hall & Macdonald, New York, 
N. Y. 
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A HIGH-HEAD CAVITATION TEST STAND FOR 
HYDRAULIC TURBINES* 


Discussion by Robert T. Knapp 


ROBERT T. KNAPP, ! M. ASCE.—The authors are to be congratulated upon 
the design and construction of a fine cavitation laboratory for the testing of 
hydraulic turbines. Their careful planning should pay dividends for a long 
time. The design of the dynamometer and the method of mounting the test 
machine are particularly noteworthy. There are many other nice things that 
can be said about this new laboratory. However, in a discussion it sometimes 
seems more fitting to raise differences in viewpoint. 

One of the possibilities sacrificed in the unique turbine mounting is that of 
visual observation of cavitation in the wicket gates or in any part of the instal- 
lation except the draft tube and the discharge edges of the runner vane. This 
may be serious in the use of the laboratory as a development tool. Another 
unfortunate limitation, in the writer’s opinion, is that the test machine cannot 
be operated as a pump and that transient conditions cannot be investigated. 
This seems unfortunate in light of the growing interest in pump-turbines. 

The writer is surprised at the weight put on the simulation of field test 
techniques. Such techniques are often dictated by necessity rather than by 
desirability. It should not be difficult, with properly selected instruments, to 
get at least equal accuracy and reliability with one operator, with or without 
an assistant, and possibly require less time to make the same run. In this 
connection, it would be interesting to know whether the laboratory has a 
planned program for the regular recalibration of the venturi meters by the 
same agencies using the same test installations that were employed in the 
original calibration. Such a program would appear necessary if the quantity 
measurements are to be maintained at the same high standard of accuracy as 
that obtainable with the other test measurements. 

The writer does not feel that it is permissible to adjust the air content 
solely on the basis of preventing the accumulation of air in the manometer 
connections. Unfortunately, there is as yet not enough known about the rela- 
tionship between air content and cavitation performance to establish a better 
criterion. It is hoped that the laboratory operators will bear this possible 
relationship in mind so as not to become too enamoured of their current 
technique. 

The samples of test results presented in Figs. 4 and 5 are very interesting. 
The authors are to be complimented on not using a “step-up” factor on their 
measured efficiencies. At the Reynolds numbers at which these tests were 
made a step-up would be questionable. The paper seems to imply that 


a. Proc. Paper 1201, April, 1957, by W. G. Whippen and G. D. Johnson. 
1. Prof., Hydr. Eng., California Inst. of Technology, Pasadena, Calif. 
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continuous operation of the turbine at “critical 0 ” would not produce signifi- 
cant cavitation damage. Further comments on this point would be appreciated. 
It would be interesting to see “sigma break” curves for a series of turbines 

of differing designs. The writer is under the impression that it is not uncom- 
mon to obtain appreciably steeper breaks. In this connection it would seem 
desirable to emphasize the authors’ warning that if emergency operation at 
very low sigma values is anticipated, the permissible limits of operation must 
be based on a complete and reliable set of model test curves covering the 
sigma values likely to be encountered. 


The writer agrees with the authors that it is most desirab!< ‘o test the ef- 
ficiency and cavitation performance of hydraulic machines « ».0'o' ype heads 
and velocities. He made the same choice in 1933-34 durin. ‘he desicn of the 


Hydraulic Machinery Laboratory of the California Institute v nology and 
has yet to regret it. However, it must be remembered that this approach 
ignores several facets of the cavitation process. For example, the vertic:! 
dimensions of large machines are great enough to affect the extent of the 
areas subjected to cavitation. To simulate these effects model tests at proto- 
type Froude numbers are necessary. But experiments show that such tests 
also fail to agree with the cavitation performance of the prototype, primarily 
because they fail to alter the properties of the test liquid to correspond to the 
changed scale. It is feared that this dilemma will continue to confront not 
only the authors, but also all other workers in cavitation laboratories for 
some time to come. 

In their concluding paragraphs the authors propose that accurate model 
tests, carried out in adequate laboratories such as the one they describe, can 
be substituted for the more expensive and the consuming field tests. In 
general, the writer agrees with this proposal. He would like to emphasize, 
however, that such model tests to be acceptable must include all of the hy- 
draulically significant features of the field installation. Furthermore, the 
model turbine must be a true model of the machine as installed in the field. 
In many cases this would require the construction of a new model turbine, 
since the model used in the design study may not be geometrically similar 
to the turbine as installed. Two common deviations from similarity are: 

(1) design changes introduced between the model test and the prototype con- 
struction, and (2) errors introduced in the manufacture and installation of the 
full-sized machine that cause it to differ from the final design drawings. In 
the opinion of the writer, if model acceptance tests are to replace field tests, 
the laboratory responsible for making these tests should also have the 
responsibility of ascertaining and certifying to the geometric similarity of 


the model and field installations, as well as to the laboratory accuracy of the 
test measurements. 


THE EFFICACY OF FLOOR SILLS UNDER DROWNED 
HYDRAULIC JUMPS# 


Discussion by Alvin J. Peterka, 
Alfred C. Ingersoll, and S. Leliavsky 


ALVIN J. PETERKA,! M. ASCE.—The conclusions drawn by Mr. Shukry 
are contrary to general practice procedures which have been developed over 
a period of years from hydraulic laboratory tests, prototype tests, and other 
prototype experiences. From a limited number of velocity measurements on 
a small hydraulic model, Mr. Shukry has concluded that a sill placed near the 
middle of the apron is in the optimum position. His conclusion is based on an 
unproved assumption that “The efficiency of any sill against bed erosion is 
indicated by the rate of adjustment of the flow to the normal distribution in 
the downstream channel.” In addition, the tests were conducted on a fixed- 
bed model where erosion could not be measured, and the conclusions are con- 
trary to test results obtained from a movable bed model of the same struc- 
ture, Figure 2. 

It is not clearly stated in the paper whether the difficulty in reconciling 
erosion test results was between several models of the same structure or 
between the models and the prototypes. The “erratic results” of the proto- 
type scour tests are not described. It is not clear, therefore, why Mr. 

Shukry abandoned the conventional movable bed tests and used velocity 
traverses to indicate scour tendencies. It is also not clear why Mr. Shukry 
has mathematically defined a scour hole since he used a fixed-bed model in 
his tests. However, certain impressions are obtained in reading this paper, 
and it is on these impressions the writer must base his discussion. 

The fact that agreement between model and prototype erosion tests was not 
found may be due to any one of several factors. It should be remembered that 
models are mechanical devices and that the interpretation of the model test 
results is often as important as the modei itself. A 3-bay model will repre- 
sent a 46-bay prototype spillway only if all factors concerning operating pro- 
cedures are considered. By proper selection of model scale and bed material 
size, it is possible to predict qualitatively where the most severe erosion will 
occur in a prototype structure. It is also possible to construct several or 
more models of the same structure and obtain identical erosion patterns. 
Contrary to the statement in the paper, it is possible, therefore, to obtain 
qualitative similarity between model and prototype and between different 
models of the same structure. On the other hand, it is not generally believed 
that an ordinary movable bed model can be used to determine quantitative 


a. Proc. Paper 1260, June, 1957, by Ahmed Shukry. 


1. Hydr. Engr., Bureau of Reclamation, U. S. Dept. of the Interior, Denver, 
Colo. 
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erosion test results that are directly transferable to the prototype structure. 
Mr. Shukry infers in his paper that this was to be expected. 

If prototype erosion was greater than that predicted by the original sec- 
tional model, unsymmetrical operation of the 46 bays could have produced al- 
most any conceivable scour pattern. Degradation of the prototype riverbed 
could also account for lack of agreement between mode! and prototype erosion 
patterns. Clear water released from a reservoir tends to pick up a silt load 
as soon as possible. The riverbed below a dam is sometimes attacked in this 
manner when the bed material is very fine. This action has nothing in com- 
mon with the erosive action encountered downstream from stilling basins. 
Prototype riverbeds have been lowered several feet in a few years as a result 
of degradation. 

Mr. Shukry states that for the Edfina barrage and others that the best loca- 
tion for the sill is in the middle of the floor. This statement is meaningless 
unless the length of the solid floor is defined in terms of the flow depth or in 
some dimensionless number. Sufficient data are not presented to allow the 
reader to compute the probable length of apron required, but it appears from 
Figure 1 that the aprons are longer than those usually required for discharges 
and heads of this magnitude. For example, it appears that a floor extending to 
Section IV or V, Figure la, would provide a sufficiently long paved floor. A 
sill at the end of the apron and riprap downstream from the sill would ordi- 
narily provide ample protection to the riverbed, since the total energy in the 
flow is relatively small. 

A sill is most effective in preventing scour when it is located so that mild 
bottom velocities strike the sill and are directed away from the riverbed ina 
positive and gentle manner. Near-zero velocity then prevails at the bed and 
water interface and optimum conditions for the prevention of scour are 
present. With the sill located in the middle of the apron, a redistribution of 
velocity occurs before the flow contacts the riverbed and possibly before a 
boundary layer is established, resulting in greater scour. Broad sloping sills - 
usually produce weaker ground rollers than steep sills, reducing the circula- 
tion and loss of bed material at the end of the apron. A properly shaped 
dentated sill also helps to produce zero velocity at the interface and a negli- 
gible ground roller. Scores of stilling basins with dentated sills have been 
tested in various laboratories and have been found to be superior to solid 
sills. 

It should also be pointed out that the action described on the Edfina apron 
should not be classed as hydraulic jump action. It appears (sufficient data 
are not given) that the Froude number of the incoming flow is very low, per- 
haps no greater than 2. In a hydraulic jump, the energy loss is only 6 percent 
of the total when the Froude number is 2. If the jump is drowned, the loss 
would be even less. The action on the apron in the Edfina tests would probably 
be better described, therefore, as a roller action where energy loss is related 
to the turbulence developed on the apron. Statements concerning the shape and 
placement of sills should therefore be applied, at best, to the structures 


tested and should not be considered to be generally true for hydraulic jump Es 
aprons. 


ALFRED C. INGERSOLL, ! M. ASCE.—In this paper Mr. Shukry has con- 
tributed greatly to the hudraulic designer’s understanding of the effectiveness 
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of barrage floor sills at various locations in minimizing downstream scour. 
The need for a clarification of the previously existing multifarious and con- 
flicting design criteria is admirably demonstrated by the tests on a movable 
bed model of the old Esna Barrage. The model tests predicted scour be- 

havior that has been largely refuted by experience on full-scale structures. 

As sound as the author’s principal conclusion - that anti-scour sills can 
be investigated by means of fixed-bed model studies - appears to be, it is 
nontheless unfortunate that no case is cited wherein either a movable-bed 
model or a prototype structure behaves clearly according to the fixed-bed 
predictions. It is true that the fixed-bed studies on the Edfina Barrage seem 
to bear out the earlier of the two sets of movable-bed studies, but the more 
recent set of tests, shown in Fig. 2, gives results diametrically opposed to 
the conclusions of the fixed-bed studies, with the mid-floor (or edge) position 
of the sill being the best by one test method and the worst by the other. It 
would seem that this all-important question cannot be truly resolved until the 
behavior of this and similiar prototype structures can be observed. 

In Eq. (6) the author appears to have evolved an expression for the shear 
stress that is independent of the viscosity or Reynolds number, and he there- 
fore proceeds to apply the Froude law to the fixed-bed model studies, com- 
paring one model test qualitatively with another, to predict which prototype 
structure will experience the least scour. The writer feels that there is a 
danger in this, and that it lies in the effect of the Reynolds number upon the 
velocity distribution (especially that downstream) and consequently upon # and 
therefore upon T>. 

Although the analogy is admittedly rough, one may seek a relation between 
8 and TR from information on smooth pipes. Inasmuch as the scale of the 
Edfina Barrage model is 1:66.7, the Reynolds number scale, when tested by 
the Froude law, is 1:(66.7)3/2 or 1:543. As shown by Rouse,(1) the variation 
to be expected in # for this change in Reynolds number is on the order of 
from 1.01 at the prototype scale to 1.03 at the model scale. 

While this magnitude of difference in itself may be negligible, it may be 
observed that it probably applies more to the downstream condition, e.g. 
Station VII, than to the upstream condition at Station I. That is, the velocity 
distribution immediately downstream from the gate is affected most by the 
gate opening and position, the viscous forces having little opportunity to come 
into play in the short distance before Station I. As the flow continues over 
the long floor and the bed-protection works, however, the viscous forces tend 
to establish a more normal velocity distribution, as is shown in Fig. 6. Thus 
when the slightly differing values of 89 are subtracted (as in Eq. 6) from a 
value of By which is the same for both model and prototype, the resulting dif- 
ference in T, could be significant. The net effect of this would be for the 
prototype to develop a higher relative bed shear stress than the model and it 
could conceivably occur that fixed-bed model studies of competing prototype 
designs would lead to an erroneous evaluation on this account. 

The curves of Fig. 7 are especially interesting to the student of open- 
channel hydraulics. The rising surface profiles for tranquil flow over a 
horizontal bed are, of course, impossible for normal channel velocity dis- 
tributions. For the marked maldistributions observed, however, the validity 
of these profiles may be confirmed by an energy analysis applied to the 
reach between Stations I and VII. 
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Considering condition (a), with discharge passed under the gate, the 
velocity distribution curve at Station VII may be analyzed to yield a value 
for the energy correction factor, 


v.? bY (1) 


of 1.32. With a depth of 4.50 meters and average velocity of 2.13 ms. per sec, 
the specific energy, 


(2) 


is 4.81 ms. 

Allowing for a reasonable loss of, say, 5 cm between Stations I and VII, 
the specific energy at Station I must be about 4.86 ms. From Fig. 7(a) the 
depth at I may be scaled at 0.68 x 4.50 or 3.06 ms. The kinetic energy must 
then be 1.80 ms. Inasmuch as the average velocity at this section is 3.14 ms. 
per sec., the velocity head is only 0.501 ms. and clearly this will require a 
value of @ of 1.80/0.501 or 3.59. 

Turning to the velocity distribution diagram for Station I of Fig. 6(a), it is 
possible to compute a value for a by scaling the diagram, provided it be 
assumed that this diagram is representative of the velocity distribution 
across the entire 10.5-meter width of one vent and pier. The value of @ so 
obtained is approximately 2.92, not sufficient to account for the high kinecux 
energy of 1.80 meters. 

The explanation for this discrepancy seems to lie in the lateral distribu- 
tion of velocity as it is badly distorted at the downstream edge of the 2.5- 
meter wide pier, just at Station I. If, then one assumes that the distribution 
diagram shown applies only to the 8-meter width of the channel upstream, 
while the average velocity, V, applies to the full 10.5-meter width down- 
stream, the value of @ becomes 4.85. This is indeed more than enough to 
account for the kinetic energy. The difference between this and the minimum 
requirement for @ of 3.59 isreadily absorbed in considering lateral 
velocity distributions between the two extreme cases postulated. 

This paper is extremely valuable in opening new vistas of research, where- 
in the elusive scour phenomenon can be studied with the fixed-bed model and 
the pitot tube. It is to be hoped that the author will continue these significant 
studies and that he will, whenever possible, relate the predictions so ob- 
tained to observations of scour on prototype structures. 
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and remodelling of the barrages shown in the authors Fig. 1, and in this 
capacity, he ventures to make some constructive remarks about Professor 
Shukry’s interesting and timely paper. 

The authors’ argument is roughly this: since it is impossible to reproduce 
in the model the true ratio of the structural dimensions of the prototype to 
the size of the particles of the granular soil, let us drop this material al- 
together, and carry out the mode) tests within rigid boundaries, deriving our 
conclusions from the observed behaviour of the flowing water. 

This is not a new approach, but it was usually believed to be particularly 
suitable for simpler structures, such for instance, as the Gully Control works 
which were the subject of the experiments of B. T. Morris and D. C. 


; Johnson(1); but whether it is also applicable to works of the size and im- 


portance of the Nile Barrages—on the safety of each of which may easily 
depend the welfare of as much as a million rural population—is a matter re- 
quiring much thought, and consideration of the various aspects of the relevant 
problems. 

Had the author been personally responsible for such designs, he would 
have certainly realised the complexity of these problems and might have 
then, very probably, abstained from veiled complaints against the alleged 
monotony of the adopted design solutions. (Which remark should not be inter- 
preted as adverse criticism, but only as an attempt to clarify the writers 
argument.) 

Let us begin with the Assiut Barrage (d in Fig. 1) the remodelling of which 
was the first among the large Nile Works, to be designed according to model 
tests completed before the drawings were prepared and actual work began. 
The author has envisaged, in his experiments, three alternatives—flow above, 
in between and below the gates—but he has not considered the fourth alterna- 
tive, i.e., flow above the undergate weir; and yet, such a weir formed (accord- 
ing to the model tests referred to) a significant part of the anti-scour-devices 
arrangement, adopted for the remodelling of this work. 

Originally, the object of this weir was a matter of gate design. In fact, the 
new regulation order, for which the Barrage was then being remodelled, 
called for a greater total height of obstruction than the old one, but the per- 
missible height for the new gates was limited by considerations pertinent to 
the housing space available behind the parapet. Hence the weir. 

But when the consulting engineers, Messrs. Coode, Wilson, Mitchell and 
Vaughan Lee, suggested metal trestles above the parapet level and in spite 
of strong opposition the suggestion was adopted (see Fig. 2),(2) the housing 
space problem was ipso facto solved, and the weir was therefore capabie of 
being dispensed with. 

This should have been, apparently, the proper solution from the hydraulic 
standpoint, since it yielded a more regular pattern of flow lines, but a 
general analysis of all the then available evidence, derived from numerous 
model tests with granular material, tended to show the merits of the under- 
gate weir in fighting tail-erosion; and the results in practice, confirmed 
entirely this conclusion, notwithstanding the fact that it had been derived 
from experiments with sand, which the author appears to condemn. 

In the pre-modelling period, the Inspector in charge lived, during the 
flood months, in constant terror of what that blessed Barrage would do next, 
for apart from a permanent scour-hole (“bayarah”), local scour in closer 
vicinity to the apron was not unfrequent, and had to be immediately levelled 
out with rubble-stone dipped from barges, of which more than 130,000 cub.met. 
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had been thrown into the water, since the barrage had been built. Little 
wonder, therefore, that when we saw that after the remodelling, the surface 
of the earthen channel behind the barrage remained almost scourless, we be- 
gan believing in miracles. 

Of course, model tests with sand are not to be interpreted too literally, but 
with a certain experience in the subject, you may very well distinguish in 
examining the information they supply, the basic and substantial, from the 
accidental and inaccurate. In fact, inconsistencies and inaccuracies in such 
experiments when they occur are not so much the result of differences in 
compactness and cohesion and other systematic causes, as the author seems 
to believe, but are chiefly due to the fact that it is difficult (if not impossible) 
to reproduce accurately in a very small size experiment, the true case his- 
tory of the phenomenon. 

The only way of dealing with such essentially accidental inaccuracies is to 
increase the scale of the tests; which means that for important works, such 
as the Nile Barrages, the scale of an experiment with sand, must not be less 
than 1/10. This explains that, immediately after the completion of the Assiut 
Barrage Remodelling, we included in the project of Mohammad Aly Barrages, 
a generously proportioned test-laboratory, in which such large scale experi- 
ments could easily be carried out. It may be unfortunate that, for administra- 
tive reasons, this laboratory was not completed, until after the second world 
war. 

Reverting to the author’s tests, it might be of interest to find out whether 
their method is sensitive enough to perceive the effect of the undergate weir 
on downstream erosion? Additional experiments to clarify this point appear 
to be, certainly, desirable. 

The writer wishes to emphasize that he has no intention whatever of ob- 
jecting against the author’s solution, but he believes that the more usual type 
of model-test with sand, must take precedence on it, and that the two systems 
can then be used in parallel, according to a consistent, self-contained pro- 
gram, aiming at mutual verification of the results. 

The next important work on the Nile, after the Remodelling of Assiut Bar- 
rage was completed, was the construction of the Mohammad Aly Barrages. 
The hydraulic lay-out of this project, for which the writer was particularly 
responsibic, incorporated two main requirements: - 


A) To take advantage of the lesson of Assiut Barrage. 

B) To take, also, into account that the backwater curve of the new barrage 
extended beyond Cairo, a large city with over a million population and 
many valuable old buildings resting on very dubious, if any, foundations. 


Thus, as the yearly flood began to subside and the water table started to 
receed, frequent instances occurred in which such buildings fell, causing loss 
of life and property. It was, therefore, essential not to widen the range of 
water level variations, i.e. reduce the additional heading-up, consequent on 
the construction of the new barrage, to a minimum. 

The success of the solution of this problem deserves, perhaps, being 
placed on record: during the highest flood after construction, the afflux on 
this barrage did not exceed 2 cm (3/4''). To achieve this result the writer 
removed (or smoothed out) all solid surfaces exposed to the direct action of 
the shock of water, or causing separation of flow from solid boundary. 

Thus, though the Remodelling Project of Assiut Barrage was taken as 
basis, the sill of the hypothetical undergate weir was joined with the 
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downstream and upstream edges of the floor, by straight lines; thus yielding 
two sloping surfaces on each side of the highest point under the gates. 

Better to explain this semi-intuitive solution, attention is called to Fig. b(3) 
in which the profiles of the three main Nile Barrages are superposed upon 
one another. In this drawing, the method of the sloping flour of the Mohammad 
Aly Barrages, as used in loco of the weir of Assiut Barrage, is quite obvious, 
and the solution seems to be evident but. . . the writer was, at the time, 
severely criticized for using it, because it was then alleged that sloping floors 
intensified downstream scour, and the designs for Gebel Aulia Dam, on the 
White Nile, which was then in course of construction, were modified accord- 
ingly. 
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The heresy, in fact, arose from a wrong interpretation of another set of 
model experiments with sand, and since subsequent practice fully vindicated 
the writers views on the matter, this problem is another subject on which 
Prof. Shukry’s experiments could throw more welcome, additional light. In 
fact sloping floors seem not to have been considered by him at all. 

It remains to analyse the conditions prevailing at the barrage of 
Nag-Hammadi (see Fig. 1-e). This is a sui generis case, because the floor 
of this barrage is built to three different levels (58.50, 59.50 and 60,50 metres 
above the conventional irrigation datum, which at some time was believed to 
be the mean sea level at Alexandria). 

As regards downstream scour, this floor arrangement is by no means a 
negligible contingency; in fact, as shown by both, the writer’s calculations of 
the probable depth of scour (see Fig. c),(4) and also by actual soundings taken 
in the channel of the river, downstream of the apron, the depth of scour de- 
creases gradually with the discharge in the part with the highest floor level, 
but in the deepest portion of the barrage, it begins to increase again as the 
river drops to its summer levels. 

This is another aspect of the scour problem for the Nile Barrages, which 
might be covered by Professor Shukry’s valuable tests. 
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But, the main point on which further experiments of the same type appear 
to be necessary, is that the scour-hole itself will basically affect the flow 
pattern, and therefore, to be able to draw significant conclusions from Prof. 
Shukry’s tests, we must possess, in addition to the results obtained for a 
level channel-bed, full information obtained in the same manner for scour- 
holes of various depths; viz., using available empirical information on the 
profiles of the most probable scour-holes, the solid bed, which in these ex- 
periments replaces the earthen channel, must be shaped accordingly, and 
then, tests carried out and diagrams prepared as for a level channel for vari- 
ous assumed depths of the scour-hole. This will then supply exhaustive in- 
formation on the matter and may, possibly, yield a convincing result. 

In fact, it must be realised that scouring action may be started at any time 
and for any accidental reasons, and we must be sure that if this occurs, the 
hole will be filled in naturally, or will not attain a dangerous depth, before we 
have had time enough to fill it in with rubble. 
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Broadly speaking, the Nile Barrages must not be visualised as granite 
castles built on a rocky foundation, but rather, as islands made of friable 
material, which are held in precarious equilbrium, afload on an immense 
sea of almost oily, semi-liquid misture of water and very fine sediment. 
Granted that this description is an obvious exaggeration, from the designers 
standpoint the exaggeration is in the right direction. 

Every precaution must, consequently, be taken to ensure stability of 
granular material, and departures from established practice must be care- 
fully considered. 

The writer, therefore, concludes that it is too early to dispense with the 
standard type of model-tests with sand altogether, replacing them by 
Professor Shukry’s interesting method; but this last method must be used to 
interpret, and sometimes explain, the older type of such experiments. 


RE FERENCES 


. “Hydraulic Design of Drop Structures for Gully Control,” Trans. ASCE., 
Vol. 108, pp. 887 and foll. 


. From S. Leliavsky, “Irrigation and Hydraulic Design,” Vol. II, 1957, Ed. 
Chapman & Hall Ltd., London, Fig. 82. 


3. Leliavsky, loc. cit., Fig. 123. 


. Fig. c is computed from Figs. 170 to 172 in the same book referred to in 
Reference 2, but in Vol. I, pages 218-220. 


hg 


ASCE 1456-25 


STILLING BASIN EXPERIENCES OF THE CORPS OF ENGINEERS# 


Discussion by Alvin J. Peterka 


ALVIN J. PETERKA, | M. ASCE.—The paper by Mr. Berryhill is an excel- 
lent summary of “how things have worked out.” With a few exceptions, it ap- 
pears that design procedures based on hydraulic model tests have provided 
relatively trouble-free operation. 

The writer wishes to raise a question, however, concerning the proper 
depth of tail water, Dg, for a stilling basin containing baffle piers. It has 
been common practice in some offices to reduce the conjugate depth, Dg, by 
10 to 20 percent because of the restraining effect of the baffle piers. Certain 
economies result; the basin floor may be set at a higher elevation reducing 
excavation costs, and the training walls may be correspondingly reduced in 
height. In some cases these savings result in a lowered basin cost. In other 
cases, however, the savings, in terms of other costs, are minor. 

In either case, the writer believes that full conjugate depth or more (about 


10 percent more) should be used when setting the floor elevation, for the fol- 
lowing reasons: 


1. The safety of the entire structure may hinge on whether the baffle piers 
remain fully effective after they have been subjected to the wear of abrasion 
and/or the erosion caused by cavitation. Rounded baffle piers are considera- 
bly less effective than square-edged piers and only nominal damage may 
render the baffle piers ineffective. With the baffle piers ineffective, the apron 
is too high and the tail water too low to produce the proper jump action. 
Dangerous riverbed erosion may result. Difficulties typical of a high apron 
are illustrated by Lock and Dam No. 1, Missouri River and Bull Shoals Dam. 


2. Due to degradation of the river channel after the dam is constructed, 
the tail water elevation becomes progressively lower over a period of years, 
reducing the factor of safety against jump sweepout and lowering the effi- 
ciency of the jump. On the other hand, an apron designed for conjugate depth 
will operate properly even if the baffle piers are not replaced. 


3. Experience has shown that actual tail water elevations are often lower 
than the computed values, used for design or model testing of the basin. This 
may have been true for Kanopolis Dam, which was found by the author to be 
operating at 0.7 Dg. Unless the tail water curve is known to be correct, 
extra tail water depth should be provided in the design. 


4. In the temporal sense, the actual tail water depth lags the depth shown 
on the tail water curve during periods of increasing discharge. On some 
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rivers, it may take several hours for the tail water to stabilize after the dis- 
charge is set. When it is necessary to open spillway gates rapidly, the basin 


must be capable of performing satisfactorily with tail water elevations below 
the stable values. 


5. A slightly drowned jump, in the writer’s opinion, produces a better all- 
around energy dissipator than one at conjugate depth. The slight loss in effi- 
ciency (energy dissipation) is more than offset by the smoother downstream 
water surface conditions and the reduction in splash and surging. A slight 
excess of tail water also moves the toe of the jump against the sloping face of 


the spillway (if one is present), stabilizing the jump and producing smoother 
overall operation. 


On the Center Hill, Conchos, Norfork and Fort Gibson Dams, Mr. Berryhill 
describes damage caused by bed meterial or other debris circulating on the 
apron. The writer has also witnessed damage of this type on other structures 
and feels that abrasion damage is more wide spread than generally believed. 
Every effort should, therefore, be made, particularly by those engaged in 
hydraulic model testing to eliminate “dead” areas in basins where debris may 

~ become trapped. Insofar as possible, basins or buckets should be made self- 
cleaning and the bed material just downstream from the end sill should be 
stabilized or reduced in elevation to prevent bed material from being swept 
into the basin during periods of necessary unsymmetrical operation. It has 
been common practice on model tested basins to use an end sill which 
creates a ground roller and moves bed material upstream against the end 
sill to prevent undermining of the apron. The writer believes that strong 
ground rollers should be avoided by providing end sills with flatter slopes 
or dentated or notched sills which tend to neutralize the upstream motion. 
With continued motion of the bed material against the end sill, some bed 
material is bound to find its way into the basin and become the source of 
abrasion damage. 

In discussing the damage at Lucky Peak Dam, Mr. Berryhill indicates that 
perhaps the damage was similar to that found in laboratory tests on concrete 
specimens subjected to the impingement of high velocity jets. In the case of 
the laboratory tests, the damaged area was in the form of an annular ring 
adjacent to the impingement area. The writer believes that at Lucky Peak 
Dam the damage resulted from an entirely different cause since the damage 
was located some distance from the impinging jet. It is very possible that a 
vortex formed at the upstream corners of each slide gate, similar to those 
which form behind radial gates on a spillway. At small gate openings, the 
low pressure core of the vortex was directed toward the floor, as mentioned 
by Mr. Berryhill, and deflected inwardly by the gate slots. Cavitation then 
occurred in the vortex core near the floor. Downstream, a high pressure 
area collapsed the cavities and damage resulted as shown in the photographs. 
The extreme damage further downstream was probably a progressive result 
of the damage which began closer to the gate. The eroded areas in Figure 13a 

indicate that the low pressure area originated behind the gate. The vortex 
theory should be investigated at the first opportunity in a hydraulic model. 
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DIVIDED FLOW THROUGH A DIVERGENT INLET CONDUIT* 


Stavros Tsakonas,** A.M. ASCE 
(Proc. Paper 1492) 


SYNOPSIS 


This investigation is concerned with the determination of the characteris- 
tics of divided flow phenomena. Two models have been analytically and ex- 
perimentally investigated. One is a conduit with a constant cross-section 
having an opening perpendicular to its length. The other is a conduit witha 
divergent inlet having a side opening at the end of the divergent. The analy- 
sis is developed by means of the free stream line theory, assuming two- 
dimensional irrotational flow of an incompressible fluid. The experiments 
were conducted under approximately two-dimensional flow conditions. The 
theory, together with the experiments, has shown quantitatively the influence 
of the velocity distribution in the inlet on the characteristics of the divided 
flow phenomena. 


INTRODUCTION 


Efficient diversion of a portion of a main flow for auxiliary purposes is 
vitally important in designing large scale hydraulic systems. Experiments 
have shown that a unique angle exists for a branch conduit which will keep 
losses at a minimum. A conduit with a slit or hole in its side produces a 
free jet stream at a certain natural angle to the main flow. This angle cor- 
responds to the angle of least loss for the case of a branch conduit. 

The objective of this investigation was to determine the parameters affect- 
ing this natural angle and other flow characteristics and to study their influ- 
ence on the jet diverted through the slit opening. 

Earlier experimental work(1,2) and analytic studies(3,4,5) have shown that 
the natural angle and other flow characteristics such as the contraction coef- 
ficient, width of jet etc. depend on two factors. These are (a) the geometry of 
the junction and (b) the ratio (r) of the discharge diverted into the branch to 


Note: Discussion open until May 1, 1958. Paper 1492 is part of the copyrighted 
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the total discharge in the main conduit beyond the junction point. Previous 
scale model tests conducted by the Bureau of Reclamation(6) have shown 
large discrepancies with the results obtained in small laboratory tests and 
also have differed with prototype measurements. 

At the suggestion of Professor A. Craya, the present study was initiated to 
determine the influence of the velocity distribution in the upstream sections 
of the main conduit on the natural angle and other characteristics of the flow 
diverted through the slit. No known experimental or analytical work had 
considered this effect in the past. 

The problem of a general velocity distribution in a main conduit preceding 
the junction point is neither theoretically tractable nor experimentally feasi- 
ble. The mathematical difficulties of the problem are increased greatly be- 
cause the use of the complex variable technique is not possible in a rotational 
flow. In addition the problem of altering the velocity distribution at the inlet 
of the conduit by laboratory means is difficult.(7) However, the use of a 
divergent inlet conduit offers a means of altering the velocity distribution. 
The divergent inlet conduit is advantageous in that a small angle of divergence 
produces considerable changes in the velocity distribution, as shown by 
Nikuradge.(8) Ordinarily such velocity distribution changes would require an 
inlet length of 20 to 50 feet before the actual experimental section. The in- 
fluence of the angle of divergence itself, i.e., in the case of irrotational flow 
on the jet characteristics can be determined analytically. This influence has 
been shown to be small for small angles of divergence. 

The theoretical part of this investigation dealt with two models. One is a 
conduit with a constant cross-section and a side opening, while the other is a 
conduit with a divergent inlet having a side opening at the end of the sloping 
part. Both cases were investigated as problems in two-dimensional irrota- 
tional flow of an ideal incompressible fluid by means of the free streamline 
theory. The experimental part of this investigation dealt with the same two 
cases of straight and divergent inlet conduits under approximate two- 
dimensional flow conditions. 

Comparison of theoretical and experimental results obtained for the 
straight conduit model served to determine the degree of the validity of the 
free streamline theory together with the previously stated assumptions. The 
second model with the divergent inlet conduit was used to determine the in- 
fluence of the non-uniform velocity distributions. 

If the deviation between the theory and the experiments in the case of the 
divergent channel is greater than the differences observed between the theory 
and the experiments in the case of the straight conduit, this deviation may be 
attributed to the influence of the velocity distribution. This was assumed 
since all other parameters (ratio of discharges, geometry of the opening and 
property of the fluid) were kept constant and only the velocity distribution 
was changed. 


In summary, the present investigation was concerned with: 


1. the theoretical determination of the natural angle of the jet and other flow 
characteristics for divided flow through a conduit with a constant cross- 


section as well as for a divergent inlet conduit, based on the free streamline 
theory of irrotational flow; 


2. the experimental determination of the natural angle of the jet in both 
models. The angles used were: = 0°, @= 2045'48", = 6925" 43”; 
3. the magnitude of the influence of the velocity distribution at the inlet on 
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the characteristics of the jet diverted through the side opening. This also 
required measurement of velocity distributions in the experiments. 

The experimental results showed agreement with the free streamline 
theory when the velocity distribution at the inlet was uniform as in the case 
of the straight conduit. However, when the velocity distribution was non- 
uniform as in the divergent conduit, the experiments showed a deviation from 
the theory, but they remained similar in their general aspects. 

The experiments thus showed that the velocity distribution in a inlet con- 
duit has a definite influence on the flow pattern of a divided flow field. 


Analysis 


The derivation of a theoretical solution for the flow characteristics of the 
problem of a divided flow through a divergent inlet conduit will be considered 
in the following analysis. 

The inflow through the divergent channel A", A' , N and P was divided into 
a flow through a slit (PQ) and a straight conduit (NBQ) as shown in Figure 1. 
By successive transformations the physical plane was transformed to the 
hodograph plane (Figure 2), then to the logarithmic plane (Figure 3), and 
finally, by the Schwartz - Christoffel transformation, to the { -plane (Figure 
4) where physical properties of the flow pattern could be identified easily. 

By considering the physical plane with a source of strength (q) at the point 
of intersection of the sides (AP and AN) of the divergent channel, the result- 
ing idealized flow is such that the field of interest is bounded by streamlines 
along which the velocity is either constant in direction or constant in magni- 
tude. The line AO (O is a stagnation point) is the boundary which divides the 
flow field into an upper part discharged through the side opening (PQ) and a 
lower part conveyed through the main conduit (NBQ). The point N is another 
stagnation point of a different nature since it does not split the flow field. 


The hodograph plane represents the complex conjugate velocities of the 
flow field, hence 


W=-u-iv, (1-1) 


where u and v are the velocity components of the point (x,y) in the z-plane. 
The logrithmic plane is related to the W-plane by 


U = log W. (1-2) 


The origin is chosen at the point P. This requires the magnitude of the 


velocity at P to be equal to unity as well as to the points Q and C. This in- 
volves no loss of generality. 


The ¢ -plane is the result of the Schwartz-Christoffel transformation: 


wat “We, (1-3) 


where P is a complex constant and the points Q, P and N have been chosen at 
=0,+~, -1, respectively. 
Equation (1-3) may be written in a partial fraction form as 


dU 1 A B C 
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a = 
a-Ww)(a-v) 


(1-5) 


Vv) (vy -W) 


= H-w)(B -w) 


a-W)(v-W) 


Using t = /€ and spliting the terms of Equation (1-4) into complex conjugate 
forms and then integrating them yields: 


iva ivy tidy ivo VE tivw 


The condition Up = 0 at t = t+ ~ has been used to eliminate the constant of 
integration. 

The constants of the transformations have been determined by relating the 
¢ -plane to the U-plane, since on traversing the boundary in the U-plane the 
change of arguments of U is known. 

Using Equation (1-6) in the following form 


i Va 


+i 


ivw 


in order to describe the singular points O, A and N leads to 


Cc 
i Vo 


=1, A $ B 


=k (1-7) 


However, it must be kept in mind that the constants A, B and C of the 


transformation are purely imaginary since P can be readily shown to be 
imaginary.(10) 


Combination of Equations (1-5) and (1-7) yields 


BL 1)(a-w) 

Cy. -w)(a- 1) 
k B a-W l= T) (o) 


where vy = 1 is due to the transformation. 
The substitution of “a =aand /W = w in Equation (1-8) leads to: 


(1-8) 


+ a? + apwe - (a* +a +1) 


(1-9) 
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Equation (1-10) is quadratic in w and yields the solution 


2 {k(a-1)-w] 


(1-11) 
where the negative sign is the proper one satisfying the relation uw >a, im- 
posed by the transformation. (See Figure 4) 
The complex potential function (Z) in the ¢ -plane for a source of strength 
q at £ = -a@and sinks of strength q' and q" at £ = -Band ¢ =¥ respectively 


is given by 


2=- Flor +2 roe(t+p) +2 (1-12) 


The complex conjugate velocity and complex potential are related by 


dZ dz d (1-13) 


Setting W= O at the stagnation point (f = - w) leads to 


a-W Y¥+W 


Al 


where r = = with r being the ratio of the discharge diverted into the branch 


to the discharge into the main conduit before the junction point. 

Substitution of Equation (1-7) for the values of A, B and C in Equation (1-6) 
and use of Equation (1-2) gives the velocity at any point in the field by the fol- 
lowing: 


k k 
-iva + i i vw 
As a result of this Equation (1-13) leads to 


where 
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Equation (1-16), which relates the £ -plane to the physical plane through a 
differential form, becomes 


C+w 1-17 


after using Equation (1-14). At this stage of the development the parameters 
a@,W, and ¥ are still undetermined. An additional relation between these 


parameters is evident from the geometry of the configuration shown in 
Figure 1. 


which, by means of Equation (1-17), yields 


(a-s)(a+7) tane 


1-18 


Combining the expression of the velocity at B, found from Equation (1-15) for 
= -B, with 


leads to 


k k 
h! l-r ( Ja JB 


This relates the width of the diverted jet to the height of the conduit. 
Use of Equation (1-19) together with (1-14), changes (1-18) into the form: 


(1-20) 


K 


k 


This integral relation will determine the parameter Y for givena,B,w. 
Using the intermediate relation 


k k 


the width of the slit is given by Equation (1-21) after the integration of 
Equation (1-17) from £ = Oto =+~. 
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The position of the stagnation point of the flow pattern is given by: 


The contraction coefficient “c” of diverted flow is defined as 


(1-23) 
: 


where ¢ and h" are given by Equation (1-21) and Equation (1-19), respectively. 


The angle of the diverted jet at < is given by setting = ¥ in Equation 
(1-15) 


Since W. = e and Iv, | = 1, 
then 


-1 -1 - 
tan V@ i2k tan —— - i2 tan 
Y 


Taking the principal values of the multivalued functions, the angle of the 
impinging jet is given by 


a 2 a -1 
@ = 2k tan —= = 2ktan + 2 tan (1-24) 
JY /F ° 


The analysis developed thus far permits determination of all characteris- 
tics of the flow pattern providing the ratio of discharge and the geometry of 
the configuration are known. 

For simplicity the following method of computation was adopted in this 
investigation. 

Assuming values of @ > 1 and 0< W <1, mw was determined by Equation 
(1-11), 8 by Equation (1-9) and ¥ by the integral relation of Equation (1-20). 
The flow characteristics and the corresponding boundaries were obtained as 
follows: 


6 by Equation (1-24), ¢/h by Equation (1-21) , C by Equation (1-23) , 
h''/h' by Equation (1-19), S by Equation (1-22) and r by Equation (1-14). 
This computation procedure has the advantage of eliminating all but one of the 
integral relations between the parameters. The integral relation which must 
be used is given by Equation (1-20). 


| | 
| 
| 
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2 - Solution for a Straight Conduit 


Using the expression evolved for the general case of the divergent inlet 


channel, the flow characteristics for the straight conduit (k = :. 0) will be 
developed. 
For the case of k = 0 Equation (1-11) gives 


or 


The former case requires “= a@ = 8 =w = 1 and is a special case which will 
be discussed later. The latter case (“@ =a) requires # = 1, thus allowing 
Equation (1-15) to be reduced to 


y « ivw (2-1) 
Jt + idw 


When ¢ = -B= - 1 and ¢ = - a, Equation (2-1) becomes, respectively, 


1 -VW 


3 


Wee [Ya] * 


Using the latter expressions, h and h' can be determined as 


be 


A Va - 


5 1- 
Since h = h', Equation (2-2) takes the form 


(2-3) 


after using Equation (1-14). Once the values for « and w are chosen, the 
value of ¥ may be obtained from Equation (2-3). 
The angle of diverted jet is determined by Equation (1-24) as 


2 tan 


y 
and 
and 
(2-4) 
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The corresponding ratio of discharges is given by Equation (1-14) as 


_@-1 Y+w (2-5) 
a-w 


rs. 


The width of the slit, upon integration of Equation (1-14)(10) leads to: 


1 


"(1+ )(a+y) 


+w)( ¥ +1)loga=( Y=w)( log ¥+2 a= Wa] 
(2-6) 


The corresponding contraction coefficient given by 


Cs n( @-1){ ¥+w) 


(@+W)( %+1)loga-( Yaw)(@ -L)log Y+2 JW SY (a-i)n 


(2-7) 


and the ratio of the width of the jet to the height of conduit is determined by 
Equation (1-19) as 
h" _ (@- 1)(Y%+w) 


TT 5 e 2-8 


Two limited cases (r = 1 anda@ > ~ with O< r< 1) were investigated and 
are described as follows: 


1) Case r = 1, 


then W = 1 and Y= ka -( V@~ 1)" (2-9) 
(V@+1)* 


Ifa >1, 


then lim 
a 


= 0.611 : (2-10) 


1 
+n Uel 


90° lim +0 and C= 
ei. 


This limiting case corresponds to a flow through an orifice on an infinitely 
large container. 


— co 


then lim, = 3, @= 60 , lim 


independently of the values of ¥ and w, 


and lim 
a 


- 
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This case (r = 1,@> ~) corresponds to the flow pattern through an orifice, 
the dimensions of which may be considered infinitely large in respect to the 
height of the conduit. The jet is diverted through an angle of 60° without any 
contraction. 


2) Casea» nwo and 0<r<il, 


then parameters Y and w are determined, respectively, by 


1-J/o )? 


——ele-r and 

1+ Jw 
and the characteristics of the flow field are determined by Equations (2-4), 
(2-6), (2-7) and (2-8). Noticing that w decreases to zero faster than 7, it fol- 
lows from the characteristic relations of the flow pattern that 6 = 0, 
a = 0 and C = 0 independently of the values (/ h', when w = 0. 

The case @>-« corresponds to an orifice with an infinitely large opening 
compared to the height of the conduit. The jet is diverted without any contrac- 
tion and the quantity h"/h' takes a maximum value for each ratio of discharge. 
The jet emerges through the free surface in the vicinity of the downstream 
edge of the slit which extends an infinite distance from the upstream leading 
edge. 

Results for a straight conduit solution as well as for its limiting cases, 
are plotted in Figures 5,6 and 7. Although it does not seem possible to show 
the identity of the expressions developed in the present analysis with the cor- 
responding expressions developed by Craya(3) and McNown,\°) the plotted re- 
sults show complete agreement with Craya’s results and good agreement with 
McNown’s results. The small discrepancy in the latter case is probably due 
to the fact that, in reducing the results to the same basis, use was made of 
values obtained from the graphs plotted by McNown. 


3 - Solutions for a Small Angle of Divergence 


After the successive transformations of £ = - 0 and /O = t the integral 
relation (1-20), which determines the parameter 7, leads to 


(28) 2 k k 2 
2 t= va ate VE 
(t°-a) (t+ y) t-1 tand(a=- +Y) 


(3-1) 


Ja 


Upon integration and substitution for 8B by Equation (1-9) this takes the 
form 


(0-4) 


f(a, k). (3-2) 


Having these two functions f(a, Y , w , k) and d(a@, ¥,w, k) which are 
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analytic functions of k in some common region (R) with respect to all 
parameters (@, »,w) and the independent variable k, we actually postulate 
that only on some surface S(@, ¥ , w), which belongs to the regions of analy- 
ticity, does the equality (3-2) hold. The functions may be expressed as power 
Series in k in the region of common analyticity: 


; 2 


Then on the surface (S), the equality (3-2) implies that coefficients of like 
powers of k are equal 


AC as,%,W) a,7,w) (a) 
A(a,Y,w)=Bla, Y,W) (b) (3-3) 
) (c) 


Equations a, b, c, etc. define several surfaces ina, Y,w space. If any two of 
the surfaces differ they may or may not intersect. If they do intersect, Equa- 
tion (3-2) would apply only to restricted values of one of the parameters, 
namely those satisfying the equation of the intersection. However, this con- 
tradicts the postulate based on physical grounds that no restriction exists on 
any of the parameters in the region P of the analyticity. 

If the surfaces do not intersect, no values of the parameter exist which are 
common to both surfaces. Hence, there is no surface for which (3-2) holds. 
This also contradicts the postulate that S(a, Y , W) must exist. 

Consequently, the remaining equations in (3-3) must be identically satis- 
fied or equivalent to (3-3a) which is the only one that specifies S(a, ¥, w). 

Equation (3-1), expressed in k powers with attention being concentrated on 
only the first term,(10) leads to: 


2 
(a= (3-4) 


(VB 


B,Y¥,W) k0,(a, By )*... 
(@- p)(p +7) 


On the other hand were Equation (1-9) written as 


2 2 
p(a‘tatl) -a -yw 


Equation (3-4) then obtains its final form: 


la 2 
VW 


k 

2 
1 
k 


» Y,W )+kd,( 2, ¥ @ 


(3-6) 
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Therefore by the argument stated above, the surface S(a, ”, w) to which 
Equation (3-1) holds is given by 


(3-7) 
lt Y 


This is the same relationship as Equation (2-3) found in the analysis for the 
case of k = 0, without using the integral relation (1-20). 


It can be readily found that the integral relation should include the case of 
k = 0 by means of Equation (3-5), since 


k 


and the right hand side of Equation (3-1) does not vanish. 

The discussion in the preceding paragraphs eliminates the need for evalu- 
ating the integral in Equation (3-1), which was attempted both by a closed- 
form method and by numerical computations. The closed-form method was 
not tractable, while numerical computations were arduous, difficult and did 
not maintain the required accuracy. 

By means of Equations (1-9), (1-11) and 3-7) the values of 8, u and ¥ can 
be determined for assumed values of @> 1 and 0< W <1. The flow charac- 
teristics then can be computed in close-form except for the ratio of the slit 


opening to the height of conduit ((/h' ) and the contraction coefficient (c). 
The investigation of some limiting cases leads to: 


1) whenu> ox, 
Equations (1-11), (1-9) and 3-7) yield, respectively, 


a- (2-2) 
f 
b 
+V@ +1 -w 
and 


[2Va@ +( Va+1) Vo] Vo (c) 
Va+1 


For the case and assumed values of 0<w <1, the values ofa, 
and ¥ can be computed by Equation (3-8 a, b, c) afterwhich the flow charac- 
teristics 6, r, h"/h' can be determined by the general expressions given in 
section 1. This case physically corresponds to a conduit with an extremely 
large opening in its side. The width of the diverted jet is smaller than that 
attained in the case of straight conduit. Furthermore, it is noticeable that, 


as the size of the slit opening increases and tends to infinity, the width of the 
jet diverted through the opening approaches zero. 


2) When >~ and a= 
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Equations (3-8) leads to: w > oB> 0.3333 and Y >0 ( ¥ does not approach 
zero as quickly as w) 
For this case 9 = 0, h"/h' = O and r= 0. The entire discharge is conveyed 
through the main conduit. 


3) Whena > 1, 

Equation (1-10) leads to “ = @ = 1, independent of the value of w. Then 
Equation (1-9) gives B the value of 1 for any value of w, while Equation (3-7) 
or (3-8c) yields Y= /V@ . However, Equation (1-14) has meaning only if 
Ww = 1, which leads to: 


r=l and h"/nt — 0 


This case corresponds to the physical case of an infinitely large container 
with a small side opening. The jet is diverted through an angle of 90° and 
the angle of divergence has no influence on the jet’s characteristics. 

Up to this stage of the analysis no approximation has been made. The flow 
characteristics r, 9, h''/h' have been determined in closed-form expressions 


and the results are plotted in Figure 7 for the cases of = a and @= a5 


In order to reach a conclusion on the influence of the angle of divergence 
in the case of an irrotational solution, the study of the physical characteris- 
tics (/h' must be considered. At this stage it is assumed that k is small, 
i.e., the angle of divergence (¢) is smaller than 5 = 6926'. Equation (1-21) 


can be written in an expanded form, where only the two first terms are taken 
into account(10) as 


=o pat + i2ck )al +... , 
° 


(3-9) 
= Cl, + i2ckI, 
where 
) an FF n Te 
and 


k k 
n( Va - /F 1+ 


The first term (I,) of Equation (3-9) is the value of //h' for k = 0 and its 
coefficient (c) reduces to 


n( Q+Y 


—— (if 
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The integrand of the second term Ig is a combination of the integrand of I; 


and the function @(¢) which is bounded. Its maximum value is attained as 
=Va@. Asa result 


1 1 4 
max = tan - tan va 


Therefore, the following inequality holds: 
fore} oo 

= 


where 


k _ 


The actual value of the second term of Equation (3-9) is much smaller than 
I'y. It can be shown by calculation that I'p = (1.5 = 3%) 1,- Since I, <1, the 
use of the second term (Ig) in Equation (3-9) with the first term (Ij) will alter 
the values of (/h' to the second or third decimal point. 

Therefore, for the case of a small angle of divergence and fora <9 which 
covers the range of practical applications, it can be assumed that the lines of 
(/h' remain almost invariable as k varies from 0 to = 

Using this fact the theoretical results obtained from hr Versus @ plane 
(Figure 7) can be transferred to 6 versus r plane (Figures 10 and 11) which 
has been used for comparison with experimental work. 


4 - Experimental Investigation 


The main features of the equipment and instrumentations as well as the 
experimental procedure are discussed briefly in this section. The entire ex- 
perimental work was conducted in the Fluid Laboratory at Columbia Univer- 
sity. The fluid used in the experimental work was water. The open channel 
in which the experiments were conducted was 21 feet long, 15 inches wide and 
37 inches deep. It is part of a circuit formed by a reservoir, pump, constant 
head tank, entrance section, two exit tanks with measuring weirs and a re- 
turning channel leading to the reservoir or to a volumetric measuring tank. 

The actual test section was made by inserting two carriers (Figure 8) and 
a partially movable bottom into the channel so that it was divided into four 
parts: (1) an entrance chamber, (2) a main conduit, (3) an upper chamber 
and (4) the remainder of the channel. The lifting mechanism under the par- 
tially movable bottom provided a maximum angle of convergence and diver- 
gence of some seven degrees in the inlet section of the main conduit. 

The measuring equipment consisted of a pitot-tube with a traversing 
mechanism, standard point gages, a protractor, a micromanometer, a set of 


tubing to feed sufficient quantity of dye into the test section and a pair of tri- 
angular weirs to measure the two discharges. 
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The pitot-tube consisted of a streamlined stem and three removable tips 
to permit reading of the differential pressures at three different stations of 
the inlet conduit. 

The protractor consisted of a quadrant which was divided into one degree 
divisions, a carriage which bore the pointer and a micrometer head used for 
small adjustments. The protractor was used for measuring the angle of the 
diverted jet. This was done by centering the pointer in the jet flow and by 
making it parallel to the jet. 

The main problem in the experimental work was to regulate and control 
the discharge of water passing through the slit and the main conduit, to read 
the angle of the diverting jets and to determine the fluid’s velocity at different 
sections in the inlet conduit. 

Regulation of the discharge was achieved by calibrated weirs of known dis- 
charge coefficients. Main readings were taken of the two discharges, the 
angle of the side jet, and the velocity distribution in the upstream inlet sec- 
tion. In order to read the angle of the jet, a sufficient quantity of dye was in- 
serted to make the jet visible. A picture also was taken of each test. An- 
other important feature of this investigation was the velocity distributions 
used for comparing the results of straight and divergent inlet conduits. Dif- 
ferential pressure between the pitot-tube and hydrostatic wall tap was mea- 


sured by means of the pitot-tube and micromanometer. The velocity distri- 
butions was obtained by 


‘ 
Vo = 10.40A p ft.“/sece* 


where A p = differential pressure in inches of water 
and P = density in slug / ft.3, 
Velocity distributions were taken at three different stations, 2 7/16, 6 3/8 

and 12 5/16 inches from the leading edge of the opening. The results were 
plotted in a non-dimensional form y/d versus u/Um (see as an example 
Figure 9) where 

y = distance from the central line of the conduit 

d = 1/2 of the corresponding height of the conduit 

u = velocity at the point under consideration 

Um = Maximum velocity at the section under consideration 


Table I gives values of hj, dj, 1;, and the corresponding area (Aj) for the 
various cases under consideration. 


Experimental work was carried out for a straight conduit as well as fora 


divergent inlet channel with angles of divergence ? = ri = 2°45' 48" and 


d = a = 6925' 40". Several slit openings and various discharge ratios were 


used for each geometry. 

In the straight conduit case the velocity distribution was almost uniform 
except at the boundary later. When an angle of divergence was used, the de- 
viation of the velocity distribution from uniform became more and more pro- 
nounced as the angle was increased. In most cases the velocity distribution 
was nearly symmetric, the exceptions being cases in which separation was 
imminent. 


The plot of the experimental results on charts representing the theoretical 
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Table I 


Angle of 
divergence 


Scheme Station 1 Station 3 


Station 2 
= 12 5/16" 


= & 3/8" 


=1.905" 
Ay 1.905 


0.9525" 


> = 1.905" h, = 1.905" 


d, = 0.9525" ad, = 0.9525" 


2 


| A, = 198hft. 


= 1.9852" 
hy 1.9852 


1.5918" 1.3001" 


= 2° 518" = 0.8926" 


0.1859fte- 


0.796" d, = 0.6506" 


2 


A, = 0.1658ft.* | A. = 


= 
1.6303" 


1.1867" h, 0.5177" 


= 6°25143" 


Q 


0.593)" d. = 0.2588" 


> 
u 


0.0539ft. 


solution shows a close agreement in the case of a straight conduit with uni- 
form velocity distribution (Figure 5), but depicts a definite deviation in the 
case of a divergent conduit with non-uniform velocity distribution. As the 
velocity distribution becomes more and more non-uniform the deviation of 


the theoretical and experimental points becomes more and more pronounced. 
(Figures 10 and 11). 


CONCLUSIONS 


The free-streamline theory is the only one which offers an analytic solu- 
tion to the divided flow phenomena, but the stumbling-block for the generali- 
zation of its principles was believed to be attributed to the velocity distribu- 
tion in the inlet. The main goal of this investigation was to determine 
qualitatively the influence of the velocity distribution on the solution obtained 
by the free-streamline theory. 

Both the experimental work and the analysis established the fact that the 
upstream velocity distribution is another essential factor for the determina- 
tion of the divided flow pattern by the free streamline theory. 

This conclusion has been established by comparison of the theoretical 
solution, obtained by the free streamline, with corresponding experiments for 
several different angles of divergence, which are equivalent to considering 
several different velocity distributions. 

As long as the upstream distribution remains essentially uniform, the free 
streamline solution is verified fully by the experiments. However, when the 
velocity distribution is not uniform, the characteristics of the diverted jet 


| 
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shows a systematic difference, but remain similar in general aspects. 

An angle of divergence of @= 2°45' 48", which alters the velocity distribu- 
tion slightly from uniformity, lowers the angle of diverting jet by 1 1/2 to 4 
degrees. An angle of divergence of ¢= 6925" 43", which alters the velocity 
distribution considerably, lowers the angle of the jet inclination by 2 1/2 to 
6 degrees. 

The similarity in results found between analysis and experiments suggests 
an attempt to incorporate the velocity factor to the principle of the free 
streamline theory. This may lead to a generalization of the free streamline 
theory for rotational flow, which at the present is relatively intractable from 
a theoretical stand-point. 

The analysis used in this investigation leads to a close-form solution for 
the flow characteristics in a straight conduit case as well as in the case of a 


divergent inlet conduit with a small angle of divergence (<5). However, 


for any angle of divergence £.<¢<e close-form expressions are obtained for 


the natural angle, the ratio of discharges and the ratio of the width of the jet 
to the height of the conduit. 

Figures 5 through 10 show the variation of the flow characteristics for 
different slit openings and discharge ratios. These variations are discussed 
in Reference 11 as a function of the pressure variation. 

The analysis developed herein, based on the free streamline theory and 
applied to a small angle of divergence of the inlet conduit, is expected to be 
useful for design purposes or/and as a guide for the interpretation of experi- 
mental data for similar geometrical configurations. 
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